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Realistyczna grafika 3D 


Wstęp 


Niniejsza książka jest pierwszą pozycją na polskim rynku wydawniczym, 
poświęconą w całości trójwymiarowej grafice komputerowej. Zawiera podstawy 
teoretyczne grafiki 3D i przykłady jej praktycznego wykorzystania. 


Jeszcze do niedawna przeciętni użytkownicy komputerów osobistych nie mogli 
tworzyć realistycznych obrazów trójwymiarowych, takich jak np. na okładce 
tej książki. Maszyny będące w ich posiadaniu nie miały odpowiedniej mocy 
obliczeniowej ani możliwości graficznych. 


Powszechna dostępność coraz szybszych komputerów sprawia, że rośnie grono 
miłośników zaawansowanej grafiki komputerowej zainteresowanych tajnika- 
mi jej tworzenia. To skłoniło nas do napisania niniejszej książki. 


Książka przeznaczona jest dla osób, które opanowały podstawy języka C. Ze 
względu na znaczną ilość przykładów oraz przystępny język mogą korzystać z 
niej również ci, którzy stykają się z tematem po raz pierwszy. 


Po przestudiowaniu jej treści oraz po wykonaniu przykładów Czytelnik może 
pokusić się o samodzielne zaprogramowanie dowolnych zadań grafiki trójwy- 
miarowej, takich jak tworzenie brył, ich transformacje oraz cieniowanie powie- 
rzchni oświetlanych różnokolorowymi lampami. Ułatwi mu zadanie biblioteka 
procedur znajdująca się na dołączonej dyskietce. Procedury zostały napisane 
w standardowym języku C, co umożliwia uruchamianie ich na różnych typach 
komputerów. 


Na niej również znajduje się napisany przez nas program RAY3D.EXE (pra- 
cujący pod nadzorem środowiska MS-Windows) umożliwiający tworzenie trój- 
wymiarowych scen oraz ich realistyczne przedstawianie. Rysunek na okładce 
został wygenerowany właśnie przez ten program. 


Autorzy 
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ROZDZIAŁ 1 


Matematyka trójwymiarowa 
— podstawowe pojęcia i wzory 


W tym rozdziale pragniemy przedstawić podstawowe pojęcia związane z geo- 
metrią trójwymiarową. Rozdział jest przeznaczony głównie dla mniej zaawan- 
sowanych Czytelników. Pozostałych także zachęcamy do jego przeczytania w 
celu odświeżenia wiadomości ze szkoły. Mamy nadzieję, że zawarte w nim 
informacje ułatwią im lekturę następnych rozdziałów. 


Podstawowymi pojęciami w geometrii analitycznej w przestrzeni są punkt 
i układ współrzędnych. Do opisu obiektów trójwymiarowych będziemy używali 
prostokątnego kartezjańskiego układu współrzędnych (rys. 1.1). W zależności 
od wzajemnych orientacji osi układ może być lewostronny (rys. 1.la) lub 
prawostronny (rys. 1.1b). Istnieje prosta metoda rozstrzygania, jakiego typu 
jest dany układ. Wystarczy wyobrazić sobie, że patrzymy na układ "od strony” 
płaszczyzny X-Y (czyli płaszczyzny rozpiętej na osiach X, Y). Wtedy w obrębie 
naszego pola widzenia oś X leży poziomo i jest skierowana w prawo, a oś Y 
pionowo do góry. Jeżeli oś Z jest skierowana do nas, to układ jest prawostronny. 
W przeciwnym razie układ jest lewostronny. 


Położenie punktu w układzie współrzędnych definiuje się poprzez przyporząd- 
kowanie mu trójki liczb. Liczby te nazywamy współrzędnymi punktu (rys 1.2a). 
Wyrażają one względne położenie danego punktu względem początku układu. 
Punkt przecięcia osi X, Y, Z nazywamy środkiem lub początkiem układu 
współrzędnych. Ten punkt ma współrzędne x=0, y=0, z=0, co w skrócie zapisu- 
jemy (0,0,0). 


Na rysunku 1.2b przedstawiono punkty Pl o współrzędnych (xl,yl,z1l)i P2o0 
współrzędnych (x2,y2,z2). Tworzą one wektor P o współrzędnych [x2-xl,y2- 
yl,z2-z1]. Prosimy zauważyć, że współrzędne punktów zamykamy w nawia- 
sach okrągłych, a wektorów w kwadratowych. Taką konwencją stosujemy w 
całej książce. Długość wektora P, czyli odległość punktu P2 od Pl, wyznaczamy 
ze wzoru: 
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Rys. 1.1 


P(x„y ,Z) 
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P.(GY2/Z2) 


P,(XYZ,) 


Rys 1.2b 
IPl=V(x,-x,)"+(va-y,)7+(2,—2,)7 


W niektórych wzorach będziemy dla skrócenia zapisu używać iloczynów ska- 
larnych i wektorowych par wektorów. Iloczyn skalarny dwu niezerowych 
wektorów A, B wyrażamy wzorem: 


AOB=|A | IBlcosa=xAxB +yAYyB + ZA ZB 


gdzie: 

IAI — oznacza długość wektora A, 
[BI — oznacza długość wektora B, 
a — kąt między wektorami A, B. 


Z powyższego wynika, że wynikiem iloczynu skalarnego dwu wektorów jest 
liczba. Jeżeli wektory są prostopadłe lub jeden z nich ma długość O, to iloczyn 
skalarny wynosi 0. Aby znaleźć kąt między dwoma niezerowymi wektorami, 
należy obliczyć ich iloczyn skalarny i podstawić do wzoru: 


AO0OB XĄAKp TYAĄJB * ZĄZB 
' 
IATIBI VxZtyąĄ+tZĄVXĘE+tygt+ZĘ 


Wynikiem iloczynu wektorowego pary wektorów A, B jest wektor W (patrz 
rys. 1.4): 
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Rys. 1.3 


Rys. 1.4 
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Wektor W jest prostopadły do płaszczyzny rozpiętej na wektorach A, B. Jego 
długość jest równa polu równoboku rozpiętego na tych wektorach. Zwrot 
wektora W jest taki, że trójka A,B,W ma tę samą orientację co przestrzeń, w 
której definiujemy iloczyn wektorowy. Iloczyn wektorowy wektorów równole- 
głych jest wektorem zerowym. A teraz wzór: 


W=AXB=[yĄZzp-yYgZa,XgZy-—X,Zp,XAYg-XpyA/] 


|[WI=IAIIBlIsina 


UWAGA: Powyższy wzór jest prawidłowy wyłącznie w przestrzeni trójwymia- 
rowej z lewostronnym układem współrzędnych. 


Możemy teraz przejść do opisu własności bardziej złożonych obiektów, jakimi 
są prosta i płaszczyzna. Matematycy tworzą różne postaci równań prostej lub 
płaszczyzny. My przytoczymy wyłącznie te, które są najczęściej stosowane w 
grafice komputerowej. 


Równanie płaszczyzny przechodzącej przez punkt Po(xo,yo,zo) i prostopadłej 
do wektora N o współrzędnych [A,B,C]: 


A(x-x,)+B(y-y,)+C(z-z,)=0, 
(A?+B?+C?*>0) 


Rysunek 1.5 ilustruje opisaną sytuację. Równanie to możemy zapisać używając 
iloczynu skalarnego: 


We*(P-P,) 


gdzie P(x,y,z) — oznacza punkt należący do płaszczyzny, a P-Po — wektor na 
płaszczyźnie. Interpretujemy to w następujący sposób: płaszczyzna prostopad- 
ła do wektora N, przechodząca przez punkt Po, jest zbiorem takich punktów 
przestrzeni, dla których iloczyn skalarny wektorów N i P-Po jest równy zero 
(N i P-Po są wektorami prostopadłymi). Wektor prostopadły do danej płasz- 
czyzny będziemy dalej nazywali wektorem normalnym do płaszczyzny. 


Równanie płaszczyzny można zapisać także w postaci ogólnej: 
Ax+By+Cz+D=0, (A”+B?+C?*>0) 


i w postaci normalnej: 
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PO(x0,y0,z0) 


Rys. 1.5 


NOP[x,y,z]-p=0 
N =[cosa ,cosf ,cosy] 


gdzie: 

N — wektor o długości równej 1, prostopadły do płaszczyzny, mający 
zwrot w stronę tej półprzestrzeni, w której nie znajduje się po- 
czątek układu współrzędnych (jeśli p=0, to wektor N może mieć 
dowolny zwrot); 

jo — odległość płaszczyzny od początku układu współrzędnych. 


Prostą, podobnie jak płaszczyznę, możemy opisać kilkoma równaniami. Mając 
dany jeden punkt i kierunek możemy jednoznacznie opisać prostą. Jak wynika 
zrysunku 1.6, współrzędne dowolnego punktu P na prostej wyznaczamy znając 
współrzędne wektora V i punktu Po. 


Równanie prostej przytoczymy w postaci parametrycznej: 
P(t)=P,+t* Vczyli 
X=Xp+tt*ea 
y=y,+t*b 
Z=Zy+t*c 
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X,Y,Z  — współrzędne punktu na prostej, 

Xo,Yo,Zo— współrzędne punktu Po, 

a,b,c — współrzędne wektora V, 

t — dowolna liczba rzeczywista (parametr); im większa, 
tym punkt P leży dalej od Po. 


A "POGOA020 


P(t)=P,+tV 


Rys. 1.6 


Jeżeli znamy współrzędne dwu punktów P1(X1,Y1,21), P2(X2,Y2,22) należących 
do prostej, to jej równanie możemy zapisać w innej postaci: 


lub 


P(t)=P,+t*(P,-P,), czyli 
X=X, +t*(X—X,) 
Y=Jy+t*(Y2—11) 
Z=Z,+t*(Zę—Z,) 
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Pierwszy z wzorów nazywamy równaniem kanonicznym, a drugi parametrycz- 
nym. 


Na rysunku 1.7 przedstawiono interpretację geometryczną ostatniego równa- 
nia. 


Rys. 1.7 


Zauważmy, że dobierając odpowiednio wartość parametru t, otrzymujemy 
punkty leżące między P1 i P2 , bądź poza nimi. Dla t równego O dostajemy 
współrzędne punktu Pq, a dla t = 1 punktu P2. 


Ostatnim z zagadnień teoretycznych będzie wyznaczanie przecięcia prostej z 
płaszczyzną. 


Prosta o równaniu 
P(t)=P,+t* V 
przecina płaszczyznę 
NOP-p=0 


w pewnym punkcie P (patrz rys. 1.8). Aby znaleźć współrzędne tego punktu, 
wystarczy wyznaczyć wartość parametru ti podstawić ją do równania prostej. 
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W tym celu punkt P podstawiamy do równania płaszczyzny (punkt przecięcia 
spełnia zarówno równanie prostej, jak i płaszczyzny) i otrzymujemy: 


No(P,+t* V-p)=NoP0+t*NoV-p=0 


P.(X:YorZe) 


P(t)=P,+tV 
Rys. 1.8 
a stąd: 
p—NoP, 
- nov 


Jeśli prosta i płaszczyzna są równoległe (patrz rys. 1.9), to wyrażenie w 
mianowniku ma wartość O (wynika to z własności iloczynu skalarnego — wektor 
N normalny do płaszczyzny jest prostopadły do wektora kierunkowego pro- 
stej V). 


Otrzymana zależność jest szczególnie często wykorzystywana w algorytmach 
znajdowania widocznych części obiektów trójwymiarowych i algorytmach typu 
ray tracing. 
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P.„YZ,) 


Rys. 1.9 


Z 
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ROZDZIAŁ 2 


Komputerowy opis 
obiektów trójwymiarowych 


W rozdziale 1 zapoznaliśmy Czytelnika z podstawowymi pojęciami i wzorami 
używanymi w grafice komputerowej. Opisane zostały wyłącznie najbardziej 
podstawowe obiekty, czyli punkt, prosta i płaszczyzna. Oczywiście otaczający 
nas świat jest znacznie bardziej skomplikowany. Mnogość kształtów i form jest 
zadziwiająca. Wszystkie te kształty mają pewne wspólne cechy, które możemy 
łatwo wyodrębnić i opisać. Takimi cechami mogą być na przykład wymiary lub 
kolor powierzchni. Obiekty mogą być kanciaste — prostokątny stół, pudełko 
zapałek, dom, bądź zaokrąglone — piłka, jajko, głowa ludzka. Jednak większość 
otaczających nas przedmiotów posiada obie te cechy, czego najlepszym przy- 
kładem może być zwykła szklanka. To zmusza nas do uwzględniania w modelu 
obiektu wielu różnych właściwości, komplikując go. 


Aby tworzyć realistyczne obrazy za pomocą metod grafiki komputerowej, 
musimy na samym początku stworzyć matematyczny opis interesujących nas 
obiektów. Dla niektórych obiektów istnieją proste wzory opisujące ich powie- 
rzchnie. W poradnikach matematycznych można znaleźć równania opisujące 
powierzchnię kuli, elipsoidy, walca, paraboloidy itp. Problemy zaczynają się 
pojawiać dopiero przy opisie nieregularnych powierzchni. Jak na przykład 
opisać powierzchnię liścia lub kształt chmury? Dopiero dynamicznie rozwija- 
jąca się w ostatnich lat geometria fraktalna pozwala na badanie tych nieregu- 
larnych kształtów. Stosowany tam aparat matematyczny nie jest tak oczywisty 
jak szkolna geometria. Ponadto realizacja odpowiednich algorytmów na kom- 
puterach jest czasochłonna — wymaga więc znacznych mocy obliczeniowych, 
nieosiągalnych dla komputerów osobistych. 


W tym momencie należy zwrócić uwagę na to, że naszym celem jest nie tylko 
znalezienie matematycznego opisu obiektu, ale także wykorzystanie go do 
późniejszej wizualizacji na ekranie komputera. To nakłada na nas pewne 
dodatkowe ograniczenia. Opis powinien być dokładny, ale także pozwalający 
na zastosowanie efektywnych algorytmów do przekształceń i prezentacji obie- 
któw. Podstawowym kryterium wyboru sposobu opisu obiektów 3D będzie to, 
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czy dana ich reprezentacja pozwala na maksymalnie szybkie znajdowanie nie 
zasłoniętych powierzchni. Jak uczy doświadczenie, najbardziej czasochłonną 
czynnością, jaką wykonuje algorytm wizualizacji obiektów trójwymiarowych, 
jest właśnie badanie, w jakich punktach jedna powierzchnia zasłania inną. 


Najprostszym i najczęściej stosowanym sposobem opisu brył trójwymiarowych 
jest reprezentacja szkieletowa. Do tej metody będziemy się odwoływać w 
dalszej części książki. W niniejszym rozdziale zaprezentujemy podstawowe 
elementy składające się na opis szkieletowy bryły i ich realizację na przykładzie 
odpowiednich struktur i funkcji języka C. Zaprezentujemy też zestaw funkcji 
przeznaczonych do tworzenia prostopadłościanów, kul, walców, stożków i brył 
obrotowych. 


2.1. Reprezentacja szkieletowa 


Zgodnie z zasadą mówiącą, że jeden obraz niesie więcej informacji niż tysiąc 
słów, wyjaśnianie koncepcji szkieletowego opisu brył zaczniemy od prezentacji 
rysunku pionka (rys. 2.1). 


Zauważmy, że powierzchnia pionka nie jest gładka. Jest ona zbiorem trójkąt- 
nych ścian. Każda ściana jest ograniczona przez trzy krawędzie. Każdą ścianę 
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wyznaczają trzy wierzchołki. Krawędzie łączą sąsiednie wierzchołki. To jest 
właśnie szkieletowy opis bryły. Za pomocą trzech podstawowych elementów: 
wierzchołka, krawędzi i ściany zbudowaliśmy dość skomplikowany geometry- 
cznie obiekt. Gdybyśmy próbowali znaleźć funkcje matematyczne opisujące 
taką powierzchnię, okazałoby się, że są niezwykle skomplikowane, a używanie 
ich do wizualizacji wymagało by bardzo szybkiego komputera. 


W jaki sposób można pokazać taką bryłę na monitorze komputerowym? Naj- 
prostszą metodą jest wyznaczenie dla każdego trójkąta jego jasności i naryso- 
wanie go. W zależności od położenia źródła światła otrzymamy inny rozkład 
oświetlenia ścian. Przykładowy wygląd pionka przedstawia rys. 2.2. 


Z pewnością Czytelnik zauważył, że otrzymaliśmy "kanciasty" pionek. Jest to 
pewna ułomność szkieletowego opisu bryły. Można ją jednak niemal całkowicie 
skompensować, dzieląc bryłę na więcej trójkątnych ścian i stosując przy 
wyznaczaniu jasności każdego punktu powierzchni odpowiednie wzory inter- 
polacyjne. Wzory zostaną zaprezentowane w podrozdziale 5.5, a efekty ich 
użycia prezentuje rys. 2.3. 
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Rys. 2.3 


2.2. Elementy reprezentacji szkieletowej 


Jak już wspomnieliśmy, reprezentacja szkieletowa bryły jest realizowana za 
pomocą trzech elementów: wierzchołka, krawędzi i ściany. W pamięci kompu- 
tera zbiory odpowiednich elementów możemy przedstawić jako tablice lub listy. 
Struktura listowa jest nieco bardziej skomplikowana w realizacji niż tablica, 
jednak pozwala na swobodniejsze operowanie elementami. Podstawową zaletą 
listy jest możliwość dodawania i usuwania dowolnej, nie znanej na początku 
liczby elementów — liczba elementów na liście może się dynamicznie zmieniać. 
Tej zalety nie posiadają tablice — tworząc tablicę musimy z góry znać liczbę jej 
elementów. Po zapełnieniu tablicy dodanie kolejnego elementu nie jest możli- 
we. My wybraliśmy listy jako podstawowy sposób organizacji danych w przy- 
kładach prezentowanych w książce i w programie demonstracyjnym. 


W następnych punktach opiszemy struktury realizujące odpowiednie elementy 
i funkcje służące do ich łączenia. Wszystkie przykłady są napisane w języku C. 
Dla zwiększenia ich czytelności (szczególnie, gdy chodzi o przedstawianie 
połączeń między strukturami) będziemy używać także bloczków z opisanymi 
polami. Każdy bloczek reprezentuje jedną strukturę. Na niektórych rysunkach 
pojawią się pola z trzema kółkami. Oznaczają one składowe struktur nieistotne 
w danym przykładzie. Definicje wszystkich opisywanych przez nas struktur, 
deklaracje stałych i nagłówki funkcji znajdują się na dyskietce w pliku ray.h. 
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Kod źródłowy wszystkich opisanych funkcji, a także funkcji nieopisanych, ale 
używanych w programie demonstracyjnym, znajduje się na dyskietce. 


2.2.1. Wierzchołek 


W trójwymiarowej przestrzeni podstawowym pojęciem jest punkt. W struktu- 
rach opisujących bryły, a także w innych przypadkach struktura nazwana 
przez nas VERTEX (ang. wierzchołek) pełni tę samą rolę. VERTEX służy do 
przechowywania informacji o współrzędnych punktu. Oto struktura w języ- 
ku C: 


typedef struct vert 
/* struktura opisująca wierzchołek */ 
( 


int info; 


double 
double 
double 


normx; 
normy; 
normz; 


/7 


współrzędne globalne wierzchołka 


double globalx; 
double globaly; 
double globalz; 


double 
double 
double 


double 
double 
int 
int 
int 


int y2 


int 
int 


localx; 
localy; 
localz; 


drawx; 
drawy; 


xl; 


y3; 


/7/ 


// 


współrzędne lokalne wierzchołka 


współrzędne do rysowania na ekranie 


struct vert huge *next; 
/* następny wierzchołek na liście */ 


JVERTEX; 
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VERTEX 


o 


Globalfx,y,z] 


globalx, globaly, — współrzędne wierzchołka wyrażone w ukła- 
globalz dzie globalnym. 

globalx, globaly, — współrzędne wierzchołka wyrażone w ukła- 
globalz dzie lokalnym. 


Układ globalny możemy utożsamiać z układem współrzędnych opisanym na 
początku rozdziału 1. Pojęcia układów lokalnego i globalnego zostaną szczegó- 
łowo wyjaśnione w następnym rozdziale. 


drawx, drawy — współrzędne rzutu perspektywicznego wie- 
rzchołka. Także pojęcie rzutu perspektywi- 
cznego zostanie wyjaśnione w następnym 
rozdziale. 


next — wskaźnik do struktury opisującej następny 
wierzchołek na liście wierzchołków. Rys. 
2.4 przedstawia sposób łączenia struktur 
VERTEX w listę jednokierunkową. 


Pozostałe pola są w tym momencie nieistotne i zostaną opisane później. 
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Local 
Global 


Rys. 2.4 


2.2.2. Krawędź 


Krawędzie opisują struktury EDGE (ang. krawędź). Każda krawędź jest ogra- 
niczona dwoma wierzchołkami. Krawędzie, podobnie jak wierzchołki, są połą- 
czone w listy za pomocą pól next. Informacja zapisana w tych strukturach jest 
używana wyłącznie podczas rysowania brył w postaci szkieletowej. 


EDGE 


Vert2 


typedef struct ed /* opis krawędzi */ 


t 
VUERTEX huge *vertl; // wskaźnik do końca krawędzi 
VERTEX huge *vert2; // wskaźnik do drugiego 
końca krawędzi 
int info; 
struct ed huge *next; 
// nastepna krawędź 


) EDGE; 


*vertl, *vert2 — wskaźniki do struktur opisujących wierz- 
chołki będące końcami krawędzi. 
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2.2.3. Ściana 


Struktura FACE opisuje ścianę, Jak pamiętamy, ściany są trójkątami. Wynika 
stąd, że do opisu danej ściany wystarczy podać trzy wierzchołki bryły należące 


do niej. 


typedef struct fc /* opis ściany */ 


( 


VERTEX huge 
VERTEX huge 


VERTEX huge 


// wskaźnik do pierwszego z trzech pkt. 
*vertl; 

// wskaźnik do drugiego z trzech pkt. 
*vert2; 

// wskaźnik do trzeciego z trzech pkt. 
*vert3; 


double vnormxn; 

double vnormyn; 

double vnormzn; 

double VNOoLmc; 

int vnorminface; 

int info; 

int phong; // = ON oznacza wygładzenie ściany 


FEATURE  *feature; 


// cechy ściany 


struct fc huge *nextAF; 


// lista aktywnych ścian 


struct fc huge *next; 


) FACE; 


*next 


vnormxn, 
vnormyn, 
vnormzn 


// następna ściana 


— wskaźnik do struktury opisującej następną 
ścianę na liście ścian. 


— współrzędne wektora normalnego (prosto- 
padłego) do ściany. Wektor normalny ma 
długość 1 i jest obliczany na podstawie ilo- 
czynu wektorowego wektorów utworzonych 
z dwu krawędzi danej ściany. Ilustruje to 
rys 2.0. 


Pozostałe pola tej struktury zostaną opisane później. 


Na rysunkach 2.6 i 2.7 przedstawiono przykładową ścianę i strukturę połączeń 
pomiędzy opisującymi ją elementami FACE, EDGE, VERTEX. 
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2.2.4. Bryła 


Struktura OBJECT zawiera szkieletowy opis danej bryły. Składa się ona z list 
wierzchołków, krawędzi i ścian, których elementy łączą się ze sobą w odpo- 
wiedni sposób. Odpowiednie pola struktury OBJECT wskazują na początki 
tych list. 


Rys. 2.5 
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VERTEX 4 


OBJECT 


Name 


Number Faces 


Ust Verices 


typedef struct obj // struktura opisująca bryłę 
i 
char *name; // nazwa nadana przez użytkownika 
int number; // unikalny numer obiektu 
int type; 
int n_vertices; // liczba wierzchołków 
int n_edges; // liczba krawędzi 
int n_ faces; // liczba ścian 


FEATURE *l features; 
VERTEX huge *l vertices; 
// wsk. do pierwszego el. listy wierzchołków 

EDGE huge *l edges; // — or — ;, ,-krawędzi 
FACE huge *]| faces; // — or = , ,>Ścian 
AXES *axeS; 

// osie ukł. lok., do którego należy ten obiekt 
JOBJECT; 


*name — tekst zawierający nazwę bryły, 

*| _vertices — wskazuje na początek listy wierzchołków, 
*] _edges — wskazuje na początek listy krawędzi, 

*] faces — wskazuje na początek listy ścian, 

*axes — wskaźnik do struktury AXES opisującej 


układ lokalny, w którym jest umieszczona 
bryła. Struktura AXES zostanie opisana w 
następnym rozdziale. 
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Rys. 2.8 


Na rys. 2.8 przedstawiono przykładowy obiekt składający się z dwóch ścian. 
Wszystkie elementy tego obiektu zostały ponumerowane (EDGE +1 oznacza 
krawędź numer 1, EDGE +2 — krawędź numer 2, VERTEX +41 — wierzchołek 
numer 1 itd.). Rysunek 2.9 ilustruje, jakie struktury danych opisują obiekt z 
rys. 2.8. Dla zachowania czytelności nie umieszczono na nim połączeń między 
ścianami i krawędziami a wierzchołkami. 


2.3. Funkcje 


W podrozdziale 2.2 opisaliśmy realizację opisu szkieletowego brył jako zestaw 
czterech struktur utworzonych w języku C. Teraz przedstawimy funkcje, za 
pomocą których można tworzyć wspomniane struktury, łączyć je w listy 
i usuwać z tych list. Wszystkie funkcje znajdują się na dyskietce. Mamy 
nadzieję, że Czytelnik wybaczy nam nieco pobieżny opis niektórych funkcji. 
Wydaje się nam, że zbyt szczegółowe opisywanie sposobu działania programu 
daje skutek przeciwny do zamierzonego — opis staje się nieczytelny, a odbiorca 
zaczyna gubić się w natłoku informacji. Lepszy efekt daje w takim przypadku 
samodzielne przanalizowanie interesującego nas fragmentu kodu. 


EE, 
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2.3.1. Tworzenie struktur 


Pod pojęciem tworzenia struktur rozumiemy oczywiście przydzielanie pamięci 
operacyjnej potrzebnej do przechowywania informacji zawartej w danej stru- 
kturze. Standardową funkcją realizującą to zadanie jest CALLOC (n, dłu- 
gość_bloku), zwracająca wskaźnik do bloku pamięci o rozmiarze n*długość_blo- 
ku bajtów. Za jej pomocą można zerezerwować pamięć dla n struktur, każda 
składająca się z długość_bloku bajtów. W naszych przykładach będziemy się 
posługiwać funkcją FARCALLOC, pozwalającą na przydzielenie bloku pamięci 
o rozmiarze większym niż 64 KB, czego nie można zrealizować za pomocą 
CALLOC. Taki sposób zarządzania pamięcią jest stosowany w komputerach 
typu PC. Czytelnik posiadający dostęp do innego typu komputera będzie mógł 
w poniższych przykładach zastąpić FARCALLOC zwykłą funkcją CALLOC. 
Należyjednak pamiętać o usunięciu modyfikatorów huge. Funkcja new_vertex, 
przedstawiona poniżej, rezerwuje pamięć dla n struktur VERTEX: 


VERTEX huge *new vertex( int n ) 
t 

return((VERTEX huge * )farcalloc(n,sizeof (VERTEX) )); 
) 


Podobną budowę mają funkcje rezerwujące pamięć dla struktur opisujących 
krawędzie i ściany: 


EDGE huge *new edge( int n ) 

return((EDGE huge * )farcalloc(n,sizeof (EDGE) )); 
RACE huge *new face( int n ) 

return((FACE huge * )farcalloc(n,sizeof(FACE))); 
OBJECH *new _object( void ) 

return((OBJECT * )farcalloc(l,sizeof (OBJECT))); 
) 


Ostatnia funkcja tworzy opisaną wcześniej strukturę OBJECT, zawierającą 
opis szkieletowy bryły. Teraz możemy napisać funkcje dołączające struktury 
opisujące wierzchołki, krawędzie i ściany do bryły OBJECT. 


a 
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łinclude "ray.h" 
VERTEX huge *add vertex to object( VERTEX huge *vert, 
OBJECT *object , double x , double y , double z ) 
l 
if( ( vert=add vertex( vert, Xx, y, Z ) ) == NULL ) 
return NULL; 


vert->next=object->l vertices; 
object->l vertices=vert; 
return( vert ); 


) 


Pierwszym argumentem tej funkcji jest wskaźnik *vert. Jeśli nie jest równy 
NULL, to znaczy, że wcześniej został zarezerwowany blok pamięci i funkcja 
add_vertex nie musi rezerwować pamięci. Struktura VERTEX, opisująca wie- 
rzchołki, jest dołączana na początku listy wierzchołków bryły. Funkcja zwraca 
wskaźnik do struktury albo NULL, gdy zabrakło pamięci. 


Funkcja add_vertex, przedstawiona niżej, inicjuje współrzędne wierzchołka 
zgodnie z argumentami x, y, z. Jeżeli wskaźnik vert = NULL, to jest wywoły- 
wana funkcja new_vertex, przydzielająca blok pamięci dla jednej struktury. 
Jeżeli vert <> NULL, to znaczy, że wcześniej w programie została zarezerwo- 
wana pamięć dla tej strukturyi jedynym zadaniem tej funkcji jest zainicjowanie 
pól współrzędnych. Jak się później okaże, takie rozwiązanie jest bardzo korzy- 
stne. Na przykład tworząc opis sześcianu, z góry wiemy, że potrzeba sześciu 
wierzchołków. Można więc zarezerwować jeden duży blok pamięci, w którym 
zmieści się sześć struktur opisujących wierzchołki, zamiast sześciokrotnie 
przydzielać pamięć dla kolejnych struktur. Nasze rozwiązanie ma dwie zalety: 
po pierwsze oszczędza czas (mniej wywołań funkcji CALLOC), po drugie 
oszczędza pamięć. 


Ostatnie stwierdzenie wymaga pewnych wyjaśnień. Otóż funkcja CALLOC jest 
najczęściej realizowana w taki sposób, że dzieli pamięć na bloki 16-bajtowe. 
Najmniejszą porcją danych, jaką może zwrócić taka funkcja, jest właśnie 16 
bajtów. Jeżeli tworzymy strukturę o długości 10 bajtów, to wiąże się to z 
przydzieleniem bloku o długości 16 bajtów. Tracimy w takim przypadku 6 
bajtów. Wynika stąd, że w celu uzyskania pamięci dla stu struktur 10-bajto- 
wych bardziej efektywne jest przydzielenie jednego bloku 1000-bajtowego, niż 
stu bloków 10-bajtowych. W pierwszym przypadku potrzeba 1000 bajtów, a w 
drugim 100*16 = 1600 bajtów. 
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include "ray.h" 
VERTEX huge *add vertex( VERTEX huge *vert , 
double x , double y , double z ) 
( 
if( vert == NULL ) 
t 
if( ( vert=new vertex( 1 ) ) == NULL ) 
return NULL; 
) 
vert->globalx=Xx; 
vert->globaly=y; 
vert->globalz=z; 
/* inicjowanie współrzędnych lokalnych */ 
/* na początku są takie same, jak globalne */ 
vert->localx=xX; 
vert->localy=y; 
vert->localz=z; 
/* inicjowanie wsp. do rysowania */ 
/* na początku są takie same, jak globalne */ 
vert->drawx=xX; 
vert->drawy=y; 
return( vert ); 


Nie opisujemy sposobu działania funkcji dołączających struktury EDGE 
i FACE, opisujących krawędzie i ściany do OBJECT, ponieważ działają na tej 
samej zasadzie co add_vertex_to_object. 


łinclude "ray.h" 
EDGE huge *add edge _to object( EDGE huge *edge, OBJECT 
*object , VERTEX huge *vertl , VERTEX huge *vert2 ) 
( 

if( edge == NULL ) 

( 

if( ( edge=new edge( 1 ) ) == NULL ) 
return NULL; 


edge->vertl=vertl; 
edge->vert2=vert2; 
edge->next=object->l _edges; 
object->l edges=edge; 


NZZNNNNNNNNZZNNZZEZZZZZZZZNENEZNZZ ZZ 
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return( edge ); 


) 


FACE huge *add_face_ to object( FACE huge *face , OBJECT 
*object , VERTEX huge *vertl , VERTEX huge *vert2 , 
VERTEX huge *vert3 , FEATURE *fe , int phong ) 
( 

if( face == NULL ) 

t 

if( ( face=new _face( l ) ) == NULL ) 
return NULL; 
> 


face->vertl=vertl; 
face->vert2=vert2; 
face->vert3=vert3; 
face->feature=fe; 
face->phong=phong; 
// nie inicjujemy wsp. wektora normalnego do ściany 
face->next=object->l_faces; 
object->l_faces=face; 
return( face ); 


) 


OBJECT 
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Rys. 2.10 
OBJECT 


DEYOM 


c) 


Rys. 2.10 
Parametry fe i phong zostaną opisane w następnych rozdziałach. 


Dołączenie nowego wierzchołka bryły następuje na początku listy wierzchoł- 
ków. Tę operację ilustruje rys. 2.10. 
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Podobnie jest realizowane dołączanie krawędzi i ścian w funkcjach add_ed- 
ge_to_object, add_face_to_object. 


2.3.2. Usuwanie struktur 


Funkcjami uzupełniającymi dla funkcji tworzących struktury są funkcje słu- 
żące do zwalniania pamięci zajmowanej przez te struktury. Przedstawimy 
teraz procedury, za pomocą których można kasować zarówno konkretne stru- 
ktury, jak i całe listy. 


Kasowanie jednej struktury VERTEX lub jednego bloku struktur VERTEX 
(jeżeli pamięć dla wszystkich struktur została przydzielona za pomocą jedno- 
krotnego wywołania funkcji new_vertex(n)) jest realizowane za pomocą funkcji 
delete_vertex. Natomiast kasowanie listy wierzchołków (przy założeniu, że 
tworząc kolejne struktury używaliśmy funkcji new_vertex(1)) wykonuje proce- 
dura delete_vertices. 


łinclude "ray.h" 

// zwalnia pamięć zajętą przez strukturę opisującą 

// wierzchołek vert, jeśli vert nie wskazuje na NULL. 

// zwraca wskaźnik do następnika skasowanego wierzchołka. 
VERTEX huge *delete_vertex( VERTEX huge *vert ) 


l 
VERTEX huge *v next; 


if( vert == NULL ) 
return NULL; 
v_next=vert->next; 
farfree( vert ); 
return( v_next ); 


) 


// zwalnia pamięć zajętą przez listę wierzchołków 
// wskazywaną przez vert 
void delete _vertices( VERTEX huge *vert ) 


t 
if( vert != NULL ) 
l 
while( vert->next != NULL ) 
vert=vert->next; 
delete_vertex( vert ); 
> 
> 
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Pary funkcji delete_edge, delete_edges i delete_fece, delete_faces, a także 
funkcja delete_feature mają podobną budowę. 


Zwalnianie pamięci przez struktury opisujące bryłę wykonuje funkcja dele- 
te_object. Najpierw usuwane są listy wierzchołków, krawędzi i ścian, a potem 
sama struktura OBJECT. 


łinclude "ray.h" 
void delete _object( OBJECT *object ) 


t 
if( object == NULL ) 
return; 
delete_faces( object->l faces ); 
delete _edges( object->l edges ); 
delete _vertices( object->l vertices ); 
delete features ( object->l features ); 
if( object->name != NULL ) 
farfree( object->name ); 
farfree( object ); 
> 


2.3.3. Tworzenie brył 


Wiemy już w jaki sposób można tworzyć opis brył w pamięci komputera, 
składając je z trzech podstawowych elementów, tj. wierzchołków, krawędzi 
i ścian. Zaprezentujemy teraz kilka przykładowych procedur generujących 
automatycznie wybrane bryły. Procedury te otrzymują parametry określające 
wymiary obiektów, nazwę nadaną przez użytkownika itp. i na tej podstawie 
tworzą odpowiednie struktury. Na koniec pokażemy, w jaki sposób można 
stworzyć rysunek obrazujący daną bryłę, wykorzystując informację zawartą w 
tych strukturach. Wszystkie opisane niżej funkcje zostały włączone do progra- 
mu demonstracyjnego, tak więc Czytelnik może się zapoznać z efektem ich 
działania. Jednak przed uruchomieniem tego programu zalecamy przeczytanie 
instrukcji obsługi zamieszczonej w rozdziale 7. Kod źródłowy zarówno opisa- 
nych, jak i pozostałych procedur znajduje się w pliku obiekty.c. 


Na wstępie przedstawiamy najkrótszą funkcję, służącą do tworzenia prosto- 
kątnego wycinka płaszczyzny. Funkcja add_plane tworzy strukturę opisującą 
prostokąt równoległy do płaszczyzny X-Z. Nadaje mu numer : number. Krawę- 
dzie mają długość a i b (patrz rys. 2.11). Jeżeli podamy argument name = 
NULL, to nazwa obiektu jest tworzona automatycznie. Funkcja zwraca 
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wskaźnik do utworzonej struktury lub NULL, gdy brak pamięci. Ponieważ 
wszystkie obiekty tworzone przez opisywane funkcje posiadają wyłącznie 
trójkątne ściany, prostokąt musi składać się z dwu ścian. 


Pewnych wyjaśnień wymagają parametry feature i phong. Feature jest 
wskaźnikiem do struktury FEATURE, opisującej cechy powierzchni bryły. 
Definicja struktury FEATURE znajduje się w pliku ray.h. Dodatkowo funkcja 
add_plane wywołuje funkcję new_feature. Jej zadaniem jest utworzenie stru- 
ktury FEATURE. Kod źródłowy tej funkcji wygląda następująco: 


include "ray.h" 

// Utworzenie struktury opisującej parametry 

// powierzchni ściany 

FEATURE *new feature (double ktR ,double ktG, double ktB, 
double krR,double krG, double krB, 
double kz, int m ) 


FEATURE *fe; 


if( ( fe=( FEATURE * )farcalloc(l , sizeof(FEATURE) )) 
== NULL ) 

return NULL; 
fe->number=0; 
fe->ktR=ktR; 
fe->ktG=ktG; 
fe->ktB=ktB; 
fe->krR=krR; 
fe->krwmwaG=krG; 
fe->krB=krB; 
fe->kz=kz; 
fe->m=n; 
fe->next=NULL; 
return fe; 


) 


Dołączenie struktury otrzymanej w wyniku działania tej funkcji do obiektu 
realizuje funkcja add_feature_to_object. 


łinclude "ray.h" 
void add_feature_to object( OBJECT *object, FEATURE *fe ) 
t 

fe->next=object->l_ features; 

object->l features=fe; 
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Znaczenie poszczególnych parametrów zostanie wyjaśnione w rozdziale 5. 
Natomiast phong decyduje o tym, czy powierzchnia budowanej bryły ma być 
wygładzana metodą Phonga w trakcie wizualizacji. Ten temat zostanie dokład- 
niej omówiony w punkcie 5.5. 


include "ray.h" 
OBJECT *add plane( char *name, int number, double a, 
double b, FEATURE *feature, int phong) 
t 
static int plane number = 0; 
// numer kolejny tworzonej struktury 
char *text = "plane"; 
/* numery od 0 do 999 */ 
/* UWAGA numer musi mieć mniej niż 3 litery */ 
char plane number text[4]; 
/* tablica czterech wskaźników do wierzchołków */ 
/* używana do wiązania wierzchołków ze ścianami 
i krawędziami */ 
VERTEX huge *vert[4]; 
OBJECT  *plane; 
FEATURE  *fe; 
VERTEX huge *vtab; 
EDGE huge *etab; 
FACE huge *ftab; 


if( ( plane=new object() ) == NULL ) 
return NULL; 
if( ( plane->name=name ) == NULL ) /* nadanie nazwy */ 
l 
if( ( plane->name=new text( 9 ) ) == NULL ) 
/* nowa , unikalna nazwa */ 
l 
delete object( plane ); 
return NULL; 
) 


lstrcat( plane->name , text ); 
itoa( plane number , plane number text , 10 ); 
lstrcat( plane->name , plane number _ text ); 

) 

plane_number++; 

plane->number=number; 

plane->n _vertices=4; 

plane->n_edges=4; 
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plane->n_faces=2; 
plane->axes=NULL; 


fe=new feature( feature->ktR , feature->ktG , 
feature->ktB ,feature->krR , 
feature->krG , feature->krB 
feature->kz , feature->m ); 


, 


if( fe == NULL ) 
t 
delete object( plane ); 
return NULL; 
) 
add feature to object( plane , fe ); 


// rozmiary prostokąta 
a/=2; 
b/=2; 


// wierzchołki 
vtab=new vertex( 4 ); 
if( vtab == NULL ) 
( 
delete object( plane ); 
return NULL; 


) 


// tworzymy 4-elementową listę wierzchołków, 
// nadając im odpowiednie współrzędne 

vert [0] = add vertex to object(vtab++, plane, —a,0, —b); 
vert [1] = add vertex to object(vtab++, plane, a,0, -b); 
vert [2] add vertex _to object(vtab++, plane, a,0, b); 
vert [3] = add vertex to object(vtab++, plane, —a,0, b); 


// krawędzie 

etab=new edge( 4 ); 

if( etab == NULL ) 

( 
delete object( plane ); 
return NULL; 
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// tworzymy 4-elementową listę krawędzi, 

// każda krawędź łączy dwa wierzchołki 
add edge to object( etab++, plane , vert[0] , vert[1] ) 
add _edge_to object( etab++, plane , vert[1] , vert[2] ); 
add_edge to object( etab++, plane , vert[2] , vert[3] ) 
add edge to object( etab++, plane , vert[3] , vert(0] ) 


// ściany 

ftab=new face( 2 ); 
if( ftab == NULL ) 
t 


delete _object( plane ); 
return NULL; 
> 
// tworzymy 4-elementową listę ścian. 
// Dodając każdą ścianę specyfikujemy, 
// które wierzchołki do niej należą, 
// podając składowe tablicy vert[] 
add_face_to object(  ftab++, plane, vert[0], vert[2], 
vert[3], fe, phong ); 
add_face_to object(  ftab++, plane, vert[0], vert[l], 
vert[2] , fe , phong ); 


return( plane ); 


) 


Do tworzenia list wierzchołków, krawędzi i ścian użyto opisanych wcześniej 
funkcji. Tablica vert[] zawiera cztery wskaźniki do struktur typu VERTEX 
opisujących wierzchołki. Jeśli tworzymy strukturę EDGE, to musimy jej podać 
jako parametry dwa wskaźniki do struktur opisujących wierzchołki ogranicza- 
jące daną krawędź. Te wskaźniki są zapamiętywane właśnie w tablicy vertl]. 
W trakcie tworzenia kolejnych list na bieżąco kontrolujemy, czy nie zabrakło 
pamięci w systemie. Służą do tego fragmenty kodu w postaci: 


ftab=new face( 2 ); 

if( ftab == NULL ) 

( 
delete_object( plane ); 
return NULL; 


) 


Jeżeli w pewnym momencie zabrakło pamięci (np. ftab == NULL), to wywołu- 
jemy funkcję delete_object usuwającą utworzoną uprzednio część obiektu. 
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a/2, 0, b/2 vertj2]-(a/2, 0, b/2) 


BZ 
/ 


vert[0]-(-a/2, 0, -b/2) 


vert[1]=(a/2, 0, -b/2) 


Rys. 2.11 


Procedura generująca prostopadłościan jest nieco bardziej skomplikowana. 


Wyjaśnijmy znaczenie poszczególnych parametrów: 


name — wskaźnik do tablicy zawierającej nazwę 
bryły. Jeżeli wywołamy funkcję z parame- 
trem name = NULL, to nazwa zostanie 


utworzona automatycznie; 


number — numer kolejny bryły; 

a — długość krawędzi równoległej do osi X ukła- 
du współrzędnych; 

b — długość krawędzi równoległej do osi Z; 

h — długość krawędzi równoległej do osi Y; 

fe, phong — patrz podrozdział 5.5. 


Funkcja zwraca wskaźnik do utworzonej struktury lub NULL, gdy brak 
pamięci. Centrum prostopadłościanu znajduje się w początku układu współ- 


rzędnych. 
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łinclude "ray.h" 


OBJECT *add_cube( char *name, int number, double a, 
double b, double h, FEATURE *feature, 


int phong ) 


static int cube _number=0; 
char *text="cube"; 


/* UWAGA numer musi mieć mniej niż 3 litery */ 
char  cube_number text[4]; /* numery od 0 do 999 */ 

/* tablica 8 wskaźników do wierzchołków */ 

/* używana do wiązania wierzchołków ze ścianami 


1 krawędziami */ 
VERTEX huge *vert[8]; 
OBJECT  *cube; 
FEATURE  *fe; 
VERTEX huge *vtab; 
EDGE huge *etab; 
FACE huge *ftab; 
if( ( cube=new object() ) == NULL ) 
return NULL; 


if( ( cube->name=name ) == NULL ) /* nadanie nazwy */ 


( 


if( ( cube->name=new text( 8 ) ) = 
t 

delete_object( cube ); 

return NULL; 
) 
lstrcat( cube->name , text ); 
itoa( cube number , cube number text , 10 ); 
lstrcat( cube->name , cube _number text ); 


NULL ) 


) 
cube _number++; /* kolejny prostopadłościan 
cube->number=numbex; 
cube->n_vertices=8; 
cube->n_edges=12; 
cube->n_faces=12; 
cube->axes=NULL; 
// Utworzenie i dołączenie do obiektu 
// struktury opisującej cechy jego powierzchni. 
fe=new feature( feature->ktR, feature->ktG, 
feature->ktB,feature->krR, 
feature->krG, feature->krB, 
feature->kz, feature->m ); 


* / 
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if( fe == NULL ) 
( 
delete _object( cube ); 
return NULL; 
) 


add_feature_to object ( cube , fe ); 


a/=2; 
b/=2; 
h/=2; 


// Rezerwacja pamięci dla listy 8 wierzchołków 
vtab=new_vertex( 8 ); 
1f( vtab == NULL ) 
( 
delete _object( cube ); 
return NULL; 
> 
// tworzenie listy wierzchołków 
vert[0)] = add_vertex to object(vtab+t+, cube,-a,-h,-b); 
vert[1] = add_vertex to object(vtab++, cube, a,-h,-b); 
vert[2] = add vertex to object(vtab++, cube, a,-h, b); 
vert(3] = add_vertex to object(vtab++, cube,-a,-h, b); 
vert[4] = add vertex to object(vtab++, cube,-a, h,-b); 
vert[5] = add _ vertex to object(vtab++, cube, a, h,-b); 
vert[6] = add_vertex to object(vtab++, cube, a, h, b); 
vert[7] = add vertex to object(vtab++, cube,-a, h, b); 


// krawędzie 
// Rezerwacja pamięci dla listy 12 krawędzi 
etab=new edge( 12 ); 
if( etab == NULL ) 
t 
delete _object( cube ); 
return NULL; 
> 
// tworzenie listy krawędzi 
add_edge to object( etab++, cube , vert[0] , vert[1] 
add edge to object( etab++, cube , vert[1i] , vert[2] 
add_edge_to object( etab++, cube , vert[2] , vert[3] 
add _edge_to object( etab++, cube , vert[3] , vert[0] 
add _edge to object( etab++, cube , vert[4] , vert[5] 
add_edge_to object( etab++, cube , vert[5] , vert[6] 
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add_edge to object( etab++, cube vert[6] vert[7) ); 
add_edge_ to object( etab++, cube vert[7] vert[4) ); 
add_edge to object( etab++, cube vert[0] vert[4) ); 
add_edge_to _object( etab++, cube vert[1) vert[5) ); 
add_edge to object( etab++, cube vert[2] vert[6) ); 
add_edge to object( etab++, cube vert[3] vert[7] ); 


// ściany 
// Rezerwacja pamięci dla listy 16 ścian 
ftab=new face( 12 ); 
if( ftab == NULL ) 
t 
delete object( cube ); 
return NULL; 
+ 


// tworzenie listy ścian 


add _face_to object ( ftab++, cube, vert[0], vert[5], 
vert[1], fe, phong ); 

add face_to object ( ftab++, cube, vert[0], vert[4], 
vert[5], fe, phong ); 

add_face to object ( ftab++, cube, vert[l], vert[6], 
vert[2], fe, phong ); 

add face to object ( ftab+t+, cube, vert[l], vert[5], 
vert[6], fe, phong ); 

add _face_to object ( ftab++, cube, vert[2], vert[7], 
vert[3], fe, phong ); 

add_face to object ( ftab++, cube, vert[2], vert[6]j, 
vert[7], fe, phong ); 

add face _to object ( ftab++, cube, vert[3], vert[4], 
vert[0], fe, phong ); 

add face _to object ( ftab++, cube, vert[3], vert[7], 
vert[4], fe, phong ); 

add _face_ to object ( ftab++, cube, vert[4], vert[6], 
vert[5], fe, phong ); 

add _ face to object ( ftab++, cube, vert[4], vert[7], 
vexrt[6], fe, phong ); 

add_face_to object ( ftab++, cube, vert[0], vert[2], 
vert[3], fe, phong ); 

add face to object ( ftab++, cube, vert[0], vert[1], 


vert[2], fe, phong ); 


return( cube ); 


) 
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Na koniec postaramy się wyjaśnić, w jaki sposób można automatycznie wyge- 
nerować dowolną figurę obrotową. Nie przytaczamy tekstu procedury add_fig- 
ure realizującej to zadanie ze względu na jej rozmiary (niemal 400 linii). 
Czytelników zachęcamy do przeanalizowania jej budowy. Funkcja add_figure 
znajduje się na dyskietce w pliku obiekty.c. 


Tworzenie bryły obrotowej przeanalizujemy na przykładzie funkcji add_sphere 
budującej sferę o zadanym promieniu. Rysunek 2.12 ilustruje sposób powsta- 
wania szkieletu bryły obrotowej. 


Rys. 2.12 


Obracając krzywą wokół osi Y otrzymujemy powierzchnię bryły obrotowej. 
Ponieważ dokładny opis matematyczny otrzymanej powierzchni jest skompli- 
kowany, a algorytmy wizualizacji tego typu powierzchni mało efektywne, 
jesteśmy zmuszeni do posługiwania się przybliżonym opisem. Taki przybliżony 
opis może polegać na przykład na wyodrębnieniu grupy punktów na powierz- 
chni i podzieleniu jej na trójkąty. Zamiast ciągłej powierzchni otrzymujemy 
w ten sposób zbiór trójkątnych ścian. 


Jak wybrać wierzchołki przyszłych ścian? 


Wystarczy przybliżyć krzywą pewną liczbą punktów (im więcej punktów, tym 
lepiej) i obracać każdy z nich wokół osi Z (patrz rys. 2.13) o pewien ustalony 
kąt (im mniejszy, tym dokładniejsze przybliżenie). Po każdym obrocie tworzy- 
my nowy punkt, aż otrzymamy siatkę, taką jak na rysunku 2.14. 
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Rys. 2.13 


Rys. 2.14 
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Sąsiednie punkty na siatce łączymy w trójkątne ściany (rys. 2.14). I w taki oto 
sposób otrzymaliśmy szkieletową reprezentację bryły obrotowej. Funkcja 
add_figure realizuje właśnie ten algorytm. 


ęłinclude "ray.h" 

OBJECT *add sphere ( char *sphere_name, int number, 
double r, int horizontal sections, 
int vertical sections, int stagger, 
FEATURE *fe, int phong ) 


VERTEX huge *vert; 

OBJECT *sp,*sphere; 

double alfa,delta_ alfa; 

double x,Yy; 

int i; 

static int sphere _number=0; 

char  *text="sphere"; 

char sphere number _text[4]; /* numery od 0 do 999 */ 
/* UWAGA numer musi mieć mniej niż 3 litery */ 

VERTEX huge *vtab; 


if( sphere _name == NULL ) /* nadanie nazwy */ 
( 
if( ( sphere name=new text( 10 ) ) == NULL ) 
/* nowa, unikalna nazwa */ 
return NULL; 
lstrcat( sphere name , text ); 
itoa( sphere number , sphere _number_text , 10 ); 
lstrcat( sphere_name , sphere _number text ); 
> 


sphere_number++; 
alfa=-M _PI/2; 
delta_alfa=M PI/vertical_ sections; 
if( ( sp=new object() == NULL ) 
( 

farfree( sphere_name ); 

return NULL; 


) 

vtab=new vertex( vertical _sections+l ); 
if( vtab == NULL ) 

t 


delete _object( sp ); 
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farfree( sphere_name ); 
return NULL; 


) 


// tworzymy krzywą składającą się 
// vertical sections wierzchołków. 


/* AARA */ 
for( 1=0 ; i<=vertical sections ; i++ ) 
( 


x=r*cos( alfa ); 
y=r*sin( alfa ); 
add_vertex to object( vtabt+, sp , x „,y , 0 ); 
alfa+=delta alfa; 
ł 
vert=sp->l_vertices; 
/* BBBB */ 


// stałe CLOSED są zdefiniowane w pliku ray.h 
sphere=add figure ( vert, sphere name, number, CLOSED, 
CLOSED, horizontal sections, 
stagger, fe, phong ); 
if( sphere == NULL ) 
t 
delete object( sp ); 
return NULL; 
) 
delete object( sp ); 
return( sphere ); 


Ciąg instrukcji od /* AAAA */ do /* BBBB */ generuje listę wierzchołków 
definiującą krzywą do obracania. W tym przypadku jest to półokrąg. Na 
podstawie tej listy funkcja add_figure buduje strukturę opisującą przybliżenie 
sfery. Po wykorzystaniu lista wierzchołków jest usuwana. 


2.3.4. Rysowanie brył 


Przedstawiliśmy procedury służące do tworzenia w pamięci komputera stru- 
ktur opisujących obiekty trójwymiarowe. Naszym następnym tematem będzie 
generowanie rysunków na podstawie informacji zawartej w tych strukturach. 
Przedstawimy funkcję realizującą to zadanie i przykładowy program. 
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Funkcja DrawObject rysuje szkielet bryły, którą opisuje struktura typu OB- 
JECT. 


void Drawobject( OBJECT *object ) 


( 
VERTEX huge *vertl, 
huge *vert2; 

EDGE huge *edge; 

if ( object == NULL ) 
return; 

// Rysowanie rzutu na płaszczyźnie XOY 

/* Wskaźnik edge przesuwamy po liście krawędzi. 
Każda struktura typu EDGE znajdująca się na tej 
liście zawiera wskaźniki do struktur VERTEX, 
opisujących współrzędne wierzchołków. Z pól x2, y2 
tych struktur pobieramy współrzędne x i y danego 
wierzchołka i rysujemy linię */ 

edge=object->l _edges; 

while( edge != NULL ) 

( 
vertl=edge->vertl; 
vert2=edge->vert2; 
MoveTo( ViewDC, vertl->x2, vertl->y2 ); 
LineTo( ViewDC, vert2->x2, vert2->y2 ); 
edge=edge->next; 

> 

> 


Zanim użyjemy tej funkcji, musimy wywołać procedurę LogicalToScreen reali- 
zującą tzw. okienkowanie (ang. windowin$). Jest to proces zamiany współrzęd- 
nych punktów z logicznych na współrzędne urządzenia graficznego, którym jest 
w naszym przypadku ekran komputera. Przez współrzędne logiczne rozumie- 
my współrzędne punktów liczone w abstrakcyjnym układzie odniesienia. Takie 
współrzędne są zapisane w strukturach opisujących obiekty. Okienkowanie 
umożliwia prezentację wybranego, prostokątnego wycinka przestrzeni na 
ekranie komputera. Wybór tego wycinka odbywa się poprzez podanie współ- 
rzędnych lewego dolnego wierzchołka (Xmin, Ymin) i prawego górnego (Xmax, 
Ymax) (rys 2.15). Ekran został przedstawiony na rysunku 2.16. Współrzędne 
na ekranie są liczone w pikslach. Wysokość ekranu wynosi xScreen piksli, a 
szerokość yScreen. 
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Rys. 2.15 


0 xScreen 


yScreen 


Rys. 2.16 
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Wnikliwy Czytelnik zauważył zapewne, że wywołując procedurę okienkowania 
podajemy tylko dwie współrzędne, pomijając współrzędną z. Wynika to stąd, 
że pragnąc narysować wizerunek trójwymiarowego obiektu na dwuwymiaro- 
wym ekranie, musimy przejść od trójwymiarowych współrzędnycjh logicznych 
do dwuwumiarowych współrzędnych ekranu (rys. 2.17). Pominięcie jednej 
współrzędnej w trakcie wyświetlania obrazów brył jest najprostszą formą 
rzutowania. Inne sposoby rzutowania zostaną opisane w rozdziale 4. W tym 
przypadku pomijając współrzędną z punktów rzutujemy je na płaszczyznę 
XOY, a potem realizujemy okienkowanie. Jednak nic nie stoi na przeszkodzie, 
aby pominąć współrzędną x (rzut na płaszczyznę YOZ) lub y (rzut na płaszczy- 
znę XOZ). 


Rys. 2.17 


Funkcja LogicalToScreen wykonuje okienkowanie dla wszystkich elementów 
listy wierzchołków wskazywanej przez parametr *vert. Okienkowanie polega 
na obliczeniu współrzędnych ekranowych (pola x2, y2 każdej struktury VER- 
TEX) na podstawie współrzędnych globalnych (pola globalx, globaly w struktu- 
rze VERTEX). Nie przytaczamy wzorów na okienkowanie, ponieważ są one 
zapisane bezpośrednio w ciele funkcji LogicalToScreen. 
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void LogicalToScreen (  VERTEX huge *vert, double xScree: 
double yScreen, double Xmin, 
Double Ymin, double Xmax, 
double Ymax ) 


t 
double sx, Sy; 
SX = ( xScreen - 0 )/( Xmax - Xmin ); 
sy = ( yScreen - 0 )/( Ymax — Ymin ); 
while( vert != NULL ) 
l 
if ( sx < sy ) 
1 
vert->x2 = sx*( ( vert->globalx ) - Xmin ); 
vert->y2 = yScreen/2 - sx*( ( vert->globaly ) -— 
( Ymin + Ymax )/2 ); 
) 
else 
( 
vert->x2 = xScreen/2 + sy*( ( vert->globalx ) -— 
( Xmin + xXmax )/2 ); 
vert->y2 = yScreen — sy*( ( vert->globaly ) - 
Ymin ); 
) 
vert=vert->next; 


) 


A teraz przykłady. Wszystkie przykładowe programy były napisane w Bor- 
land C++, dlatego wykorzystujemy w nich funkcje otwierania i zamykania 
grafiki z pakietu firmy Borland. Na dyskietce znajdują kody źródłowe wszy- 
stkich przykładów, które możne skompilować za pomocą kompilatora Bor- 
land C++. Oczywiście można skompilować te programy także za pomocą innych 
kompilatorów, jednak może się to wiązać z koniecznością zmiany wywołań 
takich funkcji graficznych, jak rysowanie odcinków, otwieranie i zamykanie 
grafiki itp. 


Na początku każdego przykładu są włączane pliki ray.h, listy.c, mat.c, obie- 
kty.c, rysuj.c, zawierające definicje odpowiednich typów danych i kody 
źródłowe wszystkich funkcji, które opisujemy w tej książce. 
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Przykład 1 


Poniższy program znajduje się na dyskietce w pliku progl.c. Tworzy sześcian 
o boku 50 i wyświetla go na ekranie komputera (rys 2.18). Wynikiem działania 
programu jest kwadrat. Właśnie tak wygląda rzut sześcianu na płaszczyznę 
XOY. 


// progl.c 

łinclude "ray.h" 
łinclude "listy.c" 
łinclude "mat.c" 
include "obiekty.c" 
łinclude "rysuj.c" 


void main( void ) 


( 
OBJECT *szescian; 
FEATURE *f; 
int gdriver = DETECT, gmode, errorcode; 


// inicjalizacja grafiki 
initgraph(6gdriver, śgmode, "c:NNborlandciibgi"); 
errorcode = graphresult(); 

if (errorcode != grOok) /* błąd */ 

( 
printf("Błąd grafiki o kodzie: %sin", 
grapherrormsg(errorcode)); 
printf("Naciśnij dowolny klawisz."); 
getch (); 
exit(0); 

) 


// parametry powierzchni 
f = new feature( 0, 0, 0, 1, 1, 1, 0, 1); 
if( f == NULL ) 
exit(0); 
// utworzenie struktur opisujących sześcian 
szescian = add _cube( NULL, 0, 50, 50, 50, f, OFF ); 
if( szescian == NULL ) 
exit(0); 
LogicalToscreen( szescian->l vertices, getmaxx(), 
getmaxy(), —100, -100, 100, 100 ); 
Drawobject( szescian ); 
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// zamknięcie grafiki 
getch(); 
closegraph (); 


Rys. 2.18 
Przykład 2 


Program prog2.c rysuje rzut kuli (rys. 2.19). Okno jest teraz położone inaczej 
niż w poprzednim przykładzie, dlatego widzimy tylko ćwiartkę kuli. Jeżeli 
wywołanie procedury okienkowania będzie następujące: 


LogicalToScreen(kula->l_vertices, getmaxx(), getmaxy(), -100, 100, 100, 100);, 
to szkielet kuli zostanie narysowany w centrum ekranu. 


// prog2.c 

tinclude "ray.h" 
include "listy.c" 
$include "mat.c" 
include "obiekty.c" 
tłinclude "rysuj.c" 
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void main( void ) 
( 
OBJECT *kula; 


FEATURE *f; 
int gdriver = DETECT, gmode, errorcode; 


// inicjalizacja grafiki 
initgraph(śsgdriver, śgmode, "c:Niborlandciibgi"); 
errorcode = graphresult(); 
if (errorcode != grok) /* błąd */ 


( 
printf( "Błąd grafiki o kodzie: %sin", 
grapherrormsg(errorcode) ); 
printf("Naciśnij dowolny klawisz."); 
getch(); 
exit(0); 
) 


// parametry powierzchni 
f = new feature( 0, 0, 0, 1, 1, 1, 0, 1); 
if( f == NULL ) 
exit(0); 
// utworzenie struktur opisujących sześcian 
kula = add sphere( NULL, 0, 50, 10, 10, OFF, £f, OFF ); 


if( kula == NULL ) 
exit(0); 
LogicalToScreen( kula->l_vertices, getmaxx(),getmaxy(), 


0, 0, 100, 100 ); 
DrawObject( kula ); 


// zamknięcie grafiki 
getch(); 
closegraph(); 
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ROZDZIAŁ 3 


Przekształcenia w przestrzeni 
trójwymiarowej 


W poprzednich rozdziałach zapoznaliśmy Czytelnika z podstawowymi pojęcia- 
mi geometrii w przestrzeni i sposobami reprezentacji trójwymiarowych figur 
w pamięci komputera. Ten rozdział pragniemy poświęcić omówieniu zagadnień 
dotyczących manipulowania bryłami. Przedstawimy wzory pozwalające na 
przesuwanie, obracanie, powiększanie i pomniejszanie dowolnych obiektów. 
Wprowadzimy pojęcie układu lokalnego i podamy funkcje pozwalające na 
przejście z jednego układu współrzędnych do innego. Ważniejsze zagadnienia 
zilustrujemy odpowiednimi przykładami. 


Rys. 3.1 
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3.1. Układ globalny i układ lokalny 


Na początku zdefiniujemy pewne pojęcia, którymi będziemy się dalej posługi- 
wać. W trójwymiarowej przestrzeni kartezjańskiej możemy wyróżnić pewien 
układ odniesienia, względem którego będziemy określać położenie wszystkich 
obiektów (brył, lamp, obserwatora). Układ ten będziemy nazywać układem 
globalnym, w skrócie UG. Przyjmujemy, że jest to lewostronny układ współ- 
rzędnych, taki jak na rysunku 3.1. 


Układ globalny możemy porównać do powierzchni naszej planety. Położenie 
takich obiektów, jak góry, domy, drzewa, określane względem Ziemi, nie 
zmienia się. Także wzajemne położenie tych obiektów względem siebie nie 
ulega zmianom. Są jednak obiekty przemieszczające się względem tak umiej- 
scowionego układu odniesienia. Na przykład samochód, poruszając się po 
drodze, przesuwa się względem powierzchni Ziemi. Można znaleźć taki układ, 
w którym samochód nie porusza się. Powinien to być układ poruszający się w 
tym samym kierunku co samochód, z identyczną prędkością. Takim układem 
jest na przykład układ związany z nadwoziem pojazdu. Ten układ odniesienia 
będziemy nazywać lokalnym. Opisaną sytuację ilustruje rysunek 3.2. 


Rys. 3.2 


Zauważmy, że układ lokalny może być nie tylko przemieszczany względem 
globalnego, ale także obracany. Dzieje się tak, gdy samochód pokonuje zakręt. 
Wraz z samochodem obraca się jego układ lokalny (rys 3.3a, 3.3b). Zaistniałą 
Bytuację możemy opisać w nieco inny sposób. Związując układ lokalny 2 
samochodem powodujemy, że położenie samochodu w układzie lokalnym nie 
zmienia się. Jeżeli obrócimy układ lokalny względem globalnego (powierzchnia 
Ziemi), to obraca się także obiekt nieruchomy w układzie lokalnym (samochód). 
Przesuwając układ lokalny względem globalnego, przesuwamy także wszystkie 
obiekty znajdujące się w tym układzie. 
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Rys. 3.3 


3.2. Scena 


Zgodnie z tytułem naszej książki pragniemy przedstawić metody pozwalające 

„na symulowanie pewnych elementów otaczającego nas świata w pamięci kom- 
putera i uzyskiwanie obrazów na podstawie takiego opisu. Sceną będziemy 
dalej nazywać zbiór wszystkich elementów tworzących tę symulację. Scena jest 
więc czymś w rodzaju wirtualnego, matematycznego świata zbudowanego 
z pewnej liczby obiektów. Współrzędne tych obiektów będą reprezentowane 
w układzie globalnym. Z każdym obiektem będzie związany układ lokalny. 


Podstawowymi elementami składowymi sceny są oczywiście bryły. Każda bryła 
jest złożona ze zbioru punktów będących jej wierzchołkami. Dwa sąsiadujące 
ze sobą wierzchołki tworzą krawędź. Każda ściana jest definiowana za pomocą 
dokładnie trzech wierzchołków. Takie założenie nie ogranicza klasy brył 
możliwych do zobrazowania, ponieważ każdą ścianę (wielokąt) można zredu- 
kować do zbioru trójkątów za pomocą dostępnych algorytmów. Należy jednak 
zaznaczyć, że nie zawsze jest to zadanie trywialne. 


Kolejnym elementem sceny jest obserwator. Możemy wyobrazić go sobie jako 
aparat fotograficzny wysłany do wirtualnego świata z zadaniem uzyskania jego 
obrazów. Obserwator może być poddawany takim samym transformacjom 
(przesunięcie, obroty), jak układy lokalne lub bryły. Na scenie mogą znaleźć 
się także źródła światła — lampy. 
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Struktura SCENE służy do przechowywania informacji o scenie: 


typedef struct sc // struktura opisuje scenę, czyli zbio! 


t // wszystkich obiektów 
char *name; // wskaźnik do nazwy sceny 
AXES *] _ axes; // wskaźnik do pierwszego elementu 


// listy układów lokalnych 
LAMP *l lamps; // lista lamp 
OBSERVER *observer;  // struktura opisująca kamerę 


double R; // trzy składniki barwy oświetlenia tłż 
double G; 
double B; 


int number_axes; 
int number objects;  // liczba brył znajdujących się 
// na scenie 


int number faces; // liczba wszystkich ścian 
int number _edges; // liczba wszystkich krawędzi 
int number _vertices; // liczba wszystkich wierzchołków 
int number_lamps; // liczba lamp 
) SCENE; 


3.3. Transformacje w przestrzeni trójwymiarowej 


Podstawowymi przekształceniami w przestrzeni są przesunięcie, skalowanie 
i obroty. Wszystkie te transformacje mogą odnosić się zarówno do obiektów, 
jak i do układów współrzędnych. W kolejnych punktach przedstawimy wzory 
i macierzowe zapisy poszczególnych przekształceń. 


3.3.1. Współrzędne jednorodne 


Pojęcie współrzędnych jednorodnych często pojawia się w literaturze poświę- 
conej grafice komputerowej. W naszej książce będziemy się nim rzadko posłu- 
giwać. Jednak dla porządku postaramy się wyjaśnić, jakie zalety ma używanie 
współrzędnych jednorodnych. 


Jak wiemy, chcąc określić położenie punktu w przestrzeni, należy podać jego 
trzy współrzędne liczone względem ustalonego układu współrzędnych. Każde” 
mu punktowi o współrzędnych x, y, z możemy przyporządkować pewien wektor 
[x,y,z]. Obrót tego punktu wokół osi Z jest transformacją, powoduje zmianę jeBo 
współrzędnych. Nowe współrzędne możemy wyznaczyć na podstawie równań 
wiążących kąt obrotu ze starymi współrzędnymi: 
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x =xcosó-ysinQ$ 
y =xsin$ +yCcos6$ 
z=0 


Ponieważ taka forma zapisu jest dość nieczytelna, możemy posłużyć się zapi- 
sem macierzowym: 


cososino O 
[x ,y ,z']=[x,y,z]|-singcoso 0 
0 0 1 


Macierz po prawej stronie równania jest macierzą transformacji. Taki zapis 
ma także inną zaletę. Złożenie dwóch przekształceń (na przykład obracamy 
punkt wokół osi X, a następnie wokół osi Z) można zapisać w postaci pojedynczej 
macierzy, która jest iloczynem macierzy tych przekształceń. 


Także skalowanie można zapisać w formie równania macierzowego. Jedynym 
przekształceniem, którego nie można opisać jako iloczyn wektora przez ma- 
cierz, jest translacja, czyli przesunięcie. 


Aby umożliwić wykonanie tego zadania; należy dodać do macierzy przekształ- 
ceń jedną kolumnę i jeden wiersz, a do wektorów opisujących położenie punktu 
dodatkową współrzędną-stałą. Teraz punktowi o współrzędnych x,y,z odpowia- 
da wektor P[x,y,z,1]. Współrzędne punktu wraz z dodatkową stałą nazywamy 
współrzędnymi jednorodnymj. Dowolna transformacja może spowodować 
zmianę jedynie trzech pierwszych współrzędnych. W dalszych rozważaniach 
będziemy przyjmować, że punkt P[x,y,z',1] jest obrazem punktu P[x,y,z,1] 
(czyli punktem P po wykonaniu danej transformacji). 


3.3.2. Przesunięcie 


Przesunięcie punktu w przestrzeni z jednej pozycji do innej polega na dodaniu 
współrzędnych wektora przesunięcia V do pierwotnych współrzędnych tego 
punktu: 


X =X+X, 
Y =Y+JM 
Z =Z+Z, 
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W notacji macierzowej przesunięcie opisujemy równaniem: 


, , , 1 —_ 1 0 l 
[x JJ ;Z ; ]=[x,y,z, ] 0 0 
XJ 


PhRy,z,1] 


/ 
4 
/ VMxXvyv.Zy 1] 


Pfx,y,z,1] 


Rys. 3.4 


Powyższe wzory pozwalają nam wyznaczyć współrzędne jednego punktu. W 
jaki sposób można przesunąć całą bryłę? Oczywiście należy poddać transfor- 
macji wszystkie wierzchołki przesuwanej bryły. To zadanie realizuje następu- 
jąca funkcja: 


void move( VERTEX huge *vert, int where, double move_x, 
double move_y, double move_z ) 


( 
while( vert != NULL ) 
( 
if( where == GLOBAL ) 
( 


/* współrzędne po przesunięciu */ 
vert->globalxt=move_x; 
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vert->globaly+=move_y; 
vert->globalz+=move_z; 
) 
else 
t 
/* współrzędne po przesunięciu */ 
vert->localx+=move_x; 
vert->localy+=move_y; 
vert->localz+=move_z; 


) 


vert=vert->next; 


Vert jest wskaźnikiem do listy wierzchołków, które mają zostać przesunięte o 
wektor move_x, move_y, move_z. 


Jeżeli parametr where = GLOBAL, to modyfikowane są współrzędne globalne 
wierzchołka. Jeżeli where = LOCAL, to modyfikowane są współrzędne wierz- 
chołka liczone w układzie lokalnym. (Stałe GLOBAL i LOCAL, podobnie jak 
inne stałe, są zdefiniowana w pliku ray.h). 


Przykład 

Jeżeli utworzymy sześcian za pomocą funkcji add_cube, to przesunięcie go o 
wektor [50,50,0] można wykonać w następujący sposób: move(cube->l_vertices, 
GLOBAL, 50, 50, 0); 

3.3.3. Skalowanie 


Skalowanie jest realizowane względem początku danego układu współrzęd- 
nych, co zostało przedstawione na rysunku 3.5. 


Skalowanie opisują równania: 


? 

X =x*S, 
> _— Ó 
JJ Sy 
z =Z*5, 
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Rys. 3.5 


Korzystając z funkcji scale skalowanie możemy wykonać w układzie globalnym 
lub lokalnym (w zależności od parametru where). 


void scale( VERTEX huge *vert, int where, double scale_x, 
double scale_y, double scale_z ) 


l 
while( vert !i= NULL ) 
A 
if( where == GLOBAL ) 
( 


/* współrzędne po operacji */ 
vert->globalx*=scale_x; 
vert->globaly*=scale_y; 
vert->globalz*=scale_z; 
) 


else 
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/* współrzędne po operacji */ 
vert->localx*=scale_x; 
vert->localy*=scale_y; 
vert->localz*=scale_z; 


) 


vert=vert->next; 


) 


Także ta funkcja wykonuje transformację na wszystkich wierzchołkach znaj- 
dujących się na liście wskazywanej przez *vert. Jeżeli użyjemy tej funkcji do 
powiększania/pomniejszania brył w układzie globalnym, to musimy liczyć się 
z tym, że dana bryła zostanie przesunięta w kierunku od/do środka układu 
współrzędnych (patrz rys 3.6a). Jak łatwo zauważyć, taki efekt nie pojawi się 
przy skalowaniu bryły, której wierzchołki są rozmieszczone symetrycznie 
względem środka układu współrzędnych; innymi słowy, gdy środek ciężkości 
bryły ma współrzędne [0,0,0](rys. 3.6b). 


Rys. 3.6 
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Wprowadzenie układów lokalnych pozwala na proste usunięcie tej niedogod- 
ności. Skalowanie bryły wykonujemy w układzie lokalnym (środek ciężkości 
bryły leży w centrum tego układu), a następnie wykonujemy przejście od 
współrzędnych lokalnych do współrzędnych w układzie globalnym. W efekcie 
obiekt zostanie powiększony lub pomniejszony, ale jego środek będzie spoczy- 
wał w tym samym miejscu, co przed skalowaniem (rys 3.6c). Ostatnia operacja 
zostanie szczegółowo opisana w punkcie 3.4.2. 


3.3.4. Obroty 


Podstawowymi obrotami w przestrzeni trójwymiarowej są obroty wokół osi 
układu współrzędnych. Takie przekształcenia zmieniają tylko dwie z trzech 


współrzędnych punktu. Na przykład podczas obrotu punktu wokół osi X nie 
ulega zmianie współrzędna x obracanego punktu. 


Przyjęliśmy, że we wszystkich poniższych wzorach kąty obrotu mają dodatni 
zwrot. Co to oznacza? Dodatni zwrot kąta obrotu wokół danej osi, to zwrot 
przeciwny do kierunku ruchu wskazówek zegara, gdy patrzymy na układ 
współrzędnych w ten sposób, że oś obrotu jest skierowana zgodnie z kierunkiem 
patrzenia. Zwroty tych kątów są przedstawione na rys 3.7. 


Y 


Rys. 3.7 


i — 


66 


M. Domaradzki, R. Gembara 


Obrót o kąt $ wokół osi X układu współrzędnych (rys. 3.8a) opisują wzory: 


X '=xX 
y =ycos$-zsinó 
z =ysinQ +ZCos$ 
1 0 0 0 
, , , 0 COS SIN 0 
[x ,Y »Z ,1]=[x,y,z,l] 0 — sin pcoso 0 
0 0 0 1 


Plx,y.z,1] 


-g PDOY Z] 


Rys. 3.8a 


Punkt P[x,y,z,1] jest przekształcany na P[x,y',z,1] podczas obrotu o kąt a 
wokół osi Y (rys. 3.8b): 


x =ZzZSINa+XxXCcosa 
9 
) =Yy 
Z”=ZCOSa—XSIN a 
cosa 0 -—sina 0 


SEE 01 00 
[x',y',z',1]=[x,y,z,l] sina 0 cosa 0 
0 0 0 1 
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Plxy,z,1] 


7 Ppey,2,1] 


Rys. 3.8b 
Obrót o kąt 8 wokół osi Z układu współrzędnych (rys. 3.8c) opisują wzory: 


x =xcosB—ysin$B 
y =xsinf ycosf 


Z'=Z 
cosBsin8 00 
, , , — SI COS 00 
[x',y',z”,1]=[x,y,z,1]| "97 s5 15 
0 001 


Korzystając z poniższych funkcji należy pamiętać, że w zależności od parame- 
tru where obrót jest wykonywany w układzie lokalnym (where = LOCAL) lub 
globalnym (where = GLOBAL). 


/* Obrót każdego elementu listy wierzchołków wskazywanej 
przez vert wokół osi X o kąt rot x stopni */ 
void rot _x( VERTEX huge *vert, int where, double rot_x ) 
( 

double x_old,y old,z old; 

double c,s; 


c=cos( rot_x*M PI/180 ); 
s=sin( rot_x*M PI/180 ); 
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4 Pi y,Z,1] 


Plx„y,z,1] 


while( vert != NULL ) 

( 
if( where == GLOBAL ) 
( 


// zapamiętanie współrzędnych punktu przed obrotem 
x old=vert->globalx; 
y_.old=vert->globaly; 
z_old=vert->globalz; 
vert->globalx=x old; 
vert->globaly=y oldtc-z old*s; 
vert->globalz=y old*stz_old*c; 


) 


else 


( 


// zapamiętanie współrzędnych punktu przed obrotem 


x _old=vert->localx; 
y_old=vert->localy; 
z _old=vert->localz; 
vert->localx=x_ old; 
vert->localy=y old*c-z _old*s; 
vert->localz=y old*s+z_ old*c; 


) 


vert=vert->next; 
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> 


/J* Obrót każdego elementu listy wierzchołków wskazywanej 
przez vert wokół osi Y o kąt rot _y stopni */ 


void rot_y( VERTEX huge *vert, int where, double rot _y ) 
l 

double x old,y old,z old; 

double c,s; 


c=cos( rot_y*M PI/180 ); 
s=sin( rot_y*M PI/180 ); 
while( vert != NULL ) 
t 
if( where == GLOBAL ) 
l 
// zapamiętanie współrzędnych punktu przed obrotem 
x _old=vert->globalx; 
y_.old=vert->globaly; 
z_old=vert->globalz; 
vert->globalx=z _old*s+x old*c; 
vert->globaly=y old; 
vert->globalz=z _old*c-x old*s; 
) 
else 
t 
// zapamiętanie współrzędnych punktu przed obrotem 
x _old=vert->localx; 
y_ old=vert->localy; 
z _old=vert->localz; 
vert->localx=z_old*s+x oldtc; 
vert->localy=y old; 
vert->localz=z old*c-x oldts; 
> 


vert=vert->next; 


) 


/* Obrót każdego elementu listy wierzchołków wskazywanej 
przez vert wokół osi Z o kąt rot z stopni */ 


void rot_z( VERTEX huge *vert, int where, double rot_z ) 


OZZZZNNNZNZZZZNNN NZ EN NNNNENNNNNNNNENNNEEZENNNNNNNNNNNNENNN A 


70 


M. Domaradzki, R. Gembara 


double x_old,y old,z old; 
double c,s; 
c=cos( rot_z*M PI/180 ); 
s=sin( rot_z*M PI/180 ); 
while( vert != NULL ) 
( 
i1f( where == GLOBAL ) 
( 
// zapamiętanie współrzędnych punktu przed obrotem 
x _old=vert->globalx; 
y_old=vert->globaly; 
z_old=vert->globalz; 
vert->globalx=x old*c-y old*s; 
vert->globaly=x old*s+ty old*c; 
vert->globalz=z_ old; 
) 
else 
( 
// zapamiętanie współrzędnych punktu przed obrotem 
x _old=vert->localx; 
y_.old=vert->localy; 
z _old=vert->localz; 
vert->localx=x old*c-y old*s; 
vert->localy=x old*sty old*c; 
vert->localz=z_ old; 


) 


vert=vert=>next; 


) 


W literaturze poświęconej grafice komputerowej opisano także obroty wokół 
dowolnie wybranej osi. Taka operacja jest realizowana jako złożenie kilku 
obrotów wokół osi układu odniesienia. Logicznie odpowiada to przejściu do 
układu lokalnego związanego z osią obrotu, wykonaniu obrotu i powróceniu do 
układu globalnego. 


Nie będziemy przedstawiać explicite równań opisujących tę transformację. 
Opisywane przez nas procedury pozwalają na definiowanie nowych układów 
współrzędnych — układów lokalnych i wykonywanie na nich dowolnych trans- 
formacji, a także zamianę współrzędnych z układu globalnego do lokalnego 
i odwrotnie. Osie układu lokalnego mogą być inaczej zorientowane niż osie 
układu globalnego. Istnieje także możliwość obrotu bryły wokół osi układu 
lokalnego, co jest de facto obrotem wokół dowolnie wybranej osi. 
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3.3.5. Przykłady wykorzystania funkcji realizujących transformacje 


Przykład 1 


Spróbujmy zbudować schody. Do tego celu wykorzystamy funkcję add_cube 
tworzącą prostopadłościany. Niech nasze schody mają cztery stopnie. Każdy 
stopień ma szerokość 40, wysokość 20 i głębokość 80 jednostek. Pierwszy 
stopień będzie leżał w środku układu współrzędnych, następne będą przesu- 
nięte względem swoich poprzedników o 20 jednostek w górę i o 30 jednostek w 
prawo. Oto program realizujący nasze zadanie: 


// prog3.c 

include "ray.h" 
kinclude "listy.c" 
include "mat.c" 
łinclude "obiekty.c" 
ęinclude "rysuj.c" 


void main( void ) 
t 
int i; 
OBJECT *stopnie [4]; 
FEATURE  *f; 
int gdriver = DETECT, gmode, errorcode; 


// inicjalizacja grafiki 
initgraph(śsgdriver, śgmode, "c:Niborlandciibgi"); 


// parametry powierzchni 
f = new feature( 0, 0, 0, 1, 1, 1, 0, 1); 
if( f == NULL ) 
exit(0); 
// utworzenie czterech stopni 
for( i=0; i<=3; i++ ) 
stopnie[i] = add cube( NULL, 0, 40, 80, 20, £f, OFF); 
// przesunięcie trzech stopni 
for( i=l; i<=3; i+ ) 
move( stopnie[i]->l vertices, GLOBAL, 30*i, 20*i, 0); 


// obliczamy współrzędne do rysowania i wyświetlamy 
// kolejno cztery stopnie 
for( i=0; i<=3; i++ ) 
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( 
LogicalToScreen( Sstopnie[i] ->l_vertices, getmaxx(), 
getmaxy(), 0, -30, 100, 100 ); 
Drawobject( stopnie([i] ); 
> 


// zamknięcie grafiki i skasowanie obiektów 
getch(); 
closegraph(); 
for( i=0; i<=3; i++ ) 
delete_object( stopnie(i] ); 


Tablica stopnie[] zawiera wskaźniki do struktur opisujących 4 prostopadłościa- 
ny. Po utworzeniu wszystkie prostopadłościany znajdują się w tym samym 
położeniu (rys 3.9a). Pierwszego prostopadłościanu nie przesuwamy, drugi 
przesuwamy o wektor [30,20,0] (patrz rysunek 3.9b). Kolejne bryły przemiesz- 
czamy o wektory [60,40,0] i [90,60,0]. Wynikiem działania powyższego progra- 
mu jest rysunek 3.9c. 


Rys. 3.9a 
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Rys. 3.9b 


Rys. 3.9c 
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Przykład 2 


W poprzednim przykładzie wykorzystaliśmy tylko jedną transformację. Poniż- 
szy przykład pokazuje, w jaki sposób możemy za pomocą opisanych funkcji 
wykonywać złożenia przekształceń. Spróbujmy wykonać rysunek spiralnych 
schodów. 


Każdy nowo utworzony prostopadłościan znajduje się w centrum globalnego 
układu współrzędnych. Procedura dodawania kolejnych stopni do spiralnych 
schodów jest następująca. Tworzymy i-ty prostopadłościan o wymiarach 
80x40x10. Przesuwamy go o 100 jednostek w kierunku osi Z. Przesunięty 
prostopadłościan obracamy o wielokrotność pewnego ustalonego na początku 
kąta (np. 20 stopni), czyli o kąt i*20 stopni. Na koniec przesuwamy otrzymaną 
bryłę o 10 jednostek (wysokość stopnia) w kierunku osi Y. Przed wyświetleniem 
należy jeszcze wywołać procedurę okienkowania dla wszystkich prostopadło- 
ścianów. 


// prog4.c 

include "ray.h" 
finclude "listy.c" 
łinclude "mat.c" 
include "obiekty.c" 
łinclude "rysuj.c" 


tdefine LICZBA STOPNI 20 
łdefine KAT _ OBROTU 15 


void main( void ) 
( 
int i; 
OBJECT . *stopnie[LICZBA STOPNI]; 
FEATURE  *f; 
int gdriver = DETECT, gmode, errorcode; 


// inicjalizacja grafiki 
initgraph(Sgdriver, Sgmode, "c:NNlborlandclibgi"); 


// parametry powierzchni 
f = new _feature( 0, 0, 0, 1, 1, 1,0, 1 ); 
if( f == NULL ) 
exit(0); 
// utworzenie stopni 
for( i=0; 1<LICZBA STOPNI; i++ ) 
t 
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stopnie[i] = add _cube( NULL, 0, 40, 80, 10, f, OFF); 
if( stopnie[i] == NULL ) 
exit(0); 

> 


// przesunięcie każdego stopnia o 120 jednostek 

// w kierunku osi Z 

// obrót wokół osi Y o i*KAT OBROTU 

// przesunięcie o 20 jednostek w górę 

for( i=0; i<LICZBA STOPNI; i++ ) 

( 
move( stopnie[i]->l vertices, GLOBAL, 0, 0, 100 ); 
rot_y( stopnie[i]->l vertices, GLOBAL, i*KAT_ OBROTU); 
move( stopnie[i]->l vertices, GLOBAL, 0, 10*i, 0 ); 
77 11 

) 


// obliczamy współrzędne do rysowania i wyświetlamy 
kolejno wszystkie stopnie 
for( i=0; i<LICZBA STOPNI; i++ ) 
t 
LogicalToScreen( stopnie[i]->l vertices, getmaxx(), 
getmaxy(), -200, -30, 200, 3600 ); 
Drawobject( stopnie[i] ); 
) 


// zamknięcie grafiki i skasowanie obiektów 
getch (); 
closegraph (); 
for( i=0; 1<LICZBA STOPNI; i++ J 
delete _object( stopnie[i] ); 


Na rysunku 3.10 przedstawiono efekt działania programu. Jest to rzut wszy” 
stkich brył na płaszczyznę OXY, dlatego nie widać ich przestrzennego rozmie- 
szczenia. Aby uzyskać lepszy efekt wizualny, możemy obrócić nasze schody na 
przykład o 30 stopni wokół osi X. Realizujemy to poprzez dodanie następujace 
linii do programu w miejscu oznaczonym przez !!!: 


rot_x( stopnie[i]->l vertices, GLOBAL, -30 ); 


Otrzymamy wtedy następujący rysunek (rys. 3.11): 
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Rys. 3.10 


| 
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4, 


Rys. 3.11 


77 


Realistyczna grafika 3D 


3.4. Zmiana układu współrzędnych 


3.4.1. Układ lokalny 


W grafice komputerowej wprowadzenie pojęcia układów lokalnych pozwala na 
uproszczenie opisu skomplikowanych operacji na obiektach. Jest jeszcze jedna 
ważna korzyść wynikająca ze stosowania takiej formy opisu. Dla każdego 
układu lokalnego możemy utworzyć wiele nowych układów lokalnych wzglę- 
dem niego. Podobnie możemy postąpić z każdym z tych nowych układów etc. 
W ten sposób powstaje pewna hierarchiczna struktura zależności jednych 
obiektów od innych. Takie struktury, wbrew pozorom, spotykamy w rzeczywi- 
stym świecie nadzwyczaj często. Wystarczy spojrzeć w lustro. To nie dowcip, 
każdy człowiek może zostać opisany jako hierarchiczna struktura. W trakcie 
ruchu człowiek przesuwa się względem podłoża, ręce poruszają się względem 
tułowia (a więc także wzglęm podłoża), a zwisająca swobodnie dłoń zmienia 
swe położenie względem ręki. O ile łatwo znaleźć równanie toru ruchu w 
układzie związanym z powierzchnią Ziemi dla tułowia, nieco trudniej dla 
nadgarstka, to rozwiązanie tego problemu na przykład dla końca palca wska- 
zującego jest dość kłopotliwe, 


Z pomocą przychodzą nam układy lokalne. Zauważmy, że w układzie związa- 
nym z dłonią końcówka palca wskazującego porusza się po okręgu. Z kolei dłoń 
w układzie związanym z nadgarstkiem także porusza się po okręgu (z pewnymi 
ograniczeniami), a nadgarstek zatacza wycinek okręgu w układzie związanym 
z łokciem. Podobnie porusza się łokieć w układzie związanym z tułowiem. Za 
pomocą serii obrotów w powiązanych ze sobą układach lokalnych można opisać 
skomplikowany tor ruchu końcówki palca. Oczywiście niezbędna jest tu jeszcze 
znajomość relacji zachodzących pomiędzy współrzędnymi punktu znajdującego 
się w danym układzie względem współrzędnych tego samego punktu wyrażo- 
nych w innym (lokalnym) układzie współrzędnych. Te zależności zostaną 
opisane w ostatnim punkcie niniejszego rozdziału. 


W naszej książce i w programie demonstracyjnym nie będziemy tworzyć 
układów lokalnych względem układu lokalnego. Położenie dowolnych obiektów 
będzie opisywane zawsze wyłącznie względem układu globalnego i lokalnego. 
W tym miejscu przypomnimy definicję struktury VERTEX przechowującej 
informację o położeniu punktów. 


typedef struct vert 


/* struktura opisująca wierzchołek */ 


double globalx;  // współrzędne globalne wierzchołka 
double globaly; 
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double 
double 


globalz; 
localx; // współrzędne lokalne wierzchołka 


double localy; 


double 


JVERTEX; 


localz; 


Pola globalx, globaly, globalz przechowują wartości trzech współrzędnych 
punktu wyrażone w układzie globalnym. Natomiast pola localx, localy, localz 
odpowiadają współrzędnym tego punktu wyrażonym względem układu lokal- 
nego. Wiemy już czym jest układ lokalny i do czego może służyć. Teraz 
pokażemy, w jaki sposób można zrealizować go jako strukturę w pamięci 
komputera, wykorzystując omówione wcześniej struktury danych i funkcje. 


Struktura AXES opisuje położenie i orientację osi układu lokalnego. 


[WE 


AXES 
, Typa 
„Fi 


typedef struct ax 


( 


char 
int 
int 
int 
int 
double 
double 


/* struktura opisująca osie układu lokalnego */ 


*name; /* nazwa nadana przez użytkownika */ 
number; /* unikalny numer */ 
type; 


n vertices;/* liczba wierzchołków */ 

n edges; /* liczba krawędzi */ 

alfa; /* kąty obrotu w mierze łukowej */ 
beta; 
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double fi; 


VERTEX huge *| vertices; /* lista wierzchołków */ 
EDGE huge  *l edges; /* lista krawędzi */ 

/* obiekt zwiazany z tym ukłedem lokalnym */ 
struct obj *object; 
struct ax *next; /* następny układ lokalny */ 


struct ax *local; /* układ lokalny ( nie używane ) */ 
JAXES; 


Ponieważ na scenie może znajdować się wiele obiektów, a z każdym z nich jest 
związany lokalny układ, struktury AXES mogą być łączone w listę. Do tego 
celu jest wykorzystywane pole *next. 


L_vertices to lista czterech wierzchołków jednoznacznie opisujących położenie 
i orientację układu lokalnego w przestrzeni. Pierwsza struktura typu VERTEX 
na tej liście określa położenie jego środka w układzie globalnym. Następne trzy 
wierzchołki reprezentują końce osi układu lokalnego (patrz rys. 3.12). Aby 
obrócić lub przesunąć układ lokalny, należy wykonać pożądaną transformację 
na wszystkich elementach tej listy. Taki sposób przechowywania informacji 
pozwala na wykorzystywanie funkcji wykonujących transformacje na wierz- 
chołkach brył do przekształceń układów lokalnych. 


Lista krawędzi: |_edges pozwala na narysowanie na ekranie osi układu lokal- 
nego. Wystarczy zmodyfikować funkcję DrawObject, aby jej argumentem była 
struktura typu AXES. 


Object — wskaźnik do struktury OBJECT opisującej bryłę znajdującą się w 
układzie lokalnym. To pole przybiera wartość NULL, do układu nie dołączono 
bryły. 


Alfa,beta,fi — kąty przejścia pomiędzy układem lokalnym a globalnym. Zostaną 
one opisane szczegółowo w następnym punkcie. 


Strukturę AXES tworzy funkcja add_axes: 


AXES *add axes( char *name, int number, 
double x _local_ scale, 
double y_ local scale, 
double z_local scale ) 
Oto przykładowe wywołanie tej funkcji: 
AXES *UkladLokalny; 
SCENE  *sceną; 
scena = add _scene( "", 0.5, 0.5, 0.5 ); 
UkladLokalny = add _axes( "UL1", 1, 40, 40, 40 ); 
add _axes_to _scene( scena , UkladLokalny ); 
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Y 


U. Lokalny AXES 


Rys. 3.12 


W ostatniej linii przykładu nowy układ lokalny jest dołączany do sceny. 


3.4.2. Przekształcenia 


Wykonywanie złożonych transformacji, a także rzutowanie obiektów na płasz- 
czyznę (metody rzutowania zostaną opisane w następnym rozdziale), powoduje 
konieczność przejścia z jednego układu odniesienia do innego. Na przykład 
obracając bryłę w układzie lokalnym, zmieniamy tylko lokalne współrzędne jej 
wierzchołków. W celu uzyskania odpowiedniej zmiany współrzędnych global- 
nych (tzn. wyrażonych względem układu globalnego) należy znaleźć prze- 
kształcenie sprowadzające układ lokalny na miejsce globalnego. Po wykonaniu 
tego przekształcenia osie obu układów powinny się pokrywać. 


Czasami istnieje potrzeba określenia, jakie będą współrzędne punktów w 
układzie lokalnym na podstawie znajomości ich wartości w układzie globalnym. 
Dzieje się tak na przykład podczas generowania obrazów z punktu widzenia 
obserwatora. 


Zmianę układu odniesienia możemy przeprowadzić wykonując trzy obroty 
i jedno przesunięcie. Na rysunku 3.13 przedstawiamy sprowadzenie układu 
lokalnego do układu globalnego wykonane w czterech krokach: 


1) Przesuwamy układ lokalny do początku układu globalnego, tak aby punkty 
[0,0,0] obu układów pokryły się (rys. 3.13a). Jest to realizowane jako translacja 
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Rys. 3.13 
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o wektor V — Trans(-vx,-vy,-vz), gdzie vx,vy,vz są współrzędnymi początku 
układu lokalnego liczonymi względem układu globalnego. 


2) Obracamy układ lokalny wokół osi Y układu głobalnego. Kąt obrotu dobie- 
ramy w ten sposób, aby rzut osi X układu lokalnego na płaszczyznę ZX w 
układzie globalnym pokrywał się z osią X układu globalnego (rys 3.13b). W 
skrócie: jest to obrót o kąt a wokół osi Y w układzie globalnym. Kąt a 
wyznaczamy z wzoru: 


a = atan2 (z,,x,) 


Funkcja atan2(y,x) wyznacza kąt z przedziału <-pi, pi>. Tangens tego kąta jest 
równy y/x (patrz rys. 3.14). 


a = atan2(y,x) 


Rys. 3.14 


Osie układu lokalnego są opisane w układzie globalnym jako wektory: 
oś X jako wektor: 


[x,,y,, ZJ] 


oś Y jako wektor: 


[Xy,Yy Z] 
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oś Z jako wektor: 


l x z” y z) Z 1 
Współrzędne tych wektorów występują w naszych równaniach. 


3) Obracamy układ lokalny wokół osi Z układu globalnego w ten sposób, aby 
oś Xukładu lokalnego pokryła się z osią X układu globalnego (rys. 3.13c). Obrót 
jest wykonywany o kątB: 


B= — atan2 (y, ,Yx? + zę ) 
4) Obracamy układ lokalny wokół osi X układu globalnego o kąt g (rys. 3.13d): 


p = — atan 2 (z „cosa — x, sina, z,sina sin + 
xycos a sin B + y, cos  ) 


Po tym obrocie układy lokalny i globalny pokrywają się (rys. 3.13e). W zapisie 
macierzowym złożenie czterech przedstawionych transformacji wygląda tak: 
Z = Rot (x,6) Rot (z, 8) Rot (y, a) Trans (= v „„,—v 


yo TVz ) 


Wynikiem powyższego iloczynu jest macierz przejścia z układu globalnego do 
układu lokalnego. Jeżeli punkt P ma w układzie globalnym współrzędne 
[z,y,z,1], to w układzie lokalnym jego współrzędne będą równe: 


[x,y,z,1]=Z e [x,y,z,1] 


Taki zapis oznacza, że współrzędne punktu P najpierw mnożymy przez macierz 
translacji, potem kolejno przez macierze obrotu wokół osi Y, Z, X. Oczywiście 
nie musimy wyznaczać bezpośrednio macierzy Z. Możemy wykonywać prze- 
kształcenia kolejno, jedno po drugim. W taki sposób działa funkcja From- 
GlobalCoordsToLocal: 


// Funkcja przeliczająca współrzędne globalne na lokalne. 
// Znamy współrzędne globalne wierzchołków 

// na liście wskazywanej przez *vert, 

// a także położenie i orientację osi układu lokalnego 
// (wyrażone we współrzędnych układu globalnego). 

// Na ich podstawie obliczamy i zapamiętujemy 

// współrzędne lokalne wierzchołków. 


void FromGlobalCoordsToLocal ( AXES *axeS, 
VERTEX huge *vert ) 
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( 


double xe,ye,ze; 


// zapamiętanie położenia początku układu lokalnego 
xe=axes->l| edges->vertl->globalx; 
ye=axes->l edges->vertl->globaly; 
ze=axes->|l edges->vertl->globalz; 


CopyGlobalCoordsToLocal( vert ); 

// transformacje 

move( vert , LOCAL , -xe , -ve , -ze ); 
rot_y( vert , LOCAL , axes->alfa ); 
rot _z( vert , LOCAL , -axes->beta ); 
rot _x( vert , LOCAL , -axes->fi ); 


Jeżeli ostatnie równanie pomnożymy obustronnie przez macierz odwrotną do 
macierzy Z, to otrzymamy równanie opisujące przejście z układu lokalnego do 
globalnego, czyli operację odwrotną do opisanej powyżej: 


[x,y,z]=[x,y,z']*Z" 


Z! = Trans (v, , v,, v, ) Rot (y ,—a) Rot (z ,-8) Rot (x,-$) 


Znając współrzędne punktu P[x,y,z',1] w układzie lokalnym, możemy na 
podstawie tego równania znaleźć współrzędne x,y,z w układzie globalnym. 
Należy wykonać te same przekształcenia co poprzednio, ale w odwrotnej 
kolejności i z przeciwnymi argumentami. Oto funkcja realizująca to zadanie: 


/7/ 
/7 
/7/ 
/7/ 
/7/ 
/7/ 
// 


Funkcja przeliczająca współrzędne lokalne na globalne. 
Znamy współrzędne lokalne wierzchołków 

na liście wskazywanej przez *vert, 

a także położenie i orientację osi układu lokalnego 

( wyrażone we współrzędnych układu globalnego ). 

Na ich podstawie obliczamy i zapamiętujemy 
współrzędne globalne tych wierzchołków. 


void FromLocalCoordsToGlobal( AXES *axes , 


( 


VERTEX huge *vert ) 


double xe,ye,ze; 


// zapamiętanie położenia początku układu lokalnego 
xe=axes->l edges->vertl->globalx; 


Realistyczna grafika 3D 85 


ye=axes->l edges->vertl->globaly; 
ze=axes->l_edges->vertl->globalz; 


cCopyLocalCoordsToGlobal( vert ); 

// transformacje 

rot_x( vert , GLOBAL , axes->f1 ); 
rot_z( vert , GLOBAL , axes->beta ); 
rot_y( vert , GLOBAL , -axes->alfa ); 
move( vert , GLOBAL , xe , ye , ze ); 


W obu funkcjach układ lokalny jest zdefiniowany w strukturze typu AXES, do 
kórej wskaźnikiem jest *axes. Zauważmy, że wartości kątów obrotu w fun- 
kcjach FromLocalCoordsToGlobal i FromGlobalCoordsToLocal są pobierane z 
odpowiednich pól struktury typu AXES opisującej układ lokalny. Aby można 
było skorzystać z dowolnej z tych funkcji, należy uprzednio wywołać funkcję 
CalculateAlfaBetaFi obliczającą te kąty według podanych przez nas wzorów. 


// Obliczenie kątów przejścia z układu lokalnego 

// do globalnego. 

// Kąty alfa, beta, fi są liczone w stopniach 

// i zapisywane do struktury wskazywanej przez axes. 


void CalculateAlfaBetaFi( AXES *axeS ) 
( 

VERTEX huge *vert; 

EDGE huge *edge; 

double xe,ye,ze; 

double xx,yxX,ZX; 

/* składowe wersora x układu lokalnego */ 

double xy,yy,Zy; 

double alfa,beta,fi,s; 

double ca,sa,cb,sb,cf,sf; 

double x,y,z,xl,yl,zl; 


edge=axes->l_edges; /* to osie układu lokalnego */ 
vert=edge->vertl; 

xe=vert->globalx; 

ye=vert->globaly; 

ze=vert->globalz; 

edge=edge->next; /* pomijamy oś Z */ 
vert=edge->vert2; 
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/* oś X układu lokalnego we współrzędnych +*/ 

/* układu globalnego */ 

xx=vert-=>globalx-xe; 

yx=vert->globaly-yve; 

zx=vert->globalz-ze; 

edge=edge->next; 

vert=edge->vert2; 

/* oś Y układu lokalnego we współrzędnych */ 

/* układu globalnego */ 

xy=svert->globalx-xe; 

yy=vert->globaly-ye; 

zy=vert->globalz-ze; 

/* to jeszcze nie są kąty przejścia z układu lokalnego 
do globalnego */ 

alfa=atan2( zx , xx ); 

s=sqrt( Xxx*xXxtzx*ZzX ); 

beta=atan2( yx , s ); 


ca=cos( alfa ); 
sa=sin( alfa ); 
cb=cos( -beta ); 
/* -beta to kąt obrotu wokół osi z ukł. glob. */ 
sb=sin( -beta ); 
fi=atan2( zy*ca-xy*sa , zy*sa*tsbtxy*ca*sbtyytcb ); 
axes->alfa=alfa*180/M PI; 
axes->beta=beta*180/M PI; 
axes->fi=fi*180/M PI; 
) 


Przedstawiony powyżej sposób zmiany układu współrzędnych jest wykorzysty- 
wany także w programie demonstracyjnym. Istnieją jednak inne metody. 
Polegają one na obrotach wokół inaczej wybranych osi lub układów odniesienia. 
Najbardziej rozpowszechnioną z tych metod jest metoda kątów Eulera. 


Kąty Eulera są wyznaczane pomiędzy różnymi położeniami obracanego układu 
współrzędnych. Zakładając, że początki układów znajdują się w tym samym 
miejscu, musimy wykonać trzy transformacje. Obracamy układ lokalny wokół 
jego własnej osi Z. Położenie osi X i Y układu lokalnego zmienia się. Otrzymany 
układ obracamy wokół jego nowej osi X, a następnie ponownie wokół osi Z. 


Takie postępowanie bywa czasem nazywane obrotami względem ruchomego 
układu odniesienia. Rozwiązanie przyjęte przez nas można nazwać obrotami 
względem nieruchomego układu odniesienia, ponieważ kolejne obroty były 
wykonywane wokół różnych osi nieruchomego układu globalnego. 
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ROZDZIAŁ 4 


Rzutowanie 


Rzutowanie jest podstawowym przekształceniem w grafice komputerowej. 
Wizualizacja jakiegokolwiek obiektu przestrzennego wymaga odwzorowania 
go na płaszczyznę. Dopiero mając dane dwie współrzędne punktu (zamiast 
trzech) 


możemy narysować odcinek łączący dwa punkty na ekranie komputerowym. 
Metody rzutowania dzielimy na dwie klasy: 


m rzutowanie równoległe, 


M rzutowanie perspektywiczne (środkowe). 


4.1. Rzutowanie równoległe 


Rzutowanie równoległe zachowuje proporcje rzutowanych przedmiotów. Od- 
cinki równoległe przed rzutowaniem są także równoległe po rzutowaniu. 
Zachowane są stosunki odległości między punktami. Przykładem zastosowania 
tej metody jest rysunek techniczny. Przedstawia się na nim trzy widoki obiektu: 
widok od frontu, z góry iz boku. Każdy z tych widoków jest rzutem równoległym 
na jedną z trzech prostopadłych do siebie płaszczyzn. 


W programie demonstracyjnym rzutowanie równoległe zastosowaliśmy do 
prezentacji trzech rzutów sceny na płaszczyzny X-Y, Ń-Z, Y-Z wyznaczone przez 
osie układu globalnego. Jeżeli płaszczyzną rzutowania jest X-Y, to po prostu 
do wyświetlania punktów używamy tylko współrzędnych x i y,'a współrzędną 
z pomijamy. Taka forma wizualizacji pozwala użytkownikowi na jednoznaczne 
określanie położenia dowolnego punktu w przestrzeni. W podrozdziale 2.3.4 
przedstawiliśmy funkcję DrawObject działającą w opisany wyżej sposób. 


Niestety, ten sposób rzutowania nie pozwala na realistyczne przedstawienie 
"trójwymiarowości" obiektów. Oglądając wizerunek sześcianu na rysunku 4.la, 
nie jesteśmy w stanie stwierdzić, który jego bok leży bliżej nas. Rysunki 4.1bi 
4.1c przedstawiają dwie możliwe interpretacje rysunku 4.la. 
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Rys. 4.1 


4.2. Rzutowanie perspektywiczne 


Aby bardziej realistycznie przedstawić obiekt trójwymiarowy, należy posłużyć 
się rzutem perspektywicznym (środkowym). Ta forma prezentacji stwarza u 
odbiorcy wrażenie przestrzenności oglądanego obrazu dwuwymiarowego. Rzu- 
ty obiektów leżących dalej od płaszczyzny rzutowania są mniejsze niż rzuty 
takich samych obiektów leżących bliżej tej płaszczyzny. Jest to symulacja 
zjawiska perspektywy występującego w Świecie rzeczywistym; proste równo- 
ległe zbiegają się na horyzoncie, a przedmioty maleją w miarę oddalania się od 
nich. 


Na rysunku 4.2a przedstawiono obraz sześcianu z rysunku 4.1la po rzutowaniu 
perspektywicznym. Teraz większość odbiorców nie będzie miała wątpliwości, 
że oglądany rysunek przedstawia szkielet bryły z rys. 4.2b. Oczywiście, jak 
zwykle w takich przypadkach bywa, znajdą się oponenci twierdzący, że możli- 
wa jest także interpretacja taka, jak na rys. 4.2c. Przyjmijmy, że takie twier- 
dzenia są tylko wyjątkiem potwierdzającym regułę. 


a) b) c) 
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Rys. 4.3 


Najprostszą metodą rzutowania perspektywicznego jest umieszczenie obser- 
watora na osi Z globalnego układu odniesienia, w punkcie E o współrzędnych 
[0,0,-d]. Płaszczyzną rzutowania (rzutnią) jest w tym przypadku płaszczyzna 
rozpięta na osiach X i Y. Rzutowanie dowolnego punktu P[x,y,z] polega na 
znalezieniu współrzędnych x i y punktu P”. P' jest punktem przecięcia płasz- 
czyzny X-Y przez prostą poprowadzoną od położenia obserwatora [0,0,-d] do 
rzutowanego punktu. Na rysunku 4.3 jest przedstawiony sposób rzutowania. 
Współrzędne punktu P'” wyznaczamy z następujących zależności: 


XxX '=x d 
z+d 

= d 
4 "Jz+d 


Powyższe równania wynikają z podobieństwa trójkątów EOP' i EQP. Ponieważ 
te trójkąty są podobne, także ich rzuty równoległe na płaszczyznach X-Z i Y-Z 
muszą być do siebie podobne. Równania otrzymujemy porównując stosunki 
długości boków w rzutach obu trójkątów na wspomniane płaszczyzny. 

Jeżeli punkt, z którego chcemy oglądać interesujące nas obiekty, nie leży na 
osi Z układu globalnego, to należy obrócić cały układ współrzędnych w taki 
sposób, aby obserwator znalazł się w [0,0,-d]. Pisząc o obrocie układu współ- 
rzędnych, mamy na myśli oczywiście obrócenie wszystkich obiektów znajdują- 
cych się w tym układzie. Takie przekształcenie jest równoważne opisywanemu 
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wcześniej przejściu z układu globalnego do układu lokalnego, przy czym układ 
lokalny to ten, w którym obserwator jest położony w punkcie [0,0,-d]. 


4.3. Ostrosłup widzenia 


Opisane w poprzednim punkcie rozwiązanie zagadnienia rzutowania zawodzi, 
gdy wprowadzamy obserwatora pomiędzy obiekty znajdujące się na scenie. Na 
płaszczyznę rzutowania będą rzutowane także te obiekty, które znajdują się za 
obserwatorem i nie powinny być widoczne. Problem ten może być rozwiązany 


poprzez zdefiniowanie ostrosłupa widzenia i rzutowanie tylko tych obiektów, 
które się w nim znajdują. 


Ostrosłup widzenia jest matematyczną konstrukcją pozwalającą komputerowi 
na symulowanie sposobu działania oka bądź aparatu fotograficznego. Na 
początku objaśnijmy, w jaki sposób powstają obrazy w aparacie fotograficznym. 
Każdy przedmiot może zostać sfotografowany, ponieważ promienie odbite od 
niego rozchodzą się w przestrzeni. Niektóre z nich trafiają do obiektywu 
aparatu. W soczewkach promienie ulegają załamaniu, po czym biegną dalej aż 
padną na błonę fotograficzną. Tu powstaje dwuwymiarowy obraz stożkowego 
wycinka przestrzeni, z którego przybyły promienie. Mamy do czynienia ze 
stożkiem, ponieważ przedłużenia wszystkich promieni, które przeszły przez 
soczewki obiektywu, zbiegają się w jednym punkcie (nazywanym ogniskiem 
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obiektywu). Taki opis jest oczywiście bardzo uproszczony, jednak zastosowanie 
podobnej metody w grafice komputerowej przynosi wyśmienite rezultaty. 


W naszej książce obiekt, którego zadaniem jest symulowanie działania aparatu 
fotograficznego, nazwaliśmy obserwatorem. Rysunek 4.4 wyjaśnia, w jaki 
sposób jest skonstruowany obserwator i jakie zastosowanie ma w tym przypad- 
ku pojęcie ostrosłupa widzenia. 


Z obserwatorem związany jest lokalny układ współrzędnych rozpięty na we- 
ktorach u, w, v. Płaszczyzna rzutowania jest wyznaczona przez wektory u i w. 
Możemy porównać ją do klatki filmu o rozmiarach 2u na Źw, znajdującego się 
w wylmaginowanym aparacie fotograficznym, skierowanym przez obserwatora 
zgodnie z kierunkiem wektora v. 


Prowadząc cztery półproste od punktu O, w którym znajduje się obserwator, 
przez punkty A, B, C, D, otrzymujemy ostrosłup widzenia. Ostrosłup widzenia 
pełni w tym przypadku rolę analogiczną do pojęcia stożka widzenia aparatu 
fotograficznego. Jeżeli punkt znajduje się wewnątrz tego ostrosłupa, to rzut 
tego punktu na płaszczyznę rzutowania leży wewnątrz prostokąta ABCD. W 
przeciwnym wypadku rzut punktu leży poza prostokątem, czyli jest niewido- 
czny dla obserwatora. 


Zakładamy, że długość wektora v jest stała i równa 1. Wtedy poprzez modyfi- 
kowanie długości wektorów u i w możemy zmieniać rodzaj obiektywu w naszym 
"aparacie". Wydłużając wektory u i w powiększamy kąt rozwarcia ostrosłupa, 
w efekcie czego poszerzamy pole obserwacji. W klasycznej fotografii odpowiada 
to zastosowaniu obiektywu szerokokątnego. Postępując w ten sposób można 
osiągnąć kąt widzenia równy niemal 180 stopni, jednak powstają przy tym 
znaczne deformacje obrazu. Zawężenie pola obserwacji można osiągnąć po- 
przez zmniejszanie długości tych wektorów. Odpowiada to zastosowaniu teleo- 
biektywu do fotografowania odległych przedmiotów. 


Opisane powyżej operacje można wykonać, stosując transformację skalowania 
osi układu lokalnego obserwatora. Zauważmy, że w tym wypadku skalowanie 
dotyczy nie obiektu znajdującego się w układzie lokalnym, ale osi tego układu. 


Jak już wspominaliśmy, obserwator jest "przyczepiony" do wyróżnionego ukła- 
du lokalnego. Wykonując transformacje na tym układzie, wykonujemy je 
równocześnie na obserwatorze. Wynika stąd, że istnieje możliwość poruszania 
się i obracania obserwatora w obrębie sceny złożonej z wielu obiektów. 


W tym miejscu należy odpowiedzieć na pytanie: jak przeprowadzić rzutowanie, 
jeśli obserwator nie znajduje się na osi z i nie jest skierowany na początek 
układu współrzędnych? Istnieją co najmniej dwa rozwiązania tego problemu. 
Po pierwsze, wiemy, że istnieje możliwość zmiany współrzędnych z układu 
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globalnego na współrzędne układu lokalnego, w tym przypadku układu obser- 
watora. Można więc znaleźć współrzędne wierzchołków wszystkich obiektów 
w układzie obserwatora i zastosować opisane wcześniej rzutowanie perspekty- 
wiczne, operując na tych współrzędnych. 
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Rys. 4.5 


Po drugie, na podstawie położenia osi układu obserwatora można analitycznie 
wyznaczyć współrzędne punktów po rzutowaniu, sprawdzając jednocześnie czy 
te punktu są widoczne. Sposób rozwiązania tego zagadnienia został przedsta- 
wiony w [1]. My przedstawimy wyłącznie wzory. 


Punkt P[x ,y,z'] jest obrazem punktu P[x,y,z] w przekształceniu polegającym 


na rzutowaniu P na płaszczyznę rozpiętą na wektorach u,w. Współrzędne tego 
punktu są dane wzorami: 


EP o — 
s=—__lul 
EPOv 

EP o —* 

_ lwl 
J EPov 


gdzie: 


A r A 
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EP — wektor o początku w punkcie E i końcu w 


P, 
|u| — długość wektora u, 
|w| — długość wektora w. 


Powyższe wzory nie są określone dla EP © v = 0. Punkt P, dla którego zachodzi 
EP O v = 0, leży na płaszczyźnie przechodzącej przez środek rzutowania E 
1 równoległej do rzutni. Punkty, dla których EP 0 v < 0, to znaczy kąt między 
wektorami EP i v jest większy od 90 stopni, leżą za obserwatorem i nie są przez 
niego widziane. Jeżeli x',y' spełniają następujące warunki: 


-|u| <X'< ul, 
-|w| <y'< KIE 


to punkt P leży wewnątrz ostrosłupa widzenia i jest widoczny. 


Struktura typu OBSERVER służy do przechowywania informacji o położeniu 
obserwatora i rozmiarach ostrosłupa widzenia. Jej najważniejszym komponen- 
tem jest struktura AXES, opisująca układ lokalny powiązany z obserwatorem. 
Jest to zwyczajny układ lokalny z odpowiednio dobranymi długościami osi. Oś 
X reprezentuje w tym przypadku wektor u, a oś Y wektor w. Kierunek osi Z 
wyznacza kierunek patrzenia obserwatora, czyli pełni rolę wektora v (patrz 
rys. 4.4 i 4.5). 


Pola modu i modw przechowują długość wektorów u i w, a co za tym idzie 
informację o tym, jakiego typu obiektyw używamy. 


Wskaźnik *target wskazuje na układ lokalny, w kierunku którego jest skiero- 
wany obserwator. Obserwatorowi można przyporządkować układ lokalny tar- 
get, wtedy niezależnie od zmian wzajemnego położenia obserwatora 1 tego 
układu, obserwator będzie zawsze "patrzył" na początek tego układu. Stosowne 
transformacje układu lokalnego obserwatora wykonuje funkcja LookAtTarget. 
Kod źródłowy tej funkcji wraz z krótkim opisem parametrów znajduje się w 
pliku mat.c. 


typedef struct obs // obserwator 
( 
double modu; 
// obiektyw szerokokątny gdy modu,modw - duże (np.>l) 
double  modw; 
// teleobiektyw gdy modu,modw - małe 
AXES *position; 
// obiekt opisujący położenie i kierunek "patrzenia" 
obserwatora */ 
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AXES *target; 


// układ lokalny, na który "patrzy" obserwator 
) OBSERVER; 


Na scenie może znajdować się tylko jeden obserwator. Do utworzenia struktury 
OBSERVER służy funkcja add_observer. Parametry x, y, z określają położenie 
obserwatora, a modu i modv długości wektorów u i w. Najbardziej realistyczną 
perspektywę uzyskujemy dla: 


0.3 < modu < 0.5 
0.3 < modw < 0.5 


Jeżeli użyjemy większych wartości, to otrzymamy obrazy zdeformowane. Jeżeli 
użyjemy mniejszych wartości, to uzyskamy duże powiększenie i małą perspe- 
ktywę. Funkcja add_observer standardowo ukierunkowuje osie układu obser- 


watora, tak aby "patrzył" w kierunku centrum układu globalnego (wywołanie 
funkcji LookAtTarget). 


Jeżeli zabrakło pamięci, funkcja zwraca NULL — w ten sposób możemy 
kontrolować poprawność jej działania. 


OBSERVER *add _observer( double x, double y, double z, 
double modu , double modw ) 


( 
OBSERVER *observer; 


if( ( observer=new observer() ) == NULL ) 
return NULL; 

observer->target=NULL; 

observer->modu=modu; 

observer->modw=modw; 


if( ( observer->position=add axes( NULL, 0, modu, 
modw, 1 ) ) == NULL ) 

( 

delete_observer( observer ); 

return NULL ; 
) 

// ustalenie położenia i kierunku patrzenia obserwatora 
move (  observer->position->l vertices, 
GLOBAL, Xx, y, Z, ); 

LookAtTarget( observer , NULL , O 
return( observer ); 


0,0 ); 


) 
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Po utworzeniu obserwatora należy dołączyć otrzymaną strukturę do sceny za 
pomocą funkcji add_observer_to_scene (kod źródłowy w pliku listy.c). 


Zadanie rzutowania perspektywicznego realizuje funkcja observer_view: 


void observer view ( OBJECT *object, OBSERVER *observer, 
int x res, int y res, double xaspect, 
double yaspect ) 


double ex,ey,ez; 

double xep,yep,zep; 

double xu,yu,zu,xw,yw,Zw,XV,YyV,ZV; 
double vz; 

double modu ,modw; 

int x _res2,y res2; 

double sx,sy; 

double epv; 

VERTEX huge *vert; 

EDGE huge  *edge; 


edge=observer->position->l_edges; 
vert=edge->vertl; 

ex=vert->globalx; // położenie obserwatora 
ey=vert->globaly; 

ez=vert->globalz; 

vert=edge->vert2; 

xv=vert->globalx-ex; // współrzędne wektora v 
yv=vert->globaly=ey; 

zv=vert->globalz-ez; 

edge=edge->next; 

vert=edge->vert2; 

xu=vert->globalx-ex; // współrzędne wektora u 
yu=vert->globaly-ey; 

zu=vert->globalz-ez; 

edge=edge->next; 

vert=edge->vert2; 

xw=vert->globalx-ex; // współrzędne wektora w 
yw=vert->globaly-ey; 

zw=vert->globalz-ez; 


modu=observer->modu; 
modw=observer->modw; 
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x _res2=x_res/2; 
y_res2=y res/2; 
8x=X _res/(2*modu); 
sy=y_res/(2*modw); 
// aby nie zniekształcić obrazu w przypadku, gdy sx <> sy 
if( 8x > Sy ) 
SX=SYy; 
else 
Sy=SX; 
// jeśli piksle nie są kwadratowe, to wynikające z tego 
// zniekształcenie kompensujemy dzięki xaspect i yaspect; 
sx*=xaspect; 
sy*=yvaspect; 


// rzutowanie wszystkich wierzchołków bryły 

vert=object->l_vertices; 

while( vert != NULL ) 

( 
xep=vert-=>globalx-ex; 

// współrzędne wektora łączącego obserwatora z punktem 

yep=vert->globaly-ey; 
zep=vert->globalz=ez; 


epv=xep*xvtyep*yvtzep*zv; // iloczyn EP o V 
if( epv <= 0 ) 
vert->info=INVISIBLE; 
// rzut wierzchołka nie jest widoczny 
else 
( 
vert->info=VISIBLE; 
// rzut wierzchołka jest widoczny 
vert->drawx=x _res2+sx*( xep*xutyep*yutzep*zu )/ 
( epv*modu ); 
vert->drawy=y res2-sy*( xep*xwtyeptywtzep*zw )/ 
( epv*modw ); 
) 


vert=vert=>next; 


) 


Rzutowane są wszystkie wierzchołki, z których jest zbudowana bryła, opisy* 
wana przez strukturę OBJECT. Parametr *observer wskazuje na strukturę 
typu OBSERVER, definiującą obserwatora. 
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X_res, y_res — szerokość i wysokość ekranu (w pikslach), na którym będziemy 
wyświetlać rzut bryły. Po rzutowaniu rzuty wszystkich widocznych wierzchoł- 
ków znajdują się w obrębie prostokąta ABCD (patrz rys. 4.4). Ponieważ 
rozmiary tego prostokąta najczęściej nie pokrywają się z rozmiarami ekranu 
(wysokość i szerokość ekranu są znacznie większe niż wysokość i szerokość 
prostokąta ABCD), należy dodatkowo przeskalować współrzędne rzutów. 
Właśnie do tego celu są używane parametry x_res i y_res. 


Xaspect, yaspect określają proporcje piksla ekranu, na którym będziemy 
wyświetlać rzut bryły. W zależności od trybu graficznego, w jakim pracuje 
karta video w naszym komputerze, zmieniają się proporcje i rozmiary piksla. 
Jak wiemy, piksel jest najmniejszym elementem obrazu. Każdy obraz rastrowy 
składa się z tablicy piksli. Kształt rysowanych figur zależy więc od proporcji 
boków pojedynczego piksla. Jeżeli narysujemy kwadrat o boku 100 na 100 
piksli na ekranie, na którym stosunek wysokości piksla do jego szerokości ma 


Rys. 4.6 
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sięjak 1 do 1, tootrzymamy kwadrat o równych bokach (rys. 4.6a). Jeżeli jednak 
wysokość piksla będzie dwa razy większa niż szerokość, to otrzymamy prosto- 
kąt — rysunek 4.6b. Rysunek 4.6c ilustruje przypadek, w którym szerokość 
piksla jest większa niż wysokość. 


Parametry xaspect i yaspect służą do kompensacji takich zniekształceń obrazu. 


Funkcja observer_view oblicza współrzędne rzutów na podstawie podanych 
wcześniej wzorów. Wyniki obliczeń, czyli współrzędne x i y rzutów, są zapisy- 
wane do pól drawx i drawy w strukturze VERTEX, opisującej rzutowany 
wierzchołek. Jeśli rzut danego wierzchołka jest widoczny, to pole info otrzyma 


wartość VISIBLE. 
typedef struct vert /* struktura opisująca wierzchołek */ 
( 
int info; 
double drawx; // współrzędne wierzchołka po rzutowaniu 
double drawy; 
JVERTEX; 


Funkcja observer_view_scene wykonuje rzutowanie dla wszystkich obiektów 
znajdujących się na scenie: 


void observer _view scene( SCENE *scene, int x res, 
int y res, double xaspect, 
double yaspect ) 


AXES *axeS; 
axes=scene->l_axe3; 
while( axes != NULL ) 
t 

if( axes->object !l= NULL ) 

observer_view( axes->object, scene->observer, 
x res, y_res, xaspect, yaspect ); 
axes=axes->nexXt; 


) 


Poznaliśmy już struktury danych opisujące podstawowe obiekty. Na rysunku 
4.7 przedstawiamy przykładowe połączenia między tymi strukturami. Jest to 
reprezentacja sceny, na której znajduje się obserwator, dwa źródła światła 
(struktura LAMP zostanie opisana w punkcie 5.4) i trzy układy lokalne. W dwu 
układach lokalnych są umieszczone bryły. 


hh ..........dÓóÓdÓŚóŚó a — 


100 M. Domaradzki, R. Gembara 


Rys. 4.7 


Jak łatwo zauważyć, mając dany wskaźnik do struktury typu SCENE, możemy 
dotrzeć do informacji o dowolnym obiekcie. Taka organizacja danych umożliwia 
także tworzenie wielu niezależnych scen w pamięci komputera i wyświetlanie 
obiektów znajdujących się tylko na jednej z nich. 


Do tworzenia struktury, takiej jak na rysunku 4.6, służy zestaw czterech 
funkcji umożliwiających łączenie różnych obiektów: 


- observer_to_scene( SCENE *sc , OBSERVER *observer ) 
m Dołączenie obserwatora do sceny. 

- add_axes_to_scene( SCENE *scene , AXES *axes ) 
m Dołączenie układu lokalnego AXES do sceny. 

- add_lamp_to_scene( SCENE *scene , LAMP *lamp ) 


m Dołączenie lampy do sceny. 
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- add_object to axes( SCENE *scene , AXES *axes , 
OBJECT *object ) 


M Dołączenie bryły (object) do układu lokalnego (axes) znajdującego się na 
scenie. 


4.4. Przykłady 


W tym podrozdziale postaramy się pokazać, w jaki sposób utworzyć prostą 
scenę, jak można manipulować obserwatorem i jak wyświetlić obrazy perspe- 
ktywiczne na ekranie komputera. 


Na początku przedstawimy funkcję służącą do rysowania rzutów brył. Będzie 
to modyfikacja używanej przez nas wcześniej funkcji DrawObserverViewOb- 
ject. Zmiany polegają na tym, że współrzędne punktów pobieramy z pól drawx, 
drawy struktur VERTEX. Przedstawimy także funkcję DrawScene rysującą 
rzuty wszystkich brył znajdujących się na scenie. Przed użyciem którejkolwiek 
z tych funkcji należy wywołać observer_view, gdy rysujemy rzut pojedynczej 
bryły lub observer_view_scene, gdy rysujemy rzuty wszystkich brył znajdują- 
cych się na scenie. Nie istnieje potrzeba wcześniejszego wywoływania funkcji 
LogicalToScreen. 


// rysowanie rzutu perspektywicznego pojedynczej bryły 
void DrawObserverViewobject( OBJECT *object ) 
t 
VERTEX huge  *vertl, 
huge  *vert2; 
EDGE huge  *edge; 


if ( object == NULL ) 
return; 
/* wskaźnik edge przesuwamy po liście krawędzi. 
Każda struktura typu EDGE znajdująca się na tej liście 
zawiera wskaźniki do struktur VERTEX, opisujących 
współrzędne wierzchołków. Z pól drawx, drawy tych 
struktur pobieramy współrzędne x i y danego wierzchołka 
i rysujemy linię */ 
edge=object->l_edges; 
while( edge != NULL ) 
( 
vertl=edge->vertl; 
vert2=edge->vert2; 
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line ( vertl->drawx, vertl->drawy, vert2->drawx, 
vert2->drawy ); 
edge=edge->next; 


) 


// Rysowanie rzutów wszystkich brył 
// znajdujacych się na scenie 
void DrawScene( SCENE *scene ) 


( 
AXES  *axe3; 


axes = scene->l _ axes; 

if ( axes == NULL ) 
return; 

while( axes != NULL ) 

( 
if( axes->object 1= NULL ) 

DrawobserverViewobject( axes->object ); 

axes=axes->next; 


Przykład 1 


W rozdziale 2 zamieściliśmy program, który tworzył sześcian i wyświetlał go 
na ekranie. Ponieważ użyliśmy rzutu równoległego, otrzymany obraz był 
pozbawiony perspektywy. Poniżej prezentujemy udoskonaloną wersją tego 
programu. Spróbujmy prześledzić, w jaki sposób są tworzone i łączone struktu- 
ry opisujące poszczególne elementy trójwymiarowej sceny. 


A. Po zainicjowaniu grafiki budujemy strukturę typu SCENE. 
B. Tworzymy strukturę AXES i dołączamy ją do SCENE. 
C. Tworzymy strukturę typu OBJECT i dołączamy ją do AXES. 


D. Tworzymy strukturę OBSERVER, opisującą aparat fotograficzny umiejsco- 
wiony w punkcie o współrzędnych 100,100,200 i także dołączamy ją do sceny. 


E. Wykonujemy rzutowanie środkowe dla wszystkich brył na scenie (w tej 
chwili na scenie jest tylko sześcian). 


F. Rysujemy rzuty brył. 
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Przedstawiony powyżej sposób postępowania może posłużyć jako wzór dla 
programów generujących bardziej skomplikowane sceny — patrz następny 
przykład. Punkty A, D, E, F powinny wystąpić w każdym takim programie w 
nie zmienionej formie. 


W miejscach programu odpowiadających punktom B i C można tworzyć wiele 
różnych brył i układów lokalnych, a także wykonywać na nich transformacje. 
Istnieje także możliwość tworzenia interakcyjnych programów, pozwalających 
na przykład użytkownikowi na przemieszczanie się po scenie, obracanie lub 
przesuwanie brył itp. w czasie rzeczywistym. W tym przypadku program 
powinien obliczać na bieżąco wszystkie zmiany współrzędnych dowolnych 
obiektów po każdej "interwencji" użytkownika (np. naciśnięciu klawisza na 
klawiaturze lub poruszeniu myszką) i rysować obraz sceny po zmianie. Prosty 
program tego typu zaprezentujemy w podrozdziale 5.6. 


// prog5.c 

tinclude "ray.h" 
tinclude "listy.c" 
include "mat.c" 
tinclude "obiekty.c" 
t$include "rysuj.c" 


void main( void ) 
t 
OBJECT *szescian; 
AXES  *UkladLokalny; 
SCENE  *sceną; 
OBSERVER *kamera; 
FEATURE  *f; 
int gdriver = DETECT, gmode, errorcode; 


// inicjalizacja grafiki 
initgraph(Ssgdriver, Ggmode, "c:Niborlandciibgi"); 


// parametry powierzchni 
f = new feature( 0, 0, 0, 1, 1, 1, 0, 1); 
if( f == NULL ) 
exit(0); 


// utworzenie sceny, układu lokalnego 

// i dołączenie układu do sceny 
scena = add_scene( "", .5, .5, .5 ); 
UkladLokalny = add _axes( "", 0, 40, 40, 40 ); 
add axes_to_scene( scena , UkladLokalny ); 
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// utworzenie struktur opisujących sześcian, 

// dołączenie ich do struktury UkladLokalny 
szescian = add cube( NULL, 0, 100, 100, 100, £f, OFF ); 
add _object to _axes( scena, UkladLokalny, szescian ); 


// utworzenie obserwatora i dołączenie go do sceny 
kamera = add observer( 100, 100, 200, .4, .4 ); 
add observer to scene( scena, kamera ); 


// obliczenie współrzędnych wierzchołków sześcianu 
// po rzutowaniu perspektywicznym i narysowanie 
// rzutu bryły 
observer_view( szescian, kamera, 640, 480, 1, 1 ); 
DrawobserverViewobject( szescian ); 


getch(); 

clearviewport(); 

// przesuwamy obserwatora 

move( kamera->position->l vertices, GLOBAL,-80,-50,100); 
LookAtTarget( kamera, NULL, 0, 0, 0 ); 

observer _ view scene( scena, 640, 480, 1, 1 ); 
DrawScene( scena ); 


Rys. 4.8a 
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Rys. 4.8b 


// zamknięcie grafiki 
getch(); 

closegraph (); 

delete scene( scena ); 


) 


Program z przykładu 1 generuje obrazek taki, jak na rysunku 4.8a, czeka na 
naciśnięcie klawisza i pokazuje widok tego samego sześcianu z innego punktu 
widzenia — rys 4.8b. 


Przykład 2 


Ten przykład pokazuje, w jaki sposób można dołączyć do sceny wiele brył. 
Program generuje 20 prostopadłościanów i układa je w kształcie spiralnych 
schodów. W przykładzie 3 z rozdziału 3 budowaliśmy taką samą figurę, ale nie 
tworzyliśmy układów lokalnych, sceny i obserwatora. Program generuje 3 
widoki schodów oglądane z trzech różnych pozycji (rys. 4.9a, 4.9b, 4.9c). 


Na rysunku 4.9a obserwator znajduje się w punkcie [0,450,-300] i "patrzy' w 
kierunku początku układu współrzędnych. Przed wykonaniem rysunku 4.9b 
obserwator został przesunięty o wektor [100,170,250]. Na rysunku 4.9c obser- 
wator jest w tym samym miejscu co poprzednio, a schody zostały obrócone 0 
kąt 45 stopni wokół osi Z. 
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// progó.c 

$include *ray.h" 
$include "listy.c" 
$include "mat.c" 
finclude "obiekty.c" 
$finclude "rysuj.c" 
fdefine LICZBA STOPNI 20 
łdefine KAT_OBROTU 15 


void main( void ) 
( 
int i; 
OBJECT *stopnie; 
AXES  *ULok; 
SCENE  *tsceną; 
OBSERVER *kamera; 
FEATURE *f; 
int gdriver = DETECT, gmode, errorcode; 


// inicjalizacja grafiki 
initgraph(sgdriver, śgmode, "c:Niborlandciibgi"); 


// parametry powierzchni 
£f = new feature( 0, 0, 0, 1, 1, 1, 0, 1); 
scena = add _scene( "", .5, .5, .5 ); 
kamera = add observer( 0, 450, -300, .4, .4 ); 
add_observer_to _scene( scena, kamera ); 


// utworzenie stopni 
for( i=0; i<LICZBA STOPNI; 1++ ) 
l 
ULOK = add _axes( "", 0, 40, 40, 40 ); 
add_axes_to_scene( scena , ULok ); 
// utworzenie struktur opisujących sześcian 
stopnie = add_cube( NULL, 0, 40, 80, 10, £f, OFF 
add_object_to axes( scena, ULok, stopnie ); 


) 


// przesunięcie każdego stopnia o 120 jednostek 
// w kierunku osi Z 
// obrót wokół osi Y o i*KAT OBROTU 
// przesunięcie o 20 jednostek w górę 
ULok = scena->l_axe3; 
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1=0; 


while( ULok != NULL ) 

( 
if( ULok->object != NULL ) 
( 


move( ULok->object->l vertices, GLOBAL,0,0,100); 
rot_y ( ULok->object->l vertices, GLOBAL, 
i*KAT_OBROTU ); 
move( ULok->object->l vertices, GLOBAL, 0, 10*i,0); 
1++; 
) 
ULok = ULok->next; 


// obliczamy współrzędne do rysowania i wyświetlamy 

// kolejno wszystkie stopnie 
observer view _scene( scena, 640, 480, 1, 1 ); 
DrawScene( scena ); 


getch (); 
clearviewport(); 
// przesunięcie obserwatora 
move (kamera->position->l vertices,GLOBAL,100,170,250); 
LookAtTarget( kamera, NULL, 0, 0, 0 ); 
observer _view scene ( scena, 640, 480, 1, 1 ); 
DrawScene( scena ); 


getch(); 

// obrót schodów 
ULok = scena->l axe3; 
while( ULok != NULL ) 


i 
if( ULok->object != NULL ) 
rot_z( ULok->object->l vertices, GLOBAL, 45 ); 
ULok = ULok=>next; 

) 


clearviewport(); 
observer view scene( scena, 640, 480, 1, 1 ); 
DrawScene( scena ); 


// zamknięcie grafiki i skasowanie obiektów 
getch(); 
closegraph (); 
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Rys. 4.9a 


delete_scene ( 


4.9b 


Rys. 
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Rys. 4.9c 
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ROZDZIAŁ 5 


Algorytmy wyznaczania 
widocznych części 
powierzchni 


Metody, które opisaliśmy w poprzednich rozdziałach, mogą służyć jedynie do 
tworzenia obrazów zbiorów punktów poddanych uprzednio pewnym transfor- 
macjom i rzutowanych na płaszczyznę w układzie obserwatora. Punkty połą- 
czone w odpowiednie struktury definiują równocześnie wierzchołki, krawędzie 
1 ściany obiektów trójwymiarowych. Wizualizacja tych obiektów na ekranie 
monitora może być sprowadzona do wyświetlenia odpowiedniej liczby linii 
przedstawiających krawędzie brył znajdujących się w obrębie ostrosłupa wi- 
dzenia. Uzyskujemy w ten sposób obraz "druciany", nieczytelny, gdyż jest na 


(AT ( 
ANŻ 


ć 


A 


Realistyczna grafika 3D 111 


nim zbyt wiele figur (rys. 5.1). Aby obraz był czytelny, należy usunąć zasłonięte 
krawędzie. Można także pokazywać tylko widoczne części ścian. 


Znajdowanie i usuwanie z obrazu niewidocznych części krawędzi ma szerokie 
zastosowanie w przypadku, gdy urządzeniem wyjściowym, na którym tworzy- 
my obraz, jest ploter lub drukarka. W przypadku wizualizacji na monitorze 
ekranowym, co jest zasadniczym tematem naszej książki, lepsze efekty osią- 
gamy stosując algorytmy znajdujące widoczne części powierzchni. W tym 
rozdziale przedstawimy budowę takiego algorytmu. Jego idea została zaczerp- 
nięta z [1]. Procedura generująca pocieniowane obrazy z programu demonstra- 
cyjnego realizuje właśnie ten algorytm. Algorytm ten bywa czasami nazywany 
scanline (co oznacza przeglądanie obrazu liniami poziomymi). Zanim opiszemy 
jego strukturę, przedstawimy ogólne cechy metod tego typu. 


5.1. Wprowadzenie 


Zaawansowane metody wizualizacji możemy podzielić na dwie zasadnicze 
klasy: 


m metody przestrzeni danych, 


M metody przestrzeni obrazu. 


Czynnikiem decydującym o przynależności danego algorytmu do jednej z tych 
klas jest przestrzeń, w której operuje. 


Algorytm przestrzeni danych wykonuje obliczenia i porównania na elementach 
sceny w przestrzeni trójwymiarowej. Taki sposób daje najczęściej dokładny 
analityczny wynik (na przykład efektem jego pracy może być tablica współrzęd- 
nych końców wszystkich widocznych odcinków) i dlatego jest kosztowniejszy 
niż algorytm przestrzeni obrazu. Przez koszt algorytmu rozumiemy oczywiście 
ilość czasu i pamięci komputera niezbędnej do wygenerowania obrazu. Im 
więcej pamięci i im dłużej program liczy, tym jest kosztowniejszy. 


Metody przestrzeni obrazu badają dwuwymiarowe rzuty obiektów na płaszczy- 
znę rysunku i dopiero gdy rzuty mają części wspólne, jest rozstrzygane zasła- 
nianie w przestrzeni trójwymiarowej (badanie, który z obiektów leży bliżej 
obserwatora). Te metody stosuje się głównie na monitorach rastrowych. Do 
nich zaliczamy algorytm, według którego działa program demonstracyjny. 


Koszt obliczeń w pierwszej grupie metod zależy głównie od liczby obiektów i ich 
złożoności. W przypadku algorytmów przestrzeni obrazu koszt jest funkcją 
rozdzielczości generowanego obrazu. Im więcej punktów naekranie, tym dłużej 
pracuje program. 


_.h_h_55Qd.ŃQQ«QQQÓóÓóÓÓÓŻ[Ż- r — 


112 


M. Domaradzki, R. Gembara 


5.2. Przykładowe algorytmy 


W jaki sposób wyświetlać obraz tak, aby niewidoczne fragmenty ścian zostały 
usunięte? 


Najprostsze rozwiązanie polega na posortowaniu wszystkich ścian względem 
ich odległości od obserwatora. Tak posortowane ściany wyświetlamy w kolej- 
ności od najdalszej do najbliższej. Każdą ścianę rysujemy w odpowiednim 
kolorze. Obraz powstaje w taki sam sposób, w jaki malarz nakłada farbę na 
płótno. Jednak tak prosty algorytm działa prawidłowo tylko dla nie przecina- 
jących się wypukłych ścian. W pozostałych przypadkach zawodzi. Jak roz- 
strzygnąć, która ze ścian na rysunku 5.2 leży bliżej nas? 


Rys. 5.2 


Innym sposobem jest zastosowanie reguły "dziel i zwyciężaj”. Dzielimy reku- 
rencyjnie prostokąt zawierający rzut sceny (ekran komputera) na prostokątne 
fragmenty (rys. 5.3). Jeżeli w danym fragmencie znajdują się rzuty kilku ścian, 
to dzielimy go na cztery równe prostokąty itd. W ten sposób obraz zostanie 
podzielony na wiele prostokątnych części. Najmniejsze z nich mają rozmiary 
jednego piksla. Dalsze podziały nie mają sensu, ponieważ piksel jest najmniej- 
szym elementem obrazu, jaki można wyświetlić. W każdej z tych części jest 
niewiele rzutów ścian. Rozstrzyganie problemu zasłaniania jest prowadzone 
osobno dla każdego prostokąta — jest to efektywniejsze niż rozstrzyganie 
zasłaniania dla całego ekranu. 
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Kolejna metoda to algorytm z buforem głębokości. Rozpatrujemy rzuty ścian 
na ekran, traktowany przez nas jako tablica piksli. Dla każdego piksla, 
należącego do rzutu każdej ściany znajdującej się na scenie, obliczamy odle- 
głość trójwymiarowego przeciwobrazu tego piksla od obserwatora. Obliczoną 
wartość zapisujemy do tablicy-bufora głębokości w wierszu i kolumnie odpo- 
wiadających położeniu tego piksla na ekranie. Zapis jest wykonywany wyłącz- 
nie, gdy zapisywana wartość jest mniejsza od zawartości tablicy. Jeżeli jest 
spełniony ten warunek, to do tablicy kolorów zapisujemy kolor ściany, do której 
należy dany punkt. Tablica głębokości powinna być na początku zainicjowana 
maksymalnie dużymi liczbami, a tablica kolorów — kolorem tła. Po wykonaniu 
tych operacji dla wszystkich rzutów ścian tablica kolorów może być, punkt po 
punkcie, wyświetlona na ekranie komputera — w niej mieści się obraz sceny. 


Jest to algorytm stosunkowo łatwy do implementacji. Szybkość jego działania 
zależy w głównej mierze od liczby piksli na ekranie. Dzieje się tak dlatego, że 
przy wzroście liczby ścian na scenie maleją najczęściej ich rozmiary, a co za 
tym idzie, zmniejsza się także liczba piksli zawartych w rzutach tych ścian. Te 
zalety sprawiły, że algorytmy z buforem głębokości można realizować nie tylko 
jako program komputerowy, ale także w postaci układów logicznych, czyli 
sprzętowo. Daje to niezwykłe przyśpieszenie ich pracy. W stacjach graficznych 
Silicon Graphics są montowane specjalizowane karty graficzne zapewniające 
sprzętową realizację takiego algorytmu. W połączeniu z szybkimi procesorami 
umożliwia to tworzenie interakcyjnej grafiki komputerowej w czasie rzeczywi- 
stym, czyli tzw. Virtual Reality. Użytkownik takiego systemu ma możliwości 
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tworzenia scen z pocieniowanymi, kolorowymi przedmiotami, a po założeniu 
specjalnych "okularów" i rękawic istnieje możliwość wejścia w istniejący wyłą- 
cznie w pamięci komputera wirtualny, trójwymiarowy świat. 


Powrócimy jednak do rzeczywistości. Dostępne dla przeciętnego użytkownika 
komputery jeszcze przez kilka najbliższych lat nie będą miały takich możliwo- 
ści. Pewną wadą algorytmu z buforem głębokości jest jego "pamięciożerność”. 
Bufor głębokości jest tablicą liczb rzeczywistych 8-bajtowych (na tylu bajtach 
jest zapisywana liczba rzeczywista w Borland C++). Natomiast tablica kolorów 
powinna dla każdego piksla zawierać trzy składowe koloru: R,G,B, czyli trzy 
bajty. Daje to w sumie 11 bajtów na piksel. Zatem dla obrazu 320 na 200 piksli 
potrzeba 320*200*11 = 687,5 KB pamięci. Dla obrazu 1024 na 768 (maksymal- 
na rozdzielczość ekranu typowej karty SVGA ) należy zarezerwować prawie 
8,5 MB, nie licząc pamięci potrzebnej dla struktur opisujących scenę. 


5.3. Algorytm przeglądania liniami poziomymi (algorytm scanline) 


Ten algorytm nie ma tak wygórowanych wymagań pamięciowych jak algorytm 
z buforem głębokości, a jego szybkość działania jest porównywalna, ponieważ 
wykorzystuje pewne zależności zachodzące pomiędzy rzutami ścian, a samymi 
ścianami. Ze względu na duże rozmiary (prawie 500 linii) nie zamieszczamy 
kodu źródłowego funkcji scanline, realizującej ten algorytm. Na dyskietce w 
pliku scan.c znajduje sama funkcja scanline, a także zbiór pomocniczych 
procedur. 


Podobnie jak poprzednio, ekran, na którym będzie rysowany obraz, jest ekra- 
nem rastrowym — składa się z piksli. Możemy go porównać do prostokątnej 
tablicy. Każdym elementem takiej tablicy jest piksel (patrz rys. 5.4). Szerokość 
obrazu wynosi Xmax piksli, a wysokość Ymax. 


Ogólna wersja tego algorytmu nie wymaga, aby ściany (a co za tym idzie także 
ich rzuty) były wielokątami wypukłymi. W tej wersji ściany są trójkątami 
i mogą się wzajemnie przenikać. 


W jaki sposób działa algorytm przeglądania liniami poziomymi? Na obrazie 
wyjściowym są rzuty wszystkich krawędzi znajdujących się wewnątrz ostrosłu- 
pa widzenia. Ponieważ większość operacji będzie wykonywana na płaszczyźnie 
rzutowania - ekranie, rzuty krawędzi będziemy odtąd nazywać krawędziami. 
Przyjmujemy, że końce krawędzie mają następujące współrzędne: począ- 
tek(xd,yd), koniec(xg,yg), yd<yg. 


Na samym początku odrzucamy krawędzie poziome, a pozostałe sortujemy w 
następujący sposób: 
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Obraz wyjściowy jest przeglądany linia po linii. Dla każdej poziomej linii (linie 
od 0 do Ymax-1) tworzymy listę krawędzi LK, dla których współrzędna yd 
odpowiada numerowi tej linii. Ta lista musi być posortowana zgodnie z rosną- 
cymi wartościami xd. Gdy spotkamy krawędzie zaczynające się w tym samym 
punkcie, sortowanie odbywa się względem rosnących odwrotności współczyn- 
nika c. 


c= Ed 


Każdy element listy krawędzi składa się z pól xd, yg, ci dwóch wskaźników do 
ścian, które zawierają daną krawędź. W strukturach FACE opisujących ściany 
znajduje się pole inface, przyjmujące dwie wartości: TRUE — gdy podczas 
przeglądania algorytm działa wewnątrz rzutu tej ściany, FALSE — gdy opusz- 
cza rzut ściany. Dodatkowo ściany wewnątrz, których rzutów znajdujemy się 
w trakcie przeglądania, są łączone w listę ścian aktywnych LFA. Dla każdej 
ściany należącej do LFA pole inface ma wartość TRUE. 


Struktura LLINE służy do przechowywania informacji o rzutach krawędzi. 


typedef struct l // opisuje rzut krawędzi na ekran 
( // z takich elementow jest złożona lista krawędzi 
double xd; 
int Y9; 
double c; 


// odwrotność współczynnika kierunkowego krawędzi 
FACE huge  *facel; 

// ściany, do których należy ta krawędź 
FACE huge  *face2; 
VERTEX huge *vertl;  // końce linii 
VERTEX huge *vert2; 
struct 1 huge *next; // następna krawędź na liście 
struct 1 huge *nextAL; // następna krawędź aktywna 
struct (1 huge *prevAL; // poprzednia krawędź aktywna 

JLLINE; 


Pole *next służy do łączenia struktur LLINE w listę krawędzi. Algorytm 
generuje tyle list, ile poziomych linii zawiera obraz wyjściowy. Listy krawędzi 
w liniach obrazu, w których nie zaczyna się żadna krawędź, są puste. Listy 
krawędzi są dołączane do tablicy krawędzi. Na rysunku 5.5b przedstawiono 
strukturę tablicy krawędzi (Row Table) zbudowanej z list elementów typu 
LLINE, wygenerowaną dla rzutu ściany z rys. 5.5a. 
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Rys. 5.5 


Na początku bufor ekranu wypełniamy kolorem tła. Rzut rysowanej sceny 
przeglądamy liniami poziomymi w kierunku ich rosnących numerów (patrz rys. 
5.2). Każdą linię przeglądamy z lewa na prawo. W każdej linii tworzymy listę 
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aktywnych krawędzi LKA i operujemy tylko na jej elementach. LKA zawiera 
wyłącznie te krawędzie, które przecinają daną linię. 


W pliku scan.c, nazwy opisywanych przez nas list są nieco inne. Zmienna 
globalna row_table jest tablicą krawędzi. Zmienna ActiveLinesList wskazuje 
na początek listy aktywnych krawędzi LKA. Ponieważ LKA jest listą dwukie- 
runkową, zdefiniowano dla niej korzeń (RootAL) i wartownika (GuardAL). 


ALGORYTM SCANLINE 


— znajdź współrzędne wszystkich obiektów w układzie obserwatora, 

— rzutuj wyłącznie punkty znajdujące się w ostrosłupie widzenia, 

— utwórz LK[Ymax], 

— wyzeruj LKA, 

— linia = numer pierwszej linii, w której zaczyna się jakakolwiek krawędź, 

— dopóki linia < Ymax, powtarzaj: 

(a) dołącz do LKA listę krawędzi wskazywaną przez LK[linia], zachowując 
uporządkowanie jej elementów w kolejności rosnących wartości xd, wyczyść 
LFA; 

(b) przeglądaj LKA, zmieniając dla każdej krawędzi pole inface na przeciwne 
w ścianach, do których należy ta krawędź. Jeżeli wchodzimy do rzutu ściany 
(inface=TRUE), to dołączamy ścianę do LFA. Jeżeli wychodzimy z LFA (infa- 
ce=FALSE), do usuwamy ścianę z LFA. Jeżeli znajdziesz dwa sąsiednie ele- 


menty o polach xd, których różnica jest większa od 1, to możesz przejść do 
następnego punktu. W przeciwnym wypadku powtarzaj (b); 


(c) teraz dla każdego z punktów leżących pomiędzy przecięciem wiersza ekranu 
o numerze linia z krawędzią opisywaną przez poprzedni element LKA a 
przecięciem wiersza linia z krawędzią opisywaną przez aktualny element LKA 
znajdź najbliższą ścianę spośród ścian znajdujących się na LFA; 


(d) jeżeli jest jeszcze jakiś element LKA, to wróć do punktu (b); 
(e) usuń z LKA elementy, dla których yg = linia; 


(f) we wszystkich pozostałych elementach LKA zmodyfikuj xd = xd + ©, 
otrzymasz w ten sposób punkty przecięcia krawędzi z wierszem ekranu 0 
numerze linia+1; 


(g) posortuj LKA względem nowych wartości xd; 


(h) linia = linia + 1. 
Koniec. 
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Rys. 5.6 


Aby idea działania tego algorytmu stała się nieco czytelniejsza, spróbujmy 
przeanalizować jego działanie na przykładzie obrazu zbudowanego z rzutów 
dwóch ścian. Rysunek 5.6 ilustruje tę sytuację. 


(1) Linia p nie przecina żadnej krawędzi. LK[p], LKA, LFA są puste. 


(ii) W czasie przeglądania linii r LK i LKA zawierają po dwie krawędzie, 
LK=LKA. Obie krawędzie należą do ściany 2. W punkcie A zmieniamy zmienną 
inface na TRUE w strukturze opisującej ścianę 2 i dołączamy tę ścianę do listy 
aktywnych ścian LFA. Dla wszystkich punktów na odcinku AB jedyną widocz- 
ną ścianą jest ściana 2. Następnym elementem LKA jest druga krawędź 'ściany 
2 przecinająca się z linią r w punkcie B. Współrzędna x tego punktu jest równa 
wartości pola xd w strukturze LLINE opisującej tę krawędź. Po dotarciu do B 
wartość zmiennej inface zmieniamy na przeciwną, tzn. FALSE, a ścianę 2 
usuwamy z LFA. Ponieważ LKA nie zawiera dla tej linii więcej elementów, 
możemy zakończyć przeglądanie. 


Giii) W linii s postępujemy jak poprzednio wyłącznie na odcinku CD. Na LFA 
znajduje się wtedy tylko ściana 1. Po dojściu do D należy dołączyć do listy ścian 
aktywnych ścianę 2. Teraz, poruszającsię od D do E, należy dla każdego punktu 
rozstrzygać, która ze ścian leży bliżej obserwatora. Porównujemy odległości od 
obu ścian. To zagadnienie jest opisane poniżej. 


UWAGA. Jeżeli na scenie nie ma przenikających się ścian, to dla wyznaczania 
nie zasłanianej ściany moglibyśmy zastosować następującą zasadę: nowe obli- 
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czenia dla zbadania zasłaniania wykonujemy tylko wtdey, gdy linia przeglą- 
dająca opuszcza ścianę zasłaniającą inną. Taki sposób postępowania może 
bardzo przyśpieszyć generowanie obrazów. 


(iv) Linia t. Na odcinku GH widoczna jest ściana 1, a na odcinku IJ ściana 2. 


W punkcie (d) algorytmu scanline mamy dane współrzędne ekranowe punktu 
i listę ścian aktywnych. Oznacza to, że prosta poprowadzona od obserwatora 
poprzez ten punkt, znajdujący się na płaszczyźnie rzutowania, przecina wszy- 
stkie ściany z LFA. Na ekranie należy wyświetlić w tym punkcie kolor tej 
ściany, która jest widoczna, czyli leży najbliżej obserwatora. Aby wyznaczyć 
odległość od ściany, należy znaleźć przeciwobraz P punktu P' z płaszczyzny 
rysunku, leżący na danej ścianie. Jest to zagadnienie odwrotne do omawianego 
wcześniej rzutowania perspektywicznego. W tamtym przekształceniu punkt P 
był nazywany obrazem punktu P. Przypomnijmy wzory pozwalające na obli- 
czenie współrzędnych punktu po rzutowaniu: 


EP o — 

x = lul 
EP oOv 

EP o —— 

s lwl 
Je EP o0Ov 


Znamy także równanie ściany (n — wektor normalny do ściany): 
nop=c,  n=[x,,y,»Z,] 
równanie prostej łączącej obserwatora z punktami P' iP: 
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Uwzględniając równania ściany i prostej, otrzymujemy zależności pozwalające 
na obliczenie współrzędnych przeciwobrazu: 
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Aby skrócić czas obliczeń, ostatnie dwa wzory możemy uprościć. Przed urucho- 
mieniem algorytmu przeglądania liniami poziomymi wszystkie obiekty znaj- 
dujące się na scenie sprowadzamy do układu obserwatora (rys. 5.7). 


Rys. 5.7 


Wtedy obserwator znajdzie się w punkcie [0,0,-1], a także: 


=" — =[1,0,0] 
lul 
—"—=[0,1,0] 
lwl 
v=[0,0,1] 
Wzory na współrzędne przeciwobrazu uproszczą się do postaci: 
. c+z_ 
x x,+ty y,tz, 
x=x't 
y=y t 
z=t-l 


Mając wyznaczoną najbliższą ścianę i współrzędne przeciwobrazu punktu 
położonego na ekranie, możemy znaleźć natężenie oświetlenia ściany w tym 
punkcie. To zagadnienie będzie tematem następnego podrozdziału. 
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5.4. Modelowanie oświetlenia 


W wyniku działania algorytmu scanline otrzymujemy wszystkie niezbędne 
informacje potrzebne do wizualizacji obiektów. Dla każdego punktu obrazu 
znamy wskaźnik do najbliższej ściany. Możemy pobrać kolor ze struktury 
FEATURE dołączonej do struktury FACE opisującej tę ścianę i wyświetlić na 
ekranie piksel o takim kolorze. Jak się łatwo domyślić, otrzymany obraz będzie 
przedstawiał pokolorowane kontury obiektów. Na pewno nie będzie to wyglą- 
dać naturalnie. 


Aby realistycznie przedstawić obiekty trójwymiarowe, należy dodatkowo uwz- 
ględniać takie elementy obrazu, jak kolor przedmiotów, połyskliwość ścian, 
położenie i barwa źródła światła itp. W celu modelowania oświetlenia można 
zastosować model zaproponowany przez Bui-Tuong Phong'a. Model Phonga 
zakłada istnienie kilku współczynników charakteryzujących zdolność powierz- 
chni do odbijania światła. W zależności od tych współczynników ściana może 
być na przykład połyskliwa bądź matowa. 


Na wstępie należy zaznaczyć, że źródłem oświetlenia w tym modelu mogą być 
zarówno lampy, jak i tło. Do opisu barwy światła stosujemy model RGB. Skrót 
RGB oznacza trzy składowe światła: czerwoną, zieloną i niebieską. Przyjmijmy, 
że każdy z tych składników może się zmieniać w zakresie od 0 do 1. Dobierając 
odpowiednio wartości R,G,B można otrzymać dowolny kolor możliwy do wy- 
świetlenia. Oto przykładowe wartości składowych dla różnych barw (w kolej- 
ności [R,G,B]): 


CZARNY - [0, 0, 0) 
BIAŁY — [1, 1, 1) 
CZERWONY - [1, 0, 0] 
ZIELONY — [0, 1, 0) 
NIEBIESKI  - [0, 0, 1) 
KARMAZYNOWY — [1, 0, 1) 
ŻÓŁTY — [1, 1, 0] 
CYAN — [0, 1, 1) 


Każdy przedmiot charakteryzuje się pewną zdolnością do odbijania oświetlenia 
tła. Współczynniki ktR, ktG, ktB, należące do przedziału [0,1], charakteryzują 
tę zdolność. Oświetlenie tła powoduje, że części obiektów, do których nie 
dociera bezpośrednio światło lamp, mogą być widoczne (ich kolor nie jest 
czarny). 


Dla powierzchni matowych decydujące znaczenie o ich jasności i kolorze mają 
współczynniki odbicia rozproszonego krR, KrG, krB. Im mają większe wartości, 
tym jaśniejsza jest powierzchnia. 
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Przedmioty gładkie odbijają światło w jednym kierunku. Przy oglądaniu ich z 
innych kierunków są ciemne. Decyduje o tym odbicie zwierciadlane (kierunko- 
we). Współczynnik kz opisuje tę właściwość powierzchni. Zauważmy, że do 
opisu odbicia rozproszonego i oświetlenia tła użyto trójek współczynników: 
jeden współczynnik dla jednej składowej światła. W przypadku odbicia kierun- 
kowego o barwie światła odbitego nie decyduje kolor powierzchni, ale barwa 
światła padającego na tę powierzchnię. Dlatego w modelu używamy tylko 
jednego współczynnika kz. 


Rys. 5.8 


Na rysunku 5.8 przedstawiono model odbicia promienia światła od powierzch- 
ni. Wszystkie oznaczone wektory są znormalizowane, to znaczy, że ich długości 
są równe 1. Rozszyfrujmy oznaczenia: 


n — wektor normalny (prostopadły) do powierz- 
chni w punkcie odbicia, 

1 — _ wektor kierunku, z którego pada światło, 

z — promień światła odbitego od powierzchni, 

e — kierunek, z którego patrzy obserwator, 

0 — kąt rozwarcia stożka. 
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Wzór opisujący ten model oświetlenia przytoczymy w wersji dla światła mono- 
chromatycznego. W funkcji scanline jasność danego punktu na ścianie jest 
obliczana właśnie w ten sposób. W programie demonstracyjnym "RAY3D.EXE"' 
znajdującym się na dyskietce funkcjonują trzy podobne formuły, każda dla 
innej składowej RGB. 


— z 
I=k,I, + [kr (091) +k,(z0e)=|qzjgh 
kt It — ten czynnik określa składową oświetlenia tła, 
nol — zależy od wzajemnego położenia źródła światła i 


powierzchni, wpływając na odbicie rozproszone. Im większy kąt padania świat- 
ła, tym jaśniejsza będzie w tym punkcie powierzchnia. Kąt podania jest liczony 
względem wektora normalnego do powierzchni. Jeżeli kąt padania wynosi 0 
stopni, to znaczy, że światło pada na powierzchnię prostopadle i oświetlenie 
jest maksymalne. Jeżeli kąt padania jest bliski 90 stopni, to znaczy, że 
promienie światła podążają niemal równolegle do powierzchni, w efekcie czego 
jasność powierzchni jest minimalna. 


kz (z0e)””* — charakteryzuje odbicie zwierciadlane. W tym przypadku ważny 
jest kąt między wektorami z, e. Jeśli promień odbity od powierzchni pada 
wprost na obserwatora lub kiedy kąt pomiędzy promieniem odbitym a kierun- 
kiem patrzenia jest niewielki, to obserwator zauważy jasną plamę na powierz- 
chni. W taki sposób odbijają światło na przykład przedmioty metalowe. Przy 
odpowiednim ustawieniu się względem źródeł światła zauważamy na nich 
jasne 'blinki". Parametr m należy do przedziału [0..100]. Ten parametr służy 
do regulowania rozmiarów i rozmycia odblasków na powierzchni. Dla powierz- 
chni matowych przyjmuje małe wartości (duży kąt 6, duża jasna plama odbla- 
sku), a dla błyszczących duże (mały kąt 6). Dla idealnego zwierciadła m 
przyjmuje wartość nieskończoną. 


IZ — natężenie światła emitowanego przez 
źródło, 
d — odległość oświetlanego punktu powierzchni 


od źródła światła, 


do — pewna stała dobierana doświadczalnie w 
celu uzyskania lepszych obrazów. 


Trzy ostatnie współczynniki nie są przez nas używane w funkcji scanline. 
Zwiększają one ilość obliczeń, nie wnosząc dodatkowej jakości. Ich zadaniem 
ma być modelowanie spadku natężenia oświetlenia wraz ze wzrostem odległo- 
ści od źródła światła. Gdy pominiemy te współczynniki, użytkownik nie musi 


wy wn 
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się troszczyć o to, czy prawidłowo dobrał jasność lampy, czy lampa nie znajduje 
się za blisko lub za daleko od oświetlanego przedmiotu, w wyniku czego obraz 
mógłby zostać prześwietlony (zupełnie jak w zwykłym aparacie fotograficz- 
nym). Obiekty można ustawiać dowolnie daleko od lampy — nigdy obraz nie 
będzie "prześwietlony" lub "niedoświetlony". Oczywiście Czytelnik może samo- 
dzielnie rozszerzyć funkcję scanline, tak aby uwzględniała i ten czynnik — 
wystarczy zmodyfikować wyrażenia obliczające oświetlenie powierzchni w 
danym punkcie. 


Używana przez nas już wcześniej struktura FEATURE służy do przechowywa- 
nia informacji o współczynnikach charakteryzujących powierzchnię: 


typedef struct feat 

( 
int number; // potrzebne przy zapisywaniu do pliku 
double  ktR; 

// cztery współczynniki charakteryzujące powierzchnię 
double ktG; 

// współczynniki odbicia oświetlenia tła dla kolorów RGB 
double  ktB; 
double  krR; 

// kolor powierzchni jako 3 wsp. odbicia 

// światła rozproszonego 
double  krG; 
double  krB; 


double kz; 

// współczynniki odbicia zwierciadlanego 
int m; 
struct feat *next; 

) FEATURE; 


Struktura LAMP opisuje parametry źródła światła. Na scenie może znajdować 
się wiele źródeł światła, dlatego istnieje pole next służące do łączenia struktur 
w listę. Każda lampa jest traktowana jako punktowe źródło, wysyłające pro- 
mienie równomiernie we wszystkich kierunkach. 


typedef struct lp 
t 
char *name; 
VERTEX huge *position; // położenie lampy 
int number; 
double IZ; 


Realistyczna grafika 3D 125 


// natężenie światła 
double IR; 
// trzy składniki barwy światła: od 0 do 1 
double IG; 
double IB; 
struct lp *next; 
// wskaźnik do kolejnej lampy 
) LAMP; 


Na koniec dodajmy, że współczynnik kz może zostać zastąpiony funkcją zależną 
od kąta padania światła. Natężenie odbicia zwierciadlanego światła dla nie- 
których materiałów maleje wraz ze wzrostem tego kąta. Na przykład, patrząc 
pod dużym kątem na szybę w oknie zauważamy, że szyba nie tylko przepuszcza 
światło z zewnątrz, ale także odbija obrazy przedmiotów znajdujących się w 
pokoju. Gdy staniemy naprzeciw okna, nie zauważymy tego efektu. 


5.5. Wygładzanie powierzchni 


Bryły budujemy z trójkątnych ścian, dlatego obrazy uzyskiwane z funkcji 
scanline, rozszerzonej o wyznaczanie natężenie oświetlenia w danym punkcie, 
składają się z "kanciastych" obiektów: sześcianów, prostopadłościanów, ostro- 
słupów. Nie można uzyskać obrazów przedstawiających "gładkie" bryły: stożki, 
kule, walce. Nie pomaga w tym przypadku nawet tworzenie brył o bardzo dużej 
liczbie małych ścian. Na rysunku 5.9 widzimy "kulę" przybliżaną coraz większą 
liczbą trójkątów. 


Aby uzyskać gładki wizerunek powierzchni na podstawie jej przybliżonej 
reprezentacji złożonej z trójkątów, należy zastosować jedną z metod wygładza- 
nia powierzchni. 


Najprostszą i najbardziej ekonomiczną metodą jest obliczanie oświetlenia w 
wierzchołkach ściany (wartość średnia z jasności ścian, do których należy dany 
wierzchołek) i aproksymowaniu za pomocą uzyskanych wartości oświetlenia 
punktów wewnątrz tej ściany. Jest to metoda Gourauda. Obrazy uzyskane w 
ten sposób są znacznie bardziej realistyczne niż przy stałej intensywności 
oświetlenia ścian. Zasadniczą wadą metody Gourauda jest to, że nie można 
uzyskać naturalnie wyglądających odblasków światła na powierzchniach ma- 
jących symulować metaliczne. 


Znacznie kosztowniejszą, ale dającą zdecydowanie lepsze efekty, jest metoda 
Phong'a. W tym przypadku interpoluje się nie natężenie oświetlenia w wierz- 
chołku, lecz wektory normalne do powierzchni. 
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Rys. 5.9a 


Rys. 5.9b 
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Rys. 5.9c 


W pierwszym kroku są wyznaczane wektory normalne do wszystkich ścian 
znajdujących się w ostrosłupie widzenia (rys. 5.10). 


Potem dla każdego wierzchołka wyznaczamy wersor normalny, będący śred- 
nią arytmetyczną wektorów normalnych do ścian zawierających ten wierzcho- 
łek (rys. 5.11). 


Teraz dla każdej ściany mamy wyznaczone trzy uśrednione wektory nl,n2,n3. 
W punkcie o współrzędnych ekranowych x,y kierunek wektora n normalnego 
do powierzchni ściany wyznaczamy interpolując wektory nl,n2,n3 (rys. 5.12). 


Wektory m i np otrzymujemy z zależności: 


p Jsan J2 J 1 7) scan 
Ry ye 2 Y1-Y 
1 2 1 2 

Jscan  J3 Ją J scan 

n = lu 0, Ra 
JI I 3 J1 J 3 


Powyższe operacje są przeprowadzane na współrzędnych rzutów wierzchoł- 
ków. W obu wzorach wartości y1, y2, y3 są współrzędnymi punktów zaczepienia 
odpowiednich wektorów, a wartość Yscan jest numerem poziomej linii przeglą- 
dania. Znając wektory ni, np w podobny sposób wyznaczamy wektor n. 
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Rys. 5.10 


M=A*R*R*R, 


Rys. 5.11 
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Po tym stosujemy model oświetlenia opisany w poprzednim podrozdziale dla 
obliczenia jasności ściany w punkcie (x,y). Na rysunku 5.13 przedstawiamy 
obraz kuli uzyskany za pomocą funkcji scanline z wygładzaniem powierzchni 
metodą Phong'a. Tekst programu zamieszczamy poniżej. 


// prog7.c 


łinclude "ray.h" 
include "listy.c" 
kinclude "mat.c" 
include "obiekty.c" 
łinclude "rysuj.c" 
+kinclude "scan.c" 
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void main( void ) 


( 


/7/ 


/7/ 
/7/ 
// 
// 
/7/ 
// 


OBJECT *bryla; 

AXES  *UkladLokalny; 

SCENE  *scena; 

OBSERVER *kamera; 

LAMP  *lampa; 

FEATURE  *lsniacy; 

int gdriver, gmode=0 , errorcode; 


inicjalizacja grafiki 

gdriver = installuserdriver("svga256", NULL ); 
Na dyskietce znajduje się plik svga256 
umożliwiający otwieranie grafiki w 

256 kolorach w pakietach Turbo C i Borland C. 
Trzeci parametr funkcji 

initgraph jest ścieżą do katalogu, w którym 
znajduje się ten plik. 

initgraph(śsgdriver, śgmode, "c:Nlborlandciibgi"); 


// parametry powierzchni 


lsniacy = new feature( 0, 0, 0, .3, .3, .3, 1, 15 ); 


scena = add scene( "", .5, .5, .5 ); 
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UkladLokalny = add_axes( "", 0, 40, 40, 40 ); 

add axes_to scene( scena , UkladLokalny ); 

bryla = add_sphere( NULL,0,20,8,8,0FF,lsniacy,ON ); 
add_object_to axes( scena, /UkladLokalny, bryla ); 
kamera = add observer( 100, 100, -200, .4, .4 ); 
add_observer_to _scene( scena, kamera ); 

lampa = add _lamp( "", 0, 200, 100, -200, 1, 1, 1, 1); 
add_lamp_to scene( scena, lampa ); 


PaletaKolorow(); 
observer _view scene( scena, 320, 200, 1.2, 1 ); 
ScanLine( scena, 320, 200, 1.2, 1 ); 


// zamknięcie grafiki 
getch(); 
closegraph(); 
delete _scene( scena ); 


) 


5.6. Przykład 


Program, który przedstawimy w tym rozdziale, wykorzystuje większość opisa- 
nych przez nas funkcji. Mamy nadzieję, że po jego przeanalizowaniu Czytelnik 
będzie będzie mógł samodzielnie go rozbudować lub wykorzystać odpowiednie 
funkcje we własnych programach. 


Program umożliwia poruszanie się obserwatora po wirtualnej przestrzeni 
trójwymiarowej, w której umieszczono źródło światła i kilka brył. Obserwator 
może być przemieszczany równolegle do osi układu współrzędnych za pomocą 
następujących klawiszy: 


1, l — góra, dół, 
<, > — lewo, prawo, 


PageUp, PageDown — do przodu (w kierunku osi z), do tyłu (przeciwnie do 
osi z). Obserwator zawsze patrzy na wybrany obiekt. Naciskając TAB możemy 
zmieniać kierunek patrzenia na inny obiekt. Klawisze + i — służą do zmiany 
pola widzenia obserwatora (powiększenie/pomniejszenie oglądanego obrazu). 
Po naciśnięciu SPACE następuje wywołanie funkcji scanline, która generuje 
obraz z usuniętymi niewidocznymi fragmentami powierzchni i cieniowaniem. 
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Struktura tego programu jest bardzo prosta. Na początku tworzymy zbiór 
układów lokalnych i związanych z nimi brył. Później dołączamy do sceny 
obserwatora i źródło światła, ustalamy paletę kolorów, wykonujemy rzutowa- 
nie i wyświetlamy szkieletowy rysunek wszystkich obiektów. Potem wchodzi- 
my do pętli wyboru, w której w zależności od naciśniętego klawisza są wywo- 
ływane odpowiednie funkcje. Na przykład jeżeli obserwator ma być przesunię- 
ty, to czyścimy ekran, przesuwamy obserwatora (funkcja move), kierujemy go 
tak, aby patrzył na początek układu lokalnego wskazywanego przez cel (fun- 
kcja LookAtTarget), wykonujemy rzutowanie i rysujemy scenę. Na końcu za 
pomocą procedury delete_scene zwalniamy pamięć zajmowaną przez struktu- 
ry. 


// prog8.c 

łinclude "ray.h" 
include "listy.c" 
kinclude "mat.c" 
finclude "obiekty.c" 
kinclude "rysuj.c" 
łfinclude "scan.c" 


// klawisze 

define ESC 27 

define GORA 72 

define DOL 80 

define PRAWO77 

define LEWO 75 

define PRZOD 73 // Pageup 
define TYL 81 // PageDown 
define SPACJA 32 

define TAB 9 


4: RE WE W dE = W W 


// rozmiary ekranu wyjsciowego 
+ define XEKRAN 320 
+ define YEKRAN 200 


void main( void ) 
t 
OBJECT  *bryla; 
AXES  *UkladLokalny,*cel; 
SCENE  *scenąa; 
OBSERVER *kamera; 
LAMP  *lampa; 
FEATURE  *lsniacy,*matowy; 
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double krok = 30; 

double katwidzenia=30,xx; 

int  gdriver, gmode=0 , errorcode; 
int znak; 


// inicjalizacja grafiki 
gdriver = installuserdriver("svga256", NULL ); 
// Na dyskietce znajduje się plik svga256 
// umożliwiający otwieranie grafiki w 
// 256 kolorach w pakietach Turbo C i Borland C. 
// Trzeci parametr funkcji 
// initgraph jest Ścieżą do katalogu, w którym 
// znajduje się ten plik. 
initgraph(śsgdriver, śsgmode, "c:NNborlandcNibgi"); 


// parametry powierzchni 


matowy = new feature( 0, 0, 0, 1, 1, 1, 0, 1 ); 
lsniacy = new feature( 0, 0, 0, .5, .5, .5, 1, 15 ); 
scena = add scene( "", .5, .5, .5 ); 


// utworzenie matowej kuli położonej w centrum 
// układu współrzędnych 
UkladLokalny = add axes( "", 0, 40, 40, 40 ); 
add axes_ to scene( scena , UkladLokalny ); 
bryla = add sphere( NULL, 0, 20, 8, 8, OFF, matowy, ON ); 
add object to axes( scena, UkladLokalny, bryla ); 


// utworzenie sześcianu o matowej powierzchni, 
// położonego w punkcie (50,0,70) 
UkladLokalny = add axes( "", 0, 40, 40, 40 ); 
add _axes_to _scene( scena , UkladLokalny ); 
bryla = add cube( NULL, 0, 30, 30, 20, matowy, OFF ); 
rot_y( bryla->l vertices, GLOBAL, 30 ); 
rot _x( bryla->l vertices, GLOBAL, 45 ); 
move( bryla->l vertices, GLOBAL, 50, 0, 70 ); 
move( UkladLokalny->l vertices, GLOBAL, 50, 0, 70 ); 
add object to axes( scena, UkladLokalny, bryla ); 


// stożek o wysokości 50, promieniu podstawy 20 
UkladLokalny = add axes( "", 0, 40, 40, 40 ); 
add _axes_to_scene( scena , UkladLokalny ); 
bryla = add cone( NULL,0, 50, 20, OFF, 6, lsniacy, OFF ); 
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move( bryla->l vertices, GLOBAL, —60, 0, -10 ); 
move( UkladLokalny->1 vertices, GLOBAL, -60, 0, -10 ); 
add_object_to axes( scena, UkladLokalny, bryla ); 


//walec 
UkladLokalny = add_axes( "", 0, 40, 40, 40 ); 
add_axes_to scene( scena , UkladLokalny ); 
bryla = add_ tube ( NULL, 0, 40, 40, 25, ON, ON, 
10, OFF, matowy, OFF ); 
move( bryla->l vertices, GLOBAL, -50, 0, 70 ); 
move( UkladLokalny->1l vertices, GLOBAL, -50, 0, 70 ); 
add object _to axes( scena, UkladLokalny, bryla ); 


//sześcian 
UkladLokalny = add axes( "", 0, 40, 40, 40 ); 
add _axes_to_scene( scena , UkladLokalny ); 
bryla = add cube( NULL, 0, 30, 30, 20, lsniacy, OFF ); 
move( bryla->l vertices, GLOBAL, 100, -40, 70 ); 
move( UkladLokalny->l vertices, GLOBAL, 100, —-40, 70 ); 
add_object to axes( scena, UkladLokalny, bryla ); 


// kula o połyskliwej powierzchni 
UkladLokalny = add _axes( "", 0, 40, 40, 40 ); 
add_axes_to_scene( scena , UkladLokalny ); 
bryla = add sphere( NULL,0,20,10,12, OFF, lsniacy,ON ); 
move( bryla->l vertices, GLOBAL, 50, 70, -70 ); 
move( UkladLokalny->1l vertices, GLOBAL, 50, 70, -70 ); 
add object to _axes( scena, UkladLokalny, bryla ); 


// zmienna katwidzenia może przybierać 

// wartości od 5 do 175 stopni 
xx = tan( (double) M PI*katwidzenia/360 ); 

// utworzenie obserwatora położonego 

// w punkcie (100,100,-200) 
kamera = add observer( 100, 100, -200, xx, xx ); 

// dołączenie obserwatora do sceny 
add_observer_to scene( scena, kamera ); 

// utworzenie lampy i dołączenie jej do sceny 
lampa = add_lamp( "", 0, 200, 100, -200, 1, 1, 1, 1); 
add _lamp_to_scene( scena, lampa ); 


PaletaKkolorow(); 
setcolor( 15 ); 


Realistyczna grafika 3D 135 


observer _view_scene( scena, XEKRAN, YEKRAN, 1.2, 1 ); 
DrawScene( scena ); 
// cel - wskaźnik do struktury typu AXES opisującej 


// układ lokalny w kierunku,z którego "patrzy" obserwator 
cel = NULL; 


xx=0.4; 
while( ( znak = getch() ) != ESC ) 
t 

switch( znak ) 

t 


case GORA: clearviewport(); 
move( kamera->position->l vertices, GLOBAL, 0, 
krok, 0 ); 
LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer _view scene( scena, XEKRAN, YEKRAN, 
1.2, 1 ); 
DrawScene( scena ); 
break; 
case DOL : clearviewport(); 
move (  kamera->position->l vertices, GLOBAL, 
0, -krok, 0 ); 
LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer _view scene (scena, XEKRAN, YEKRAN, 
1.2, 1 ); 
Drawscene( scena ); 
break; 
case PRAWO: clearviewport(); 
move ( kamera->position->l vertices, GLOBAL, 
krok, 0, 0 ); 
LookAtTarget( kamera, cel, 0, 0, O ); 
observer _view scene (scena, XEKRAN, YEKRAN, 
1.2, 1 ); 
DrawScene( scena ); 
break; 
case LEWO: clearviewport(); 
move ( kamera->position->l vertices, GLOBAL, 
-krok, 0, 0 ); 
LookAtTarget( kamera, cel, 0, 0, 0. ); 
observer _view scene (scena, XEKRAN, YEKRAN, 
1.2, 1); 
DrawScene( scena ); 
break; 
case PRZOD: clearviewport(); 
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move ( kamera->position->l vertices, GLOBAL, 
0, 0, krok ); 
LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer_view scene (scena, XEKRAN, YEKRAN, 
1.2, 1); 
DrawScene( scena ); 
break; 
case TYL:  clearviewport(); 
move ( kamera->position->l vertices, GLOBAL, 
0, 0, -krok ); 
LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer_view scene( scena, XEKRAN, YEKRAN, 
1.2, 1); 
DrawScene( scena ); 
break; 
// pomniejszenie kąta patrzenia obserwatora o pięć stopni 
// w efekcie uzyskamy powiększenie obrazu 
case '—' : clearviewport(); 
katwidzenia += 5; 
if( katwidzenia > 175 ) 
katwidzenia = 175; 
xx = tan( (double) M PI*katwidzenia/360 ); 
ObserverLens( kamera, XX, XX ); 
LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer view scene( scena, XEKRAN, YEKRAN, 
1.2, 1 ); 
DrawScene( scena ); 
break; 
// zwiększenie kąta patrzenia obserwatora o pięć stopni 
// w efekcie uzyskamy pomniejszenie obrazu 
case '+' : clearviewport(); 
katwidzenia -—= 5; 
if( katwidzenia < 5 ) 
katwidzenia = 5; 
xx = tan( (double) M PI*katwidzenia/360 ); 
ObserverLens( kamera, XxX, XX ); 
LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer view scene( scena, XEKRAN, YEKRAN, 
- - 1.2, 1); 
DrawScene( scena ); 
break; 
// wywołanie funkcji scanline 
case SPACJA: clearviewport(); 
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LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer _view scene( scena, XEKRAN, YEKRAN, 
1.2, 1 ); 
ScanLine( scena, XEKRAN, YEKRAN, 1.2, 1 ); 
break; 
// zmiana obiektu, na który patrzy obserwator 
case TAB:  clearviewport(); 


if( cel == NULL ) 
cel = scena->l _ axes; 
else 


cel = cel=>next; 
LookAtTarget( kamera, cel, 0, 0, 0 ); 
observer _view scene( scena, XEKRAN, YEKRAN, 
1.2, 1); 
DrawScene( scena ); 
break; 


// zamknięcie grafiki 
closegraph (); 
delete scene( scena ); 
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ROZDZIAŁ 6 


Środowisko MS-Windows 


6.1. Wstęp 


Środowisko graficzne Microsoft Windows zrealizowane jest jako nadbudowa 
systemu operacyjnego DOS i wraz z nim stanowi wielozadaniowy okienkowy 
system operacyjny, w którym każdemu programowi umożliwiono bezpośredni 
dostęp do 16MB pamięci chronionej, a ogółowi programów współbieżne wyko- 
nywanie z wykorzystaniem pamięci wirtualnej o pojemności do 64MB. 


Programy napisane pod środowisko Windows odznaczają się podobnym do 
siebie wyglądem oraz strukturą przyjmowanych przez program komend, co 
powoduje, iż ich obsługa jest często o wiele wygodniejsza i prostsza do nauki 
niż programów napisanych pod system DOS. Użytkownik może w prosty sposób 
przełączać aktualnie wykonywane zadania oraz wymieniać między nimi dane. 
Pod środowiskiem Windows można również uruchamiać programy konwencjo- 
nalne (napisane pod DOS), lecz wówczas nie można korzystać z szerokiej gamy 
opcji dostępnych dla programów Windows-owych. System Windows bazuje na 
tzw. graficznym interfejsie użytkownika — GUI (ang. Graphical User Interface), 
zwanym również popularnie graficznym środowiskiem okienkowym. GU] jest 
bardzo szeroko rozpowszechnionym interfejsem graficznym, zaimplementowa- 
nym na wielu typach komputerów — na komputerach kompatybilnych z IBM, 
lecz pracujących pod kontrolą systemu OS/2, jest to Presentation Manager, dla 
komputerów pracujących pod kontrolą systemu UNIX jest to X-Window, dla 
stacji Sun Microsystems — NeWS itp. 


6.2. Pojęcia podstawowe 


GUI pozwala na bardzo efektywne wykorzystanie wszelkich zasobów związa- 
nych z ekranem komputera. Ekran do tej pory, jako urządzenie pomocne 
użytkownikowi przy wprowadzaniu danych, spełniał funkcję bardzo skromną 
— był jedynie swojego rodzaju echem dla znaków wpisanych z klawiatury. 


GUI przemienia ekran w bardzo rozbudowany (co nie znaczy skomplikowany) 
interfejs, który umożliwia użytkownikowi wprowadzanie danych przy pomocy 
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klawiatury bądź innych urządzeń, np. myszki. Że względu na możliwość 
używania zarówno klawiatury, jak i myszki, musimy wprowadzić rozróżnienie 
słowa kursor. 


Kursorem myszki jest obiekt, który przemieszczany jest po ekranie za pomocą 
myszki, zaś kursorem klawiatury jest znacznik w tekście przemieszczany przy 
pomocy klawiatury, informujący nas o aktualnym położeniu w tekście. Ekran 
wyświetla różne obiekty graficzne, takie jak ikony, okienka, buttony czy 
scrollbary (terminy te są objaśnione w dalszej części tego rozdziału). 


Chcielibyśmy teraz wyjaśnić podstawowe pojęcia, których będziemy używać w 
dalszych częściach naszej książki, zwłaszcza w rozdziale będącym instrukcją 
obsługi naszego programu. Pojęcia te odnoszą się również do każdego innego 
programu napisanego pod system MS-Windows, dlatego iż określają one obie- 
kty będące fundamentalnymi składnikami każdego z tych programów. Używa- 
my nazw spolszczonych lub wręcz oryginalnych (angielskich), gdyż polskie 
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odpowiedniki po prostu nie istnieją, a wprowadzanie własnych odpowiedników 
na potrzeby książki uważamy za bezcelowe. 


Okienko (ang. window) jest podstawowym obiektem służącym do komunikacji 
pomiędzy użytkownikiem a systemem dla każdej aplikacji napisanej dla Win- 
dows i jedyną drogą dostępu aplikacji do urządzenia wyświetlającego. 


Okienko zadania można zredukować do ikony, którą za pomocą myszki można 
przemieszczać po ekranie i z powrotem przekształcić w okienko. MS-Windows 
jest systemem wielozadaniowym, więc w danej chwili może być aktywne więcej 
niż jedno zadanie, stąd na ekranie może znajdować się wiele ikon zadań 
wykonywanych współbieżnie. 


Okienko jest prostokątnym obszarem na ekranie i jest identyfikowane poprzez 
caption bar. Oknem nazywamy okienko powiększone do maksymalnego dopu- 
szczalnego rozmiaru. Wszystkie zasady posługiwania się okienkami dotyczą 
również i okien. 


Na rys.6.1 przedstawione jest okienko główne naszego programu z zaznaczo- 
nymi omawianymi obiektami. 


Każde okienko zawiera pole systemowe (ang. system menu box ), po wybraniu 
którego pojawia się menu opcji tzw. systemowych, umożliwiających podstawo- 
we operacje na okienku, takie jak zamknięcie okienka, zwiększenie lub zmniej- 
szenie jego rozmiarów, przełączenie na inne aktualnie wykonujące się zadanie 
itp. 


Za pomocą belki tytułowej (ang. caption bar ), która zawiera nazwę zadania, 
można okienko przemieszczać po ekranie. 


Większość operacji wykonywanych w okienkach jest realizowana za pomocą 
menu. Pola menu (ang. menu item ) są usytuowane bezpośrednio pod caption 
bar w belce menu (ang. menu bar )i wybranie ich powoduje wyświetlenie opcji. 
Wybranie opcji wyświetlonej jako przyciemniona nie jest w danym kontekście 
możliwe — taką opcję nazywać będziemy opcją zabronioną. Wybranie opcji 
zakończonej wielokropkiem powoduje wyświetlenie okienka dialogowego (ang. 
dialog box ). Wybranie opcji zakończonej trójkątem powoduje wyświetlenie 
opcji niższego rzędu, tzw. podmenu (ang. popup menu ). Opcje, które mogą być 
odhaczane, działają jak przełączniki. Każde wybranie takiej opcji powoduje 
przełączenie jej ze stanu wyłączona do stanu włączona i odwrotnie. 


W prawym górnym narożniku okienka głównego są wyświetlane dwa trójkąty 
(ang. minimize i maximize box ). Naciśnięcie myszką na trójkąt skierowany ku 
dołowi (minimize box ) powoduje zredukowanie okienka do ikony, a skierowa- 
nego ku górze (maximize box ) przekształcenie okienka w okno, czyli wyświet- 
lenie okienka na całym ekranie. 
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Dialog box jest wyświetlany w celu dokonania rozstrzygnięć; jest to tymczasowe 
okienko, które pozwala na wybór większej ilości rozkazów zgrupowanych pod 
jedną opcją menu. Dialog box może zawierać ramki tekstowe (ang. edittext ), 
ramki listowe (ang. listbox ), klucze (ang. checkbox ), przyciski (ang. radiobut- 
ton) i komendy (ang. button). Do ramki tekstowej wprowadza się tekst z kla- 
wiatury, z ramki listowej dokonuje się wyboru elementów, klucze i przyciski 
można ustawiać, a komendy wykonywać (komenda może wywołać następny 
dialog box). 
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Surface Scrollbary 


WEED m EYETI TĄ 
Green DOZ ULAMA LAŻŹA ŻŻ AAA PA 4 


WWO BWŚS BANIE Checkb 
ZZZZZŻŻŻŻŹNAŻAŁĄJN SĄ 225 eckbox 


Sheen BRIADININ MAŁ Z 


ZY, AMO sara UZ "em Ż 


Blink |SŻŻŻŻ 


zza 
Pd Paa, 
SZARA ZAŻŁ A Ź Ź 

GA UZŻEŻ 


a 


RZA PAULO 2, ; 
BŻ AŁMA a ŻA A M Z 2 


Phong 


Rys. 6.4 


Radiobuttony i checkboxy występują najczęściej w zestawach (ang. groupbox ). 
Różnią się od siebie tylko tym, że w danym groupboxie można ustawić dowolną 
liczbę checkboxów, ale co najwyżej jeden radiobutton. 


Na rysunkach 6.2, 6.3, 6.4 pokazano dialog boxy z naszewo programu wraz z 
zaznaczonymi omawianymi obiektami. 


Niektóre okienka są wyposażone w belki przewijania (ang. scrollbar ). Ma to 
miejsce wówczas, gdy okienko jest za małe do pomieszczenia całego tekstu lub 
rysunku. 


Scrollbary mogą być pionowe albo poziome. Każdy jest zakończony parą 
strzałek skierowanych poza niego. Każde naciśnięcie lewego przycisku myszy 
na jedną z nich powoduje zmianę położenia ruchomego obiektu wewnątrz 
scrollbara o jedną jednostkę (zwiększenie lub zmniejszenie — w zależności od 
kierunku strzałki). Na rys.6.3 występuje scrolibar pionowy służący do przeglą- 
dania zawartości listboxa. Niekiedy scrolibary są używane w celu ustawienia 
wartości parametru. Taki sposób ich wykorzystania przedstawiono na rys.6.4, 
gdzie za pomocą poziomych scrolibarów można zmieniać wartości Red, Green 
itp. 


Najprostszym w użyciu (i w wyglądzie) dialog boxem jest message box pojawia- 
jący się wówczas, gdy program potrzebuje prostej odpowiedzi od użytkownika 
(np. Tak lub Nie). Message box zawiera caption bar, system menu box, 1a ogół 
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jedną lub kilka linijek tekstu, jeden lub kilka (dwa lub trzy) buttony np. OK, 
OK i Cancel, Yes i No itp. 


SaveObject 


File already exists. 
Qverwrite ? 


Rys. 6.5 


Na rys.6.5 przedstawiono dwa przykładowe message boxy z naszego programu. 


Pierwszy z nich żąda od użytkownika potwierdzenia wcześniej wydanego 
polecenia — message box ma charakter pytający, drugi - informuje użytkownika 
o zdarzeniu i również żąda potwierdzenia. 


Przemieszczanie się w dialog boxie między opisanymi obiektami może odbywać 
się za pomocą myszki albo z klawiatury. Przechodzenie z jednej do drugiej 
grupy obiektów odbywa się bądź przez naciśnięcie w obrębie grupy lewego 
przycisku myszy, bądź przez naciśnięcie klawisza Tab. W obrębie jednej grupy 
można poruszać się używając kursorów klawiatury, natomiast uzyskać wybra- 
nie obiektu, naciskając klawisz Space. Większość dialog boxów jest wyposa- 
żona w buttony OKi Cancel. Wydanie komendy OK powoduje zaakceptowanie 
i zakończenie dialogu. Wydanie komendy Cancel (lub naciśnięcie klawisza 
Esc) powoduje odrzucenie i zakończenie dialogu. 


Jeśli nazwa buttona, radiobuttona czy podobnego obiektu zawiera podkreśloną 
literę (np. na rys.6.3 button Open), to najprostszym sposobem wydania komen- 
dy jest naciśnięcie klawisza Alt wraz z tą literą. 


6.3. Tryby przetwarzania 


Podstawą dobrego wielozadaniowego systemu operacyjnego jest dobra gospo- 
darka pamięcią operacyjną. Opisywanie algorytmów i struktur danych używa- 
nych przy zarządzaniu pamięcią w systemie MS-Windows jest tutaj całkowicie 
zbędne, stąd ograniczymy się do przekazania niezbędnego minimum wiedzy 
potrzebnego do zrozumienia różnych trybów pracy systemu. Jest to o tyle 
ważne, gdyż od trybu przetwarzania, w którym Windows aktualnie pracują, 
zależy m.in. szybkość działania naszego programu oraz możliwość tworzenia 
scen trójwymiarowych o dużym stopniu komplikacji. 
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Używając terminu Windows mamy na myśli środowisko MS-Windows wersja 
3.0 bądź późniejsza. 


Windows możemy uruchomić w jednym z trzech następujących trybów: 


1) real mode — na komputerach opartych na procesorze 8086 (lub 80286, 80386 
z pamięcią operacyjną mniejszą niż 1MB). Aby móc wykorzystać niektóre 
chociaż z zalet Windows, zalecane jest w tym trybie stosowanie driver'a 
i odpowiedniej do niego płyty pamięci (LIM EMS 4.0). 


2) standard mode — na komputerach opartych na procesorze 80286 z pamięcią 
przynajmniej 1MB (lub 80386 z pamięcią operacyjną mniejszą niż 2MB). Jest 
to tryb 286-protected procesora. Windows może używać w tym trybie do 16MB 
pamięci konwencjonalnej i rozszerzonej (ang. extended ). 


3) 386 enhanced mode — na komputerach opartych na procesorze 80386 z 
pamięcią przynajmniej 2MB. Ten tryb jest najlepszy, ale i najbardziej rozbu- 
dowany. Windows używa rejestrów procesora służących do stronicowania 
pamięci w celu realizacji pamięci wirtualnej. Niepotrzebne w danej chwili 
strony (o stałej długości 4KB) Windows może przechowywać na dysku i potem 
z powrotem ładować je do pamięci operacyjnej. Ponadto w trybie tym realizo- 
wane jest wielozadaniowe wykonywanie programów klasy DOS (ze względu na 
korzystanie z trybu Virtual-86 procesora 80386). 


W opisanych powyżej trybach wszystkie uwagi odnoszące się do procesora 
80386 dotyczą również procesora 80486. 


W naszym programie, ze względu na możliwość tworzenia skomplikowanych 
brył, jak i rozbudowanych scen, konieczna staje się odpowiednia gospodarka 
pamięcią operacyjną, pozwalająca na jak najbardziej efektywne jej wykorzy- 
stanie, a także umożliwiająca maksymalnie szybkie działanie programu, stąd 
tryb przez nas zalecany to 386 Enhanced Mode. 
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ROZDZIAŁ 7 


Instrukcja użytkowania 
programu RAY 


/.1 Wprowadzenie 


Po aktywowaniu programu RAY3D.EXE na ekranie ukazuje się okno główne 
naszego programu, które pokazano na rys.7.1. W celu zmniejszenia obszaru 
zajmowanego przez program na ekranie można nacisnąć minimize box, tzn. 
przekształcić okno w okienko (terminy te zostały wyjaśnione wcześniej w 
rozdziale 6). 


Okienko główne podzielone jest na cztery mniejsze okienka — na rysunku 
oznaczone jako okno O0,...,okno 3. 


W oknie 0 rysowana jest scena — wszystkie bryły znajdujące się na niej, bez 
symboli lamp oraz obserwatora. Scena rysowana jest jako rzut perspektywicz- 
ny na płaszczyznę rzutowania w układzie obserwatora (terminy te zostały 
objaśnione w rozdziale 4). Wszystkie bryły są przedstawione w postaci szkie- 
letowej, bez usuwania niewidocznych linii. 


W górnej części okna wyświetlane są wartości współrzędnych globalnego 
układu odniesienia. Współrzędne te wyświetlane są na bieżąco, tzn. podczas 
ruchu myszką nad aktywnym jednym z okien 1,2 lub 3. 


W oknach 1,2 i 3 rysowane są rzuty sceny na płaszczyzny rzutowania odpowie- 
dnio: XOZ, XOY i ZOY. Każde z tych okien wyświetla odpowiednie osie 
globalnego układu odniesienia, które wyznaczają jednoznacznie płaszczyznę 
rzutowania dla danego okna. Każde z tych okien rysuje obserwatora oraz 
lampy, a także osie lokalnych układów odniesienia związanych z konkretną 


bryłą. 
Wszelkie operacje na scenie (przemieszczanie obiektów, poruszanie się po 
globalnym układzie odniesienia, wybieranie obiektu itp.) można wykonywać 


właśnie na tych oknach, a mówiąc ściślej — na jednym z nich, tzw. oknie 
aktywnym. Okno aktywne w prawym górnym rogu wyświetla napis Active, zaś 
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Okno 1 Okno 0 Lampa 


- DZ 
GĘ 


= RAY 


ile Oplect Stęne Screen Render | | 


Bryła wraz z osiami Okno 2 Obserwator 
lokalnego układu 


Rys. 7.1 


nieaktywne — Inactive. Na rysunku oknem aktywnym jest okno 3, zaś okna 1 
i2 są nieaktywne. Wybieranie okna na aktywne następuje przez naciśnięcie 
lewego przycisku myszy w obrębie tego okna lub przez wciskanie klawisza Tab. 


Na rysunku zaznaczono symbole lampy, obserwatora, numery poszczególnych 
okien oraz bryły i osie jej układu lokalnego. 


Po aktywnym oknie można poruszać się wzdłuż osi globalnego układu odnie- 
sienia (klawisze kursora klawiatury) oraz przybliżać się lub oddalać od sceny 
(klawisze + oraz —). Prędkość przybliżania i oddalania jest zmienna i może być 
ustawiana przez użytkownika, co zostanie opisane w następnych punktach tego 
rozdziału. 
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7.2. Tworzenie brył oraz sceny. Operacje na scenie 


Na scenie można umieścić następujące bryły: 


kulę (ang. Sphere), 

pierścień (ang. ToruS), 
prostopadłościan (ang. Cube), 
płaszczyznę (ang. Plane), 
stożek (ang. Cone), 

walec (ang. Tube), 


figurę obrotową (ang. Rotate Figure). 


W celu dodania do sceny kuli należy z menu wywołać opcję Ob- 
ject Add I Sphere. Taki zapis należy rozumieć następująco: 


z menu głównego należy wywołać opcję Object, 
po rozwinięciu podmenu wywołać opcję Add, 


po rozwinięciu następnego podmenu wywołać opcję Sphere. 


Konwencji tej będziemy używali w całej książce, znak | oznacza więc, że opcja 
menu posiada własne menu. 


Add Sphere 


Name Radius 
Radlus 
Horizontal sections [6 | 


Vertical sections [6 | 


DO Rotate sections 


ANON NSŃ 
W $ 


: 
ODA COWOZ 


Rys. 7.2 
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Po wywołaniu opcji Object I Add | Sphere na ekranie ukazuje się dialog box 
przedstawiony na rys.7.2. 


Dialog box Add Sphere zawiera następujące obiekty: 
— edittext Name pozwalający wpisać nazwę dla tworzonej kuli, 
— edittext Radius pozwalający podać promień kuli, 


— edittext Horizontal Sections pozwalający podać ilość poziomych płatów, 
— edittext Vertical Sections pozwalający podać ilość pionowych płatów, 


— checkbox Rotate Sections — efekt uzyskiwany poprzez jego włączenie ilustru- 
je rys.7.3. Na rysunku przedstawiono dwa wycinki poziomych płatów po- 
wierzchni zbudowanych z trójkątów, których wzajemne położenie zależy od 


iii 


ustawienia cneckboxa, 


Checkbox wyłączony Checkbox włączony 


Rys. 7.3 


— naciśnięcie buttona Surface powoduje wyświetlenie dialog boxa Surface — 
zostanie to dokładnie opisane pod koniec tego punktu, gdyż dotyczy wszy- 
stkich tworzonych brył. 


Na rysunku obok dialog boxa narysowany jest przekrój przykładowej kuli z 
zaznaczonymi obiektami, które można zmieniać w tymże dialog boxie. 


ZAZNA 
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W celu dodania do sceny pierścienia należy z menu wywołać opcję Ob- 
ject I Add I Torus. 


Po wywołaniu opcji Object | Add ITorus na ekranie ukazuje się dialog box 
przedstawiony na rys.7.4. 


Add Torus 


Name RadiusY 


Radius X 


Radius Y 
Radius 


74 
74 
PA 


Radius « RadiusX 


Rys. 7.4 


Dialog box Add Torus zawiera następujące obiekty: 


— edittext Name pozwalający wpisać nazwę dla tworzonego pierścienia, 

— edittext RadiusX pozwalający podać poziomy promień przekroju pierścienia, 
— edittext RadiusY pozwalający podać pionowy promień przekroju pierścienia, 
— edittext Radius pozwalający podać promień pierścienia, 

— edittext Horizontal Sections pozwalający podać ilość poziomych płatów, 

— edittext Vertical Sections pozwalający podać ilość pionowych płatów, 

— działanie checkboxa Rotate Sections ilustruje rys.7.3, 

— naciśnięcie buttona Surface powoduje wyświetlenie dialog boxa Surface. 


Na rysunku obok dialog boxa narysowany jest przekrój przykładowego pier- 
ścienia z zaznaczonymi obiektami, które można zmieniać w tymże dialog boxie. 
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W celu dodania do sceny prostopadłościanu należy z menu wywołać opcję 
Object | Add I Cube. 


Po wywołaniu tej opcji na ekranie ukazuje się dialog box przedstawiony na 
ry8.7.0. 


Add Cube 


Height 


Width Depth 


Rys. 7.5 
Dialog box Add Cube zawiera następujące obiekty: 
— edittext Name pozwalający wpisać nazwę dla tworzonego prostopadłościa- 
nu, 
— edittext Width pozwalający podać szerokość prostopadłościanu, 
— edittext Depth pozwalający podać głębokość prostopadłościanu, 
— edittext Height pozwalający podać wysokość prostopadłościanu, 
— naciśnięcie buttona Surface powoduje wyświetlenie dialog boxa Surface. 


Na rysunku obok dialog boxa narysowany jest przykładowy prostopadłościan 
z zaznaczonymi obiektami, które można zmieniać w tymże dialog boxie. 


W celu dodania do sceny płaszczyzny należy z menu wywołać opcję Ob- 
jectl Add | Plane.Po wywołaniu tej opcji na ekranie ukazuje się dialog box 
przedstawiony na rys.7.6. 


Dialog box Add Plane zawiera następujące obiekty: 


— edittext Name pozwalający wpisać nazwę dla tworzonej płaszczyzny, 
— edittext Width pozwalający podać szerokość płaszczyzny, 
— edittext Depth pozwalający podać głębokość płaszczyzny, 


— naciśnięcie buttona Surface powoduje wyświetlenie dialog boxa Surface. 
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Add Plane 
ŻA ORÓPOZŹ Ź TĄ KALPOZTZRAZZRZA Ż 


Rys. 7.6 


Na rysunku obok dialog boxa narysowana jest przykładowa płaszczyzna z 
zaznaczonymi obiektami, które można zmieniać w tymże dialog boxie. 


W celu dodania do sceny stożka należy z menu wywołać opcję Ob- 
ject | Add I Cone. 


Po wywołaniu tej opcji na ekranie ukazuje się dialog box przedstawiony na 
rys.7.7. 


Add Cone 
Name 
Aadius eight 
Height 


Sections 


e a W” m a” 
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Sections=3 


Rys. 7.7 


Dialog box Add Cone zawiera następujące obiekty: 

— edittext Name pozwalający wpisać nazwę dla tworzonego stożka, 
— edittext Radius pozwalający podać promień podstawy stożka, 

— edittext Height pozwalający podać wysokość stożka, 


Realistyczna grafika 3D 155 


- edittext Sections pozwalający podać ilość wierzchołków znajdujących się na 
podstawie stożka, 


- checkbox Bottom umożliwiający utworzenie stożka z otwartą podstawą 
(checkbox wyłączony) lub zamkniętą (włączony), 


— naciśnięcie buttona Surface powoduje wyświetlenie dialog boxa Surface. 
Na rysunku obok dialog boxa narysowany jest przykładowy stożek z zaznaczo- 
nymi obiektami, które można zmieniać w tymże dialog boxie (w tym przypadku 
jest to ostrosłup, gdyż Sections = 3). 


W celu dodania do sceny walca należy z menu wywołać opcję Object | Add I Tu- 
be. 


Po wywołaniu tej opcji na ekranie ukazuje się dialog box przedstawiony na 


rys.7.8. 
Add Tube RadiusTop 
Radius Bottom C r-2 
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Rys. 7.8 
Dialog box Add Tube zawiera następujące obiekty: 


— edittext Name pozwalający wpisać nazwę dla tworzonego walca, 

— edittext Radius Bottom pozwalający podać promień dolnej podstawy walca, 
— edittext Radius Top pozwalający podać promień górnej podstawy walca, 

— edittext Height pozwalający podać wysokość walca, 

— edittext Sections pozwalający podać ilość punktów podstawy walca, 
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— checkbox Bottom umożliwiający utworzenie walca z otwartą dolną podstawą 
(checkbox wyłączony) lub zamkniętą (włączony), 


— checkbox Top umożliwiający utworzenie walca z otwartą górną podstawą 
(checkbox wyłączony) lub zamkniętą (włączony), 


— działanie checkboxa Rotate Sections ilustruje rys.7.3, 
— naciśnięcie buttona Surface powoduje wyświetlenie dialog boxa Surface. 


Na rysunku obok dialog boxa narysowany jest przykładowy walec z zaznaczo- 
nymi obiektami, które można zmieniać w tymże dialog boxie. 


W celu dodania do sceny figury obrotowej należy z menu wywołać opcję 
Object I Add | Figure. 


Po wywołaniu tej opcji na ekranie ukazuje się dialog box przedstawiony na 
rys.7.9. 


Add Rotate Figure 


Name Seciions 


Oś obrotu DJ Bottom 
1] Iop 


ŻJ Rotate sections 


Rys. 7.9 


Dialog box Add Rotate Figure zawiera następujące obiekty: 


— okienko, w którym użytkownik definiuje krzywą, która po obróceniu wokół 
osi obrotu daje figurę obrotową. Krzywa powstaje z kolejnych punktów, 
które użytkownik dodaje poprzez naciśnięcie lewego przycisku myszy; 


— edittext Name pozwalający wpisać nazwę dla tworzonej figury obrotowej; 
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edittext Sections pozwalający podać ilość wierzchołków znajdujących się na 
podstawie figury; 


checkbox Bottom umożliwiający utworzenie figury z otwartą dolną podstawą 
(checkbox wyłączony) lub zamkniętą (włączony); 


checkbox Top umożliwiający utworzenie figury z otwartą górną podstawą 
(checkbox wyłączony) lub zamkniętą (włączony); 


działanie checkboxa Rotate Sections ilustruje rys.7.3; 


button Clear służący do wyczyszczenia zawartości okienka i wykasowania 
krzywej; 


button Undo służący do anulowania ostatniej akcji — po dwukrotnym użyciu 
Undo zmienia się na Reundo; 


button Delete służący do kasowania punktu; 


naciśnięcie buttona Surface powoduje wyświetlenie dialog boxa Surface. 


W każdym z dialog boxów do tworzenia bryły znajduje się button Surface. Jego 
naciśnięcie powoduje pojawienie się dialog boxa Surface przedstawionego na 
rys.7.10. 


Rys. 7.10 


Dialog box Surface umożliwia zdefiniowanie powierzchni aktualnie tworzonej 
bryły. Przy pomocy scrollbarów Red, Green i Blue można zdefiniować kolor 
powierzchni, który na bieżąco jest wyświetlany w okienku Kolor. Scrolibar 
Sheen służy do zmiany połyskliwości powierzchni od 0 (matowa) do 100% 
(zwierciadlana). Scrollbar Blink służy do zmiany rozmiarów odblasku światła 
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na powierzchni. Włączenie checkboxa Phong powoduje wygładzanie powierz- 
chni brył w trakcie generowania obrazu wyjściowego. 


Aby wykonać jakąkolwiek operację na obiekcie, należy najpierw go wybrać. 
Wybranie (selektowanie) bryły odbywa się na różne sposoby: 


— najprostszy sposób to naciśnięcie lewego przycisku myszy, gdy jej kursor 
znajduje się na środku lokalnego układu odniesienia bryły, którą chcemy 
wybrać, 

- można również wywołać opcję Object | Find, co zostanie opisane w dalszej 
części tego rozdziału, 

— wywołanie opcji Object I NextObject lub naciśnięcie kombinacji klawiszy 
Ctr1l+N powoduje wybranie następnej bryły na scenie; jeśli nie było dotąd 
żadnej wybranej bryły, to wybiera pierwszą utworzoną przez użytkownika. 
Opcja Object I NextObject jest zabroniona, gdy na scenie nie ma żadnej 
bryły. 

Wybranie lampy odbywa się przez naciśnięcie lewego przycisku myszy, gdy jej 

kursor znajduje się na środku kółka będącego symbolem lampy. 


Wybranie obserwatora odbywa się podobnie jak selektowanie lampy (operuje- 
my na symbolu obserwatora) lub poprzez wywołanie opcji Scene I SelectOb- 
server (opcja ta jest zabroniona, gdy obserwator jest już wybrany). 


Wybrane obiekty (w przypadku bryły również osie jej układu lokalnego) są 
rysowane innym kolorem (jaśniejszym), wszystkie inne obiekty rysowane są 
kolorem czarnym na białym tle. 


W danej chwili na scenie może być tylko jeden wybrany obiekt. 


Zaselektowany obiekt można dowolnie przemieszczać i orientować. W tym celu 
można wywołać opcję Object I Transform lub nacisnąć klawisze Ctrl+T (op- 
cja ta jest zabroniona, gdy na scenie nie ma żadnego wybranego obiektu). 
Wywołanie tej opcji powoduje wyświetlenie dialog boxa Transformation przed- 
stawionego na rys.7.11. 


Dialog box Transformation zawiera następujące obiekty: 

— groupbox Origin, a w nim: 
e radiobutton Global — wybiera globalny układ odniesienia, 
e radiobutton Local — wybiera lokalny układ odniesienia, 


— groupbox What, a w nim: 
e radiobutton Object — wybiera sam obiekt, 
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e radiobutton Axes — wybiera same osie, 

e radiobutton Object-Axes — wybiera obiekt wraz z osiami, 
następujące radiobuttony do wyboru transformacji: 

© radiobutton Translation — wybiera translację o wektor, 


e radiobutton Position — wybiera umieszczenie obiektu w punkcie 
o współrzędnych (X,Y,Z), 


e radiobutton Scale — wybiera przeskalowanie obiektu o wartości 
XY i Z wzdłuż odpowiadających im osi, 


e radiobuttony RotationX, RotationY i RotationZ — wybierają 
obrót wokół odpowiadających im osi, 


edittexty X, Yi Z — umożliwiają wpisanie wartości dla odpowiedniej trans- 
formacji. 


Transformation 


© Global 


O Local 


© Transłation © Rotation X 


O Position  ORotation Y 
O Scale O Rotation Z 


Rys. 7.11 


W zależności od rodzaju obiektu, na którym będzie wykonana transformacja 
oraz od ustawienia poszczególnych radiobuttonów w dialog boxie , możliwe są 
następujące kombinacje: 

Dla zaselektowanej bryły: 


radiobutton Position zabrania radiobutton Object, 
radiobutton Scale zabrania radiobuttony Axes i Object-Axes, 
radiobutton RotationX zabrania edittexty YiZ, 

radiobutton RotationY zabrania edittexty Xi Z, 

radiobutton RotationZ zabrania edittexty XiY. 
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Dla zaselektowanej lampy: 


- zabronione są następujące radiobuttony: Local, Axes, Object-Axes, Scale, 
RotationX, RotationY i RotationZ. 


Dla zaselektowanego obserwatora: 


— zabronione są radiobuttony: Object, Object-Axes i Scale, a ponadto: 
radiobutton RotationX zabrania edittexty Y i Z, 
— radiobutton RotationY zabrania edittexty Xi Z, 
— radiobutton RotationZ zabrania edittexty Xi Y. 


Ograniczenia te mają związek z przyjętymi przez nas zasadami w programie 
oraz są zgodne z ograniczeniami matematyki trójwymiarowej. Nie jest np. 
możliwe skalowanie osi lokalnego układu odniesienia, czy też obrót wokół osi 
X, gdy wybrało się obrót wokół osi Y. 


Przemieszczanie obiektu można także osiągnąć przy pomocy wyłącznie myszy. 
Gdy dany obiekt jest wybrany, należy nacisnąć prawy przycisk mysz. Pojawia 
się wówczas obwiednia obiektu, która w przypadku brył ma kształt prostokąt- 
nej ramki, zaś w przypadku lamp oraz obserwatora kształt kółka. Trzymając 
cały czas wciśnięty prawy przycisk myszy można przemieszczać obwiednię po 
aktywnym oknie. Zwolnienie prawego przycisku myszy wewnątrz aktywnego 
okienka jest równoznaczne z zaakceptowaniem transformacji. Opisana trans- 
formacja jest równoznaczna z następującym ustawieniem radiobuttonów w 
dialog boxie Transformation: 


— dla bryły: Global, Object-Axes, Position; 
— dla lampy: Global, Object, Position; 
— dla obserwatora: Global, Axes, Position. 


Wciśnięcie klawisza Space powoduje "centrowanie" sceny — polega to na tym, 
iż obiekt zasełektowany jest rysowany na środku okienek 1, 213. Nie ma w tym 
momencie miejsca żadna transformacja obiektów, a jedynie zmiana wyświet- 
lania sceny na ekranie (podobna sytuacja występuje wówczas, gdy użyjemy 
kursorów klawiatury w celu przemieszczenia się po okienku aktywnym). 


Zaselektowaną bryłę można skopiować wywołując opcję Object Copy lub 
naciskając klawisz Ins (opcja ta jest zabroniona, gdy nie ma żadnej wybranej 
bryły). Każde wywołanie tej opcji powoduje jednokrotne skopiowanie bryły. 


Wielokrotne kopiowanie bryły można osiągnąć poprzez wywołanie opcji Ob- 
ject I Replicate (opcja ta jest zabroniona, gdy nie ma żadnej wybranej bryły). 
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Wybranie tej opcji powoduje wyświetlenie dialog boxa Replicate przedstawio- 
nego na rys.7.12. 


Replicate 


Number 3_| 


Translation 
E Y z maj 


Rys. 7.12 


Dialog box Replicate zawiera następujące obiekty: 


— edittext Number pozwalający wpisać ilość kopii bryły; 


— groupbox Translation, a w nim edittexty X, Y i Z, gdzie użytkownik może 
wpisać wartości, które składają się na wektor [X,Y,Z], o który będzie 
przesuwana kolejna bryła podczas powielania. 


Wybraną bryłę bądź lampę można skasować ze sceny wywołując opcje Ob- 
ject | Delete dla bryły (opcja zabroniona, gdy nie ma zaselektowanej bryły) 
oraz Scene | DeleteLamp dla lampy (opcja zabroniona, gdy nie ma zasełekto- 
wanej lampy). Można również użyć klawisza Del (automatycznie rozpoznaje 
czy kasowanie ma dotyczyć bryły, czy lampy). Przed skasowaniem danego 
obiektu ze sceny program prosi o potwierdzenie, wyświetlając stosowny mes- 
sage box. Przykładowe message boxy dla kasowania bryły i kasowania lampy 
przedstawiono na rys.7.13. 


Delete Delete 


© Delete this object ? © Delete this lamp ? 


Rys. 7.13 


W celu znalezienia bryły na scenie należy wywołać opcję Object IFind lub 
nacisnąć klawisze Ctrl+F (opcja jest zabroniona, gdy na scenie nie ma żadnej 
bryły). 
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Wywołanie tej opcji powoduje wyświetlenie dialog boxa FindObject przedsta- 
wionego na rys.7.14. 


Find Object: 


axes7.plane2 
axes6.conel 
axes$5.figure6 


 OALSCNĄ ŚŚ 


AGCA 


Rys. 7.14 


Dialog box FindObject zawiera listbox wraz z pionowym scrollbarem oraz 
edittext. W listboxie pojawiają się nazwy wszystkich brył i osi ich lokalnych 
układów odniesienia obecnych na scenie. Nazwy te są nieco inne od nazw 
podanych przez użytkownika w momencie tworzenia bryły, gdyż zawierają 
ponadto numery osi. Np. podając nazwę dla kuli: Kulal faktycznie otrzyma ona 
nazwę Axes***,Kulal (*** oznacza kolejny numer osi na scenie). Otrzymana 
w ten sposób nazwa obiektu jest unikatowa i na tej samej scenie nigdy już nie 
wystąpi; zwalnia to użytkownika od podawania za każdym razem innej nazwy 
dla tworzonej bryły, np. podczas powielania. 


Nazwa wybranej bryły wyświetlana jest w edittext. Naciśnięcie buttona Find 
powoduje znalezienie bryły oraz zaznaczenie jej jako obiektu zaselektowanego. 


Zaselektowaną bryłę można nagrać na dysk wywołując opcję File |SaveOb- 
ject (opcja jest zabroniona, gdy na scenie nie ma wybranej bryły). Na ekranie 
pojawia się wówczas dialog box przedstawiony na rys.7.15. 


Save File Name As: 


TUBE.OBT 


Current Directory: cĄborlandcjexe 


Rys. 7.15 


Dialog box zawiera edittext, w którym należy podać nazwę nagrywanej bryły, 
jako domyślne rozszerzenie stosuje się OBT. Program nagrywa bryłę w bieżą- 
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cym katalogu. Można go zmienić w dialog boxach opisanych dalej, a służących 
do wgrywania brył bądź scen. Button OK może być zabroniony w sytuacji, gdy 
edittext nie zawiera żadnego znaku. 


W przypadku gdy użytkownik poda samą nazwę (bez rozszerzenia), program 
sam dołączy domyślne rozszerzenie OBT. W sytuacji gdy plik o takiej nazwie 
już na dysku istnieje, program wyświetla message box przedstawiony na 
rys.7.16, żądając potwierdzenia zapisu. 


SaveObject 


AD File already exists. 
Overwrite ? 


Rys. 7.16 


Nagraną wcześniej na dysk bryłę można wgrać do programu, czyli dołączyć do 
sceny wywołując opcję File | LoadObject. Na ekranie pokazuje się wówczas 
dialog box przedstawiony na rys.7.17. 


Open File Name: 


Files in cĄborlandcjexe 


Rys. 7.17 


Dialog box zawiera listbox, w którym wypisane są nazwy wszystkich plików z 
rozszerzeniem OBT w bieżącym katalogu. Wybrana nazwa pliku zostaje wy- 
świetlona w edittext. W nim można również podać bezpośrednio nazwę pliku — 
domyślnym rozszerzeniem jest OBT. Dialog box umożliwia również zmianę 
bieżącego katalogu oraz napędu dyskietek (odpowiednie wybranie pozycji w 
listboxie np. [-a-] zmieni napęd dyskietek na a:, zaś [..] powoduje "jeden krok 
wyżej" w ścieżce na dysku). 
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Naciśnięcie buttona Open powoduje załadowanie bryły do pamięci (o ile plik 
taki istnieje na dysku) i automatyczne dołączenie jej do sceny. Jeśli nazwa 
podana przez użytkownika jest bez rozszerzenia, wówczas program sam "do- 
kleja" rozszerzenie OBT. Button Open może być zabroniony w sytuacji, gdy 
edittext nie zawiera żadnego znaku. 


Dwukrotne szybkie naciśnięcie lewego przycisku myszy nad aktywnym oknem 
powoduje wyświetlenie dialog boxa Objectlnfo, zawierającego podstawowe 
informacje o zaselektowanej bryle. Jeśli takiej nie ma, to podwójne wciśnięcie 
przycisku myszy nie spowoduje żadnej akcji. Rys.7.18 przedstawia dialog box 
Objectlnfo dla przykładowej bryły — kuli. 


Object INFO 


Name: axes1l.Spherel 
Numberwverices: 134 
Numberedges: 396 
Number faces: 264 


Rys. 7.18 
Dialog box Objectlnfo zawiera następujące informacje o zaselektowanej bryle: 


— Name — nazwa bryły, 

— Number vertices — liczba wierzchołków bryły, 

— Number edges — liczba krawędzi bryły, 

— Number faces — liczba ścian bryły. 

Naciśnięcie buttona Surface powoduje wywołanie dialog boxa Surface, umo- 
żliwiając w ten sposób użytkownikowi zmianę parametrów powierzchni bryły 
utworzonej wcześniej. 

Do sceny można dodać lampę o dowolnej barwie i w dowolnym miejscu, 


wywołując opcję Scene | AddLamp lub naciskając klawisze Ctrl+L. Jej wy- 
wołanie powoduje wyświetlenie dialog boxa AddLamp przedstawionego na 


rys.7.19. 
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Add Lamp 


Rys. 7.19 
Dialog box AddLamp zawiera następujące obiekty: 


— edittext Name, pozwalający wpisać nazwę dla tworzonej lampy, 


— trzyscrollbary Red, Green i Blue, służące do ustawienia koloru lampy, który 
na bieżąco wyświetlany jest w okienku Kolor, 


— groupbox Position, a w nim trzy edittexty X, Y i Z, służące do wpisania 


wartości odpowiadających umieszczeniu lampy w punkcie (X,Y,Z) globalne- 
go układu odniesienia. 


Observer 


O Selected Object: axes1.cube0 


© Current: axes7.plane2 


166 M. Domaradzki, R. Gembara 


Do zmiany parametrów obserwatora, kierunku patrzenia itp. służy dialog box 
Observer, który można wywołać poprzez opcję Scene I Observer lub naciska- 
jąc klawisze Ctr1+0. Dialog box Observer jest przedstawiony na rys.7.20. 


LEI (32 


mały kęt patrzenia duży kął patrzenia 


Rys. 7.21 
Dialog box Observer zawiera następujące obiekty: 


— groupbox View, a w nim scrolibar Angle służący do zmiany kąta patrzenia 
obserwatora. Rys.7.21 przedstawia prostą scenę z małym oraz dużym kątem 
patrzenia. 


— groupbox Target, a w nim: 


© radiobutton SelectedObject, obok którego wyświetlana jest 
nazwa zaselektowanej bryły. Po wybraniu tego radiobuttona 
obserwator będzie patrzył na bryłę, której nazwa umieszczona 
jest obok. Radiobutton SelectedObject jest zabroniony, gdy na 
scenie nie ma zaselektowanej bryły; 


© radiobutton Current, obok którego wyświetlana jest nazwa 
bryły, na którą patrzył do tej pory obserwator. Po wybraniu tego 
radiobuttona obserwator będzie nadal patrzył na tę bryłę 
(sytuacja ma miejsce wówczas, gdy użytkownik zmienia tylko 
kąt patrzenia obserwatora bez zmiany obiektu patrzenia). 
Radiobutton Current jest zabroniony, gdy obserwator nie 
patrzył dotąd na żadną bryłę; 

e radiobutton None — jego wybranie powoduje, iż obserwator nie 
będzie patrzył na żadną konkretną bryłę. 


Jeśli obserwator patrzy na konkretną bryłę, to każde przemieszczenie obser- 
watora lub tej bryły powoduje zmianę kierunku patrzenia obserwatora tak, by 
zawsze patrzył na tę bryłę. Np. na rys.7.1 na początku rozdziału wyraźnie 
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widać, że obserwator patrzy na prostopadłościan (odcinek w symbolu obserwa- 
tora jest skierowany na środek lokalnego układu odniesienia bryły). 


Wywołanie opcji Scene | UnselectAl1 powoduje, iż zaselektowany obiekt (bry- 
ła, lampa bądź obserwator) przestaje być obiektem wybranym. Opcja ta jest 
zabroniona, gdy na scenie nie ma żadnego wybranego obiektu. 


Całą scenę, a więc bryły, lampy i ich wzajemne rozmieszczenie można skaso- 
wać, wywołując opcję Scene | New. Wykonanie tej operacji poprzedzone jest 
message boxem przedstawionym na rys.7.22, dając w ten spsób użytkownikowi 
możliwość anulowania wywołanej opcji. 


Rys. 7.22 


Wywołanie opisanej wyżej opcji powoduje powrót do początkowych wartości — 
obserwator w punkcie (0,0,-400) itp. 


Całą scenę można nagrać na dysk wywołując opcję File I SaveScene (opcja jest 


zabroniona, gdy na scenie niema żadnej bryły). Na ekranie pojawia się wówczas 
dialog box przedstawiony na rys.7.23. 


Save File Name As: 


FIGS.SCN 


Current Directory: cHborlandcjexejscn 


Rys. 7.23 


Dialog box zawiera edittext, w którym należy podać nazwę nagrywanej sceny; 
jako domyślne rozszerzenie stosuje się SCN. Program nagrywa scenę w bieżą- 
cym katalogu. Button OK może być zabroniony w sytuacji, gdy edittext nie 
zawiera żadnego znaku. W przypadku gdy użytkownik poda samą nazwę (bez 
rozszerzenia), program sam dołączy domyślne rozszerzenie SCN. 


Podobnie jak w przypadku nagrywania bryły, w sytuacji gdy plik o takiej 
nazwie już na dysku istnieje, program wyświetla message box przedstawiony 
na rys.7.24, żądając potwierdzenia zapisu. 
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SaveScene 


File already exists. 
Overwrite ? 


Rys. 7.24 
Nagraną wcześniej na dysk scenę można wgrać do programu i później np. 


modyfikować ją, wywołując opcję File l LoadScene. Na ekranie pokazuje się 
wówczas dialog box przedstawiony na rys.7.25. 


Open File Name: 


allfigs.scn 


Files in cĄborlandcjexejscn 


Rys. 7.25 


Obsługa tego dialog boxa jest analogiczna do obsługi dialog boxa dla opcji 
File | LoadObject opisanej wcześniej w tym rozdziale. W tym przypadku 
domyślnym rozszerzeniem plików jest SCN. 


Jeśli użytkownik postanowi wczytać z dysku nową scenę, a jednocześnie w 
programie jest już jedna scena, wówczas program wyświetla message box 
identyczny jak dla opcji Scene I New, żądając od użytkownika potwierdzenia 
skasowania starej sceny. 


Kolor tła na scenie można zmieniać, wywołując opcję Scene I Sky. Pojawia się 
wówczas dialog box przedstawiony na rys.7.26. 


Zawiera on groupbox Color, a w nim trzy scrolibary Red, Green i Blue, które, 
podobnie jak w przypadku opisanych wcześniej dialog boxów Surface i Ad- 
dLamp, służą do ustawienia koloru. Zmieniany kolor tła jest na bieżąco 
wyświetlany w okienku Kolor — początkowo ustawiony jest kolor szary. 


Ustawiony w ten sposób kolor tła jest widoczny podczas generowania obrazu 
wyjściowego, co jest tematem następnego punktu. 
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Kolor 


Rys. 7.26 


Prędkość przybliżania się do sceny (klawisz +) oraz prędkość oddalania (kla- 
wisz —) jest regulowana i może być ustawiona przez użytkownika poprzez 
wywołanie opcji Screen IScale. Wyświetla się wówczas dialog box Scale 
przedstawiony na rys.7.27, zawierający ediitext, gdzie użytkownik może podać 
nową wartość. Początkowo Scale ustawione jest na 1.4. 


Rys. 7.27 


W celu odświeżenia ekranu (narysowanie całej sceny we wszystkich oknach od 
nowa) można wywołać opcję Screen | Redraw lub nacisnąć klawisze Ctrl+W. 
Użycie tej opcji jest konieczne po wielu transformacjach obiektów, gdyż niektó- 
re obiekty mogły zostać częściowo zamazane przez inne (jest to wynik zasto- 
sowanego algorytmu rysowania obiektów — w celu maksymalnego przyśpiesze- 
nia rysowania, odnawiane są tylko te obiekty, na których operuje w danej chwili 
użytkownik). 


7.3. Generowanie obrazu wyjściowego 


Parametry obrazu wyjściowego można ustawić wywołując opcję Render |Pic- 
tureSetup. Wyświetlany jest wówczas dialog box PictureSetup przedstawiony 
na rys.7.28. 
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Picture Setup 


O 640 x 350 


O 800 x 600 O 320x200 


© | 640 x 480] O User defined 


Rys. 7.28 
Dialog box PictureSetup zawiera następujące obiekty: 


— groupbox Mode, a w nim sześć radiobuttonów służących do wybrania 
rozdzielczości (w pikslach) obrazu wyjściowego. Poszczególne radiobuttony 
mogą być zabronione w sytuacji, gdy rozdzielczość przez nie proponowana 
przewyższa aktualną rozdzielczość MS-Windows (np. na rysunku zabronio- 
ny jest radiobutton 1024x768, gdyż Windows zostało uruchomione w roz- 
dzielczości 800x600); 


— groupbox Screen, a w nim dwaedittexty Width i Height, w których wyświet- 
lana jest odpowiednio szerokość i wysokość w pikslach obrazu wyjściowego; 


— groupbox Pixel, a w nim dwa edittexty XAspect i YAspect, w których 
wyświetlany jest stosunek boków piksla. 


W przypadku wybrania jednego z radiobuttonów groupboxa Mode, za wyjąt- 
kiem UserDefined, wszystkie edittexty są zabronione. Jest to oczywiste, gdyż 
użytkownik wybierając np. rozdzielczość 640x480, automatycznie otrzymuje 
wartość Width równą 640 oraz Height — 480, więc nie miałoby sensu umożli- 
wienie zmian tych wartości. Wartości Xaspect i YAspect są wyświetlane w 
zależności od trybu pracy danej karty graficznej i nie są zależne od rozdziel- 
czości (np. dla karty SVGA pracującej w trybie VGA są zawsze równe lil, 
różne są natomiast dla tej samej karty, ale pracującej np. w trybie CGA). 
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Wybór ostatniego radiobuttona z groupboxa Mode, a więc UserDefined daje 
użytkownikowi możliwość ustawienia dowolnych wartości dla wymiarów obra- 
zu wyjściowego (nie większych jednak niż dostępne aktualnie w Windows) oraz 
dowolnych wartości określających stosunek boków piksla. W tym jedynym 
przypadku wszystkie edittexty są aktywne. 


Domyślną wartością rozdzielczości jest 320x200 (aspekty zależą od konkretnej 
karty). 


Wywołanie opcji Render | Start uruchamia proces generowania obrazu wyj- 
ściowego. Opcja jest zabroniona, gdy na scenie nie ma żadnej bryły. 


Obraz generowany jest na oddzielnym ekranie, na tle, którego kolor ustawiony 
został w dialog boxie Sky opisanym wcześniej. Przed wyświetleniem obrazu 
końcowego niezbędne są odpowiednie przeliczenia, stąd na ekranie pojawia się 
napis "Calculating. Please wait.' oraz podawany jest (w procentach) stopień 
przeliczenia. Gdy stopień ten osiągnie wartość 100%, na ekranie pojawia się 
obraz wyjściowy. Liczenie obrazu wyjściowego można przerwać, wciskając 
klawisz Esc (reakcja nań może nie być natychmiastowa), program wróci 
wówczas do okna głównego. Gdy rysowanie obrazu wyjściowego zakończy się 
(znikają napisy i kursor myszy), do okna głównego można wrócić, naciskając 
również klawisz Esc, lecz wtedy program żąda potwierdzenia, wyświetlając 
nessage box przedstawiony na rys.7.29. 


() Return to Main Window ? 


Rys. 7.29 
"Odhaczenie" opcji Render |ToFile powoduje, że obraz wyjściowy podczas 


procesu liczenia (przed wyświetleniem) będzie jednocześnie zapisywany do 
pliku w formacie TIFF (jest to format 24 bitowy). Nazwę tego pliku można 


Save File Name As: 


SCREEN.TIF 


Current Directory: cĄborlandcjexe 


Rys. 7.30 
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podać, wywołując opcję File I Export. Wywołanie tej opcji powoduje wyświet- 
lenie dialog boxa przedstawionego na rys.7.30. 


Obsługa tego dialog boxa jest analogiczna do obsługi dialog boxów SaveScene 
oraz SaveObject opisanych wcześniej. Domyślnym rozszerzeniem dla plików 
jest TIF (standard dla formatu TIFF dla IBM PO). 


Początkowo w programie nie ma żadnej nazwy dla ekranu wyjściowego w 
formacie TIFF. "Odhaczenie" opcji Render |ToFile, jeśli nie ma nazwy plik 
TIFF, powoduje automatyczne wywołanie dialog boxa Export. 


Do wygenerowania obrazu wyjściowego niezbędne jest minimum 256 kolorów. 
W sytuacji gdy użytkownik wywoła opcję Render I Start przy Windows uru- 
chomionych np. w 16 kolorach, na ekranie pojawia się message box przedsta- 
wiony na rys.7.31, wymuszający powrót do okna głównego. 


Program requires min 256 colors 
to generate Output Picture. 


Retum to Main Window. 


7.4. Zakończenie pracy z programem 


W celu zakończenia pracy z programem należy wywołać opcję File | Exit lub 
nacisnąć klawisze Ctrl+X albo użyć menu systemowego (opcja Close). Jeśli na 
scenie nie ma żadnej bryły, program dezaktywuje się natychmiast, w przeciw- 
nym razie żąda potwierdzenia, wyświetlając message box przedstawiony na 
rys.7.82. Odpowiedź pozytywna kończy program. 
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Gdy program pracuje w tle (jako zadanie drugoplanowe) i zostanie wywołana 
opcja zamykająca nadbudowę MS-Windows, wówczas pojawia się message box 
przedstawiony na rys.7.33. 


Odpowiedź pozytywna kończy pracę z programem i nadbudową MS-Windows. 


I'm still active ! 
Close ? 
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ROZDZIAŁ 8 


Przykłady tworzenia scen 


8.1. Przykład 1 


Na początek rozważmy budowę prostej sceny przedstawionej na rys.8.1.1. 


Scena składa się z pięciu brył, z których my zdefiniowaliśmy tylko trzy; kula 
oraz walec są skopiowane. 


WNCYAT 


„dr. 


Rys. 8.1.1 


ZZ EO Ą 


Realistyczna grafika 3D 175 


Pierwszą bryłą, jaką tworzymy, jest sześcian — wywołujemy z menu opcję 
Object | Add I Cube i wpisujemy parametry jak poniżej: 


Name: szescian 
width: 40 
Depth: 40 
Height: 40 


Następnie tworzymy walec, czyli wywołujemy opcję Object I Add I Tube i wpi- 
sujemy wartości: 


Name: walec 

Radius Bottom: 5 

Radius Top: 5 

Height: 40 

Sections: 10 

Bottom — nie ustawiony 
Top — nie ustawiony 


Rotate Sections — ustawiony 


Otrzymany w ten sposób walec jest figurą "świeżo” utworzoną, więc jest 
zaselektowany (wybrany) — możemy zatem dokonać na nim dowolnej operacji. 
Naciskamy klawisze Ctrl+T (lub z menu opcja Object |'Iransform) i ustawia- 
my odpowiednie radiobuttony: 


Origin: Local 

what: ObjectśsAxes 
Rotation Z 

Z: 90 


Obrócony w ten sposób walec przesuwamy przy pomocy myszy w prawo (wzdłuż 
osi X globalnego układu odniesienia) na oknie aktywnym 1 lub 2. 


Teraz tworzymy kulę poprzez wywołanie opcji Object l Add ISphere i w 
dialog boxie wpisujemy wartości: 


Name: kula 
Radius: 20 
Horizontal Sections: 10 

Vertical Sections: 10 

Rotate Sections — ustawiony 


Kula jest zaselektowana, więc myszą przesuwamy ją w prawo, podobnie jak 


przed chwilą walec. 


a W EW O Z ZZOZ Z ZEE 
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Mamy już utworzone trzy podstawowe bryły, teraz zajmiemy się ich powiela- 
niem. 


Selektujemy walec i wywołujemy dialog box Replicate poprzez opcję Ob- 
ject | Replicate i wpisujemy wartości: 


Number: 1 
X: -80 
Y: 0 
Z: 0 


Walec zostanie skopiowany jeden raz i nowy zostanie umieszczony po lewej 
stronie sześcianu (to samo mogliśmy osiągnąć kopiując bryłę przy pomocy opcji 
Object I Copy lub klawisza Ins i następnie przemieszczając otrzymany w ten 
sposób nowy obiekt w lewo). 


Rys. 8.1.2 
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Podobnie postępujemy z kulą, a więc po zaselektowaniu jej wywołujemy opcję 
Object | Replicate i wpisujemy wartości: 


Number: 1 
X: -160 
Y: 0 
Z: 0 


Przedstawiona na rysunku na początku tego przykładu scena jest już gotowa, 
pozostaje nam jedynie odpowiednie zorientowanie obserwatora w celu obejrze- 
nia naszej konstrukcji z każdego punktu. 


Selektujemy bryłę będącą w centrum konstrukcji, czyli sześcian i wywołujemy 
opcję Scene I Observer, gdzie ustawiamy niewielki kąt patrzenia (ok. 30) oraz 
jako Target wybieramy Selected Object. 


Teraz oddalamy się od sceny przy pomocy klawisza —, aby z daleka zobaczyć 
konstrukcję i usytuowanie obserwatora na scenie. Selektujemy obserwatora 
i dowolnie przemieszczamy go po scenie; będzie on zawsze patrzył na wybrany 
przez nas sześcian. Przykładowe rysunki z różnego punktu patrzenia obserwa- 
tora są przedstawione na rys.8.1.2. 


Uwagi : 


1. W przypadku kuli i walca dobrze jest ustawić w dialog boxie Surface checkbox 
Phong, obraz otrzymany w ten sposób będzie bardziej naturalny. 


2. Gdy podczas np. przemieszczania brył obraz ulegnie częściowemu zatarciu, 
należy użyć opcji Screen | Redraw (lub nacisnąć klawisze Ctrl+W). 


3. Zamiast myszą obiekty można przemieszczać przy użyciu opcji Ob- 
ject|Transform (Translation lub Position w globalnym układzie odniesienia). 


8.2. Przykład 2 


A teraz prześledźmy budowanie sceny bardziej skomplikowanej od omówionej 
przed chwilą. Zbudujmy ludzika przedstawionego na rys.8.2.1. 


Pierwszą rzeczą, jaką utworzymy, jest tułów, który jest figurą obrotową. 
Wywołujemy więc opcję Object I Add | Figure i definiujemy krzywą o kształcie 
zbliżonym jak na rys.8.2.2. 
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Rys. 8.2.1 


Add Rotate Figure 


Name Sections 


Bottom 


lop 


DJ Rotate sections 
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Następnie tworzymy głowę, która jest kulą o parametrach: 


Name: glowa 
Radius: 25 
Horizontal Sections: 10 

Vertical Sections: 10 

Rotate Sections — ustawiony 


Tak utworzoną kulę przemieszczamy przy pomocy myszy na oknie aktywnym 
nr 2 wzdłuż osi Y globalnego układu odniesienia. W celu precyzyjnego uloko- 
wania głowy na tułowie możemy przybliżyć się do sceny, naciskając klawisz +. 


Następnie tworzymy oczy — wywołujemy dialog box do tworzenia kuli i wpisu- 
jemy parametry: 


Name: oko 

Radius: 5 
Horizontal Sections: 10 

Vertical Sections: 10 

Rotate Sections — ustawiony 


Nadal na aktywnym oknie nr 2 przemieszczamy oko wzdłuż osi Y na żądane 
miejsce, a następnie kopiujemy obiekt (klawisz Ins lub opcja Object I Copy). 
Otrzymaną w ten sposób parę oczu przemieszczamy wzdłuż osi Z na oknie 
aktywnym nr 3: jedno oko przesuwamy w lewo, drugie zaś w prawo; rzut na 
tym oknie to en face naszego ludka. 


Następnie tworzymy nos, który jest stożkiem o parametrach: 


Name: nos 
Radius: 4 
Height: 25 
Sections: 10 
Bottom — nie ustawiony 


Tak otrzymaną bryłę obracamy wokół osi X lokalnego układu odniesienia — w 
tym celu w dialog boxie Transformation (Ctr1+T lub Object! Transform) 
ustawiamy odpowiednie radiobuttony: 


Origin: Local 

what: ObjectśsAxes 
Rotation X 

X: 90 


Obrócony w ten sposób nos przemieszczamy wzdłuż osi Y na oknie 2 na żądane 
miejsce. 


FRZENNNNZENNNNNNNNNNNNNNNNNNNENNENEZEZENZEZ z OOOĆ 
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W tym uproszczonym rysunku głowy pozostaje nam tylko dodać usta, które są 
prostopadłościanem o parametrach: 


Name: usta 
width: 3 
Depth: 15 
Height: 3 


Na oknach 2, a potem 3 odpowiednio pozycjonujemy obiekt za pomocą myszy. 


Nogi (uda) bedą walcami o parametrach: 


Name: nogal 

Radius Bottom: 10 

Radius Top: 10 

Height: 60 

sections: 10 

Bottom — nie ustawiony 
Top — nie ustawiony 


Rotate Sections — ustawiony 


Otrzymaną bryłę przemieszczamy na oknie nr 2, następnie przemieszczamy 
wzdłuż osi Z w lewo na oknie nr 3, kopiujemy i przemieszczamy w prawo. 


Następnie tworzymy kolana, z których każde jest kulą o parametrach: 


Name : staw 
Radius: 10 
Horizontal Sections: 10 
Vertical Sections: 10 

Rotate Sections — ustawiony 


Z tak utworzoną kulą postępujemy podobnie jak przed chwilą z walcem, czyli 
przemieszczamy na oknie 2, potem na oknie 3 przemieszczamy w lewo, kopiu- 
jemy i nową kulę (drugie kolano) przemieszczamy w prawo. 


Następną bryłą jaką utworzymy będzie walec o parametrach: 


Name: noga2 

Radius Bottom: 8 

Radius Top: 10 

Height: 60 

Sections: 10 

Bottom — nie ustawiony 
Top — nie ustawiony 


Rotate Sections 
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— ustawiony 


181 


Postępujemy z nim analogicznie jak z poprzednimi częściami nogi. 


Z nóg pozostało nam już tylko utworzenie stóp. Każda z nich będzie prosto- 
padłościanem o parametrach: 


Name: 
width: 
Depth: 
Height: 


stopa 
30 
17 

- 


Podobnie jak inne części nogi i te przemieszczamy na oknach 2 i 3, kopiujemy 


itd. 


Przechodzimy do ostatniej części, a więc do rąk. Za kulki w ramionach mogą 
posłużyć kule będące kolanami ludzika, więc możemy je skopiować, a następnie 
jak części nóg odpowiednio przemieścić na oknach 2i3. 


Ręka będzie składała się z dwóch walców, z których pierwszy ma parametry: 


Name: 

Radius Bottom: 
Radius Top: 
Height: 
Sections: 
Bottom 

Top 

Rotate Sections 


a drugi: 


Name: 

Radius Bottom: 
Radius Top: 
Height: 
Sections: 
Bottom 

Top 

Rotate Sections 


rekal 

6 

8 
40 

10 
— nie ustawiony 
— nie ustawiony 
— ustawiony 


10 

— nie ustawiony 
— nie ustawiony 
— ustawiony 


Jako staw łokciowy posłuży nam kula o parametrach: 


Name: 
Radius: 


staw 
7 


Horizontal Sections: 10 


vertical Sections: 


Rotate Sections 


10 
— ustawiony 
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Każdy z trzech przed chwilą wymienionych obiektów przemieszczamy i kopiu- 
jemy analogicznie jak części nóg — wynika to wszystko z symetrii tworzonej 
sceny. 


Otrzymany w ten sposób ludek stoi na baczność, aby zasymulować ruch (np. 
bieg) należy odpowiednio zmienić ustawienie nóg i rąk — może być to np. 
ustawienie górnej części jednej z nóg pod kątem (obrót o kąt 45 wokół osi Z 
lokalnego układu odniesienia) i odpowiednie potem poprzemieszczanie pozo- 
stałych brył, aby tworzyły zwartą całość. 


Odpowiednie zorientowanie różnych części nóg i rąk pozostawiamy w tym 
miejscu inwencji użytkownika — jako przykład może służyć nasza scena. 


Możemy teraz spojrzeć na utworzonego ludka z każdego punktu sceny, prze- 
mieszczając po niej obserwatora. Selektujemy tułów ludka i wywołujemy dialog 
box Observer, w którym jako Target wybieramy Selected Object oraz ustalamy 
niewielki kąt patrzenia (np. 25). Dzięki takiemu ustawieniu parametrów 
obserwator zawsze będzie patrzył na ludzika (konkretnie na środek lokalnego 
układu odniesienia figury obrotowej będącej tułowiem). Teraz wybieramy 
obserwatora np. poprzez opcję Scene |SelectObserver, oddalamy się od 


Rys. 8.2.3 
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sceny i przemieszczamy po niej dowolnie obserwatora. Otrzymane w ten sposób 
przykładowe spojrzenia na ludka przedstawiamy na rys.8.2.3. 


Uwagi: 


1. We wszystkich tworzonych w tej scenie bryłach radzimy w parametrach 
powierzchni ustawić checkbox Phong. 


2. Celowo nie podawaliśmy żadnych konkretnych parametrów powierzchni, 
takich jak kolor czy połyskliwość brył, gdyż omawiany przykład był jedynie 
pokazem budowania obiektów oraz ich wzajemnej orientacji na scenie, zaś 
dopiero teraz otwierają się możliwości przed użytkownikiem, takie jak koloro- 
wanie powierzchni, ustalanie barwy lamp oraz tła. 


3. Gdy utworzona figura obrotowa (tułów) okaże się nieproporcjonalna do 
reszty ludka (np. za duża), możemy ją zmniejszyć, wykonując transformację 
Scale (w dialog boxie Transformation). 
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Dodatek A 


Skrót Specyfikacji Tag Image File Format (TIFF) 5.0 


1. Wstęp 


W celu umożliwienia użytkownikowi zapisu otrzymywanych obrazów w forma- 
cie TIFF podajemy jego skrótową specyfikację. 


Opisany format plików ułatwia wymianę obrazów zapisanych cyfrowo pomię- 
dzy różnymi aplikacjami, a także systemami komputerowymi. Pola TIFF były 
definiowane głównie z myślą o zastosowaniach DTP. Okazało się jednak, że 
format ten może być użyteczny dla innych aplikacji. TIFF jest niezależny od 
systemów operacyjnych, systemów plików, procesorów itp. Obrazy są zapisy- 
wane w plikach, które definiujemy jako ciągi bajtów, numerowanych od O do 
N. Maksymalna długość pliku wynosi 2132 bajtów. W systemach DOS, UNIX 
i OS/2 zaleca się stosowanie rozszerzeń "TIF", zaś dla komputerów Macintosh 
typem pliku jest "TIFF". 


1.1. Struktura plików TIFF 


Poszczególne pola są identyfikowane przez niepowtarzalne etykiety (ang. tags). 
Pozwala to na umieszczanie pól w pliku TIFF lub ich pomijanie, w zależności 
od potrzeby. Plik TIFF składa się z trzech głównych części: 


1) na początku krótki nagłówek (ang. file header), 
2) katalog wszystkich pól w pliku, 
3) dane dla tych pól. 


Nagłówek 
0 Porządek bajtów 
2 wersja 
4 offset do 
6 


zerowego 1LFD 


Realistyczna grafika 3D 185 


IFD 


A Ilość wejść 

A+2 Wejście do pola 0 

A+14 Wejście do pola 1 

A+26 

A+2+B*12 offset do następnego IFD 
Pole 0 

4 Etykieta (Tag) 

X+2 Typ 

X+4 Długość 


offset do danej wskazywanej 
przez to pole 


AJ 
+ 
co 


1.2. Nagłówek pliku 


Bajty 0-1: Pierwsze słowo w pliku TIFF określa porządek bajtów: 
"IT" (hex 4949) — kolejność zapisu bajtów od mniej 
znaczących do bardziej znaczących dla 16 
i 32-bitowych liczb integer. 
"MM" (hex 4D4D) — kolejność zapisu bajtów od bardziej 
znaczących do mniej znaczących dla 16 i 32-bitowych 


liczb integer. 

Bajty 2-3: Numer wersji TIFF. Wartość 42 (2A). Nie odnosi się do 
aktualnej wersji TIFF. 

Bajty 4-7: To długie słowo zawiera offset do pierwszego IFD w 


pliku. IFD (ang. Image File Directory) może być 
umieszczony w dowolnym miejscu w pliku. Offset jest 
liczbą parzystą. Termin offset oznacza przesunięcie 
danego miejsca w pliku w stosunku do początku tego 
pliku. Pierwszy bajt ma zawsze offset równy O. 
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1.3. Image File Directory 


Zawiera 2-bajtowy licznik ilości pól. Pola umieszczone są bezpośrednio za 
licznikiem. Na końcu każdego IFD jest offset do następnego IFD (lub O, jeśli 
takiego nie ma). Każde pole jest ciągiem 12 bajtów: 


Bajty 0-1: Zawierają nagłówek pola. 

Bajty 2-3: Typ pola. 

Bajty 3-4: Długość pola. 

Bajty 8-11: Offset — parzysta liczba określająca położenie danych 


dla tego pola. 


Pola muszą być umieszczone w kolejności rosnących numerów etykiet. Jeśli 
długość danej jest mniejsza bądź równa 4 bajty, to dla zaoszczędzenia miejsca 
jest ona umieszczana zamiast offsetu. Wielkość danej jest zdeterminowana 
przez jej typ: 


l — BYTE — 8 bitowa liczba całkowita bez znaku. 

2 — ASCII — ciąg bajtów przechowujący tekst, zakończony 0. 

3 — SHORT — 16 bitowa liczba całkowita bez znaku. 

4 — LONG — 32 bitowa liczba całkowita bez znaku. 

5 — RATIONAL  —dwie liczby LONG, pierwsza oznacza licznik, druga 


mianownik ułamka. 


Długość danej typu ASCII jest długością tekstu + 1 (końcowe 0). Dla niektórych 
pól istnieje możliwość wybrania jednego z dwóch typów SHORT lub LONG. W 
każdym pliku TIFF może być więcej niż jeden IFD. Każdy IFD definiuje 
"podplik"'. Można w ten sposób tworzyć "podobrazy" umieszczone na części 
obrazu głównego. 


1.4. Definicje 


Pojęcia używane w opisie TIFF'a można zdefiniować w następujący sposób: 
Obraz (ang. image) jest prostokątną tablicą piksli (punktów obrazu), z których 
każdy jest opisany przez jedną lub więcej próbek (ang. samples). W obrazie 
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monochromatycznym mamy jedną próbkę na piksel; dane obrazu kolorowego 
zawierają trzy próbki na piksel (RGB). 


Klasy obrazów w formacie TIFF 


. Wstęp 


TIFF był projektowany, aby ułatwić życie sprzedawcom skanerów i twórcom 
programów DTP, a także użytkownikom tych produktów, poprzez zmniejszenie 
ilości formatów zapisu obrazu. Oczekiwano, że będzie nim zainteresowana 
niewielka grupa osób. Rozmiar sukcesu, jaki odniósł TIFF, był zaskoczeniem 
dla jego twórców. Pojawił się jednak pewien problem. TIFF miał być formatem 
elastycznym. Elastyczność tę osiągnięto kosztem prostoty. W związku z tym, 
że trudno jest napisać program obsługujący wszystkie opcje (prawdopodobnie 
jeszcze takiego nie stworzono), twórcy TIFF'a doszli do wniosku, że należy 
uszczuplić nieco jego możliwości, tak aby otrzymać standard możliwie prosty 
do zaimplementowania. Wybrano tylko te opcje, które były najczęściej używa- 
ne. Była to oczywiście decyzja arbitralna, jednak od tej pory łatwiej jest 
zapisywać obrazy, tak by mogły być odczytane przez różnorodne aplikacje. 


Zdefiniowano cztery klasy dla różnych typów obrazów: 


1) klasa B — obrazy 1 bitowe (ang. bitlevel images) 

2) klasa G — obrazy o szarej skali jasności (ang. grayscale images) 
3) klasa P — obrazy kolorowe (ang. palette color images) 

4) klasa R — obrazy RGB (ang. RGB full color images) 


Dla uproszczenia w dalszej części opisu będziemy używać skrótów: "TIFF B, 
"TIFF G", "TIFF P", "TIFF R". Mówiąc o wszystkich klasach będziemy pisać 
"TIFF X". 


2.1. Wymagania ogólne dotyczące TIFF X 


Jeśli w opisie wystąpią opcje, to znaczy, że do zapisu obrazu można użyć 
dowolnej z nich, jednak programy odczytujące muszą obsługiwać wszystkie z 
nich. Szczególną uwagę należy zwrócić na zalecenia. Istnieje bowiem możli- 
wość, że w przyszłości klasy TIFF zostaną tak zdefiniowane, że będą zawierać 
wyłącznie opcje obecnie zalecane. 
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2.2. Pola o wartości standardowej 


Nie muszą być zapisywane, jeśli ta wartość jest odpowiednia. Program odczy- 
tujący obraz powinien oczekiwać na dowolną wartość w takim polu (niekonie- 
cznie standardową). 


2.3. Porządek bajtów "MM" i"ll" 


Odczytujący obraz musi być przygotowany do prawidłowego obsłużenia obu 
porządków. Zadziwiająco często obrazy przekraczają granice między PC a 
komputerami zawierającymi procesory firmy Motorola. Można by wymusić na 
twórcach oprogramowania stosowanie określonej kolejności zapisu bajtów, 
jednak zaczną pojawiać się problemy z dużymi obrazami RGB. Ponadto odwra- 
canie kolejności bajtów w czasie wczytywania obrazu skanerem może spowo- 
dować zwolnienie procesu czytania lub nawet zatrzymanie jego mechanizmu. 


2.4. Wielokrotne podpliki (ang. subfiles) 


Program odczytujący musi być przygotowany na to, że plik może zawierać 
więcej niż jeden obraz, chociaż potrzebny jest tylko pierwszy z podobrazów. 
Zapisując obraz lub kilka obrazów w jednym pliku w formacie TIFF, należy 
pamiętać, aby koniecznie wstawić cztery zerowe bajty po ostatnim IFD (ozna- 
cza to, że ten podobraz był ostatni w pliku). Jeśli zapisujemy kilka podobrazów, 
to pierwszy w pliku musi być obrazem o pełnej rozdzielczości. Pozostałe 
podobrazy mogą być zapisane w dowolnej kolejności. 


2.5. Pola 


Poniżej opisujemy pola, które muszą być umieszczone w każdym pliku TIFF 
X. Opis każdego pola zaczyna się od jego nazwy. Tag to liczba spełniająca rolę 
etykiety (w nawiasach jest podana wartość hexagonalna), a N oznacza ilość 
spodziewanych wartości. 


Typy danych: 
BYTE — 8-bitowa liczba całkowita bez znaku. 
SHORT — 16-bitowa liczba całkowita bez znaku. 
LONG — 32-bitowa liczba całkowita bez znaku. 
RATIONAL — dwie liczby LONG, pierwsza oznacza licznik, druga 


mianownik ułamka. 


Realistyczna grafika 3D 189 


2.5.1. NewSubfileType 


Tag = 254 (FE) 
Type = LONG (zalecane) 
N =l 


Rodzaj danych w podpliku. Pole to składa się z 32 
bitów. Nie używane bity to 0. Bit O to najmniej 
znaczący bit. 


Aktualnie zdefiniowane wartości: 


Bit0 = 1, jeśli obraz jest wersją o zredukowanej rozdzielczości 
innego obrazu tego samego pliku TIFF (w przeciwnym 
wypadku 0). 
Bitl = 1, jeśli obraz jest stroną z wielostronicowego obrazu. 
Bit2 = 1, jeśli obraz definiuje maskę przezroczystości dla innego 
obrazu w tym samym pliku. 
Przykład: 


Cztery obrazy w jednym pliku TIFF: obraz o pełnej rozdzielczości, obraz o 
zredukowanej rozdzielczości i dwie maski przezroczystości dla tych obrazów. 
Każdy z tych obrazów będzie miał inną wartość pola NewSubfileType. 


Wartość standardowa: O. 


2.5.2. ImageWidth 


Tag = 256 (100) 

Type = SHORT lub LONG (zalecany jest LONG — coraz 
większe obrazy) 

N =l 


Szerokość (liczba kolumn) obrazu w pikslach. 


2.5.3. ImageLength 


Tag = 257 (101) 

Type = SHORT lub LONG (zalecany jest LONG — coraz 
większe obrazy) 

N =l 
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Długość (wysokość) obrazu w pikslach. 


2.5.4. RowsPerStrip 


Tag = 278 (116) 
Type = SHORT lub LONG (zalecany LONG) 
N = l] Wartość pola to liczba linii w pasku. Obrazy są 


podzielone na paski, aby przyspieszyć dostęp do 
danych. Ilość pasków w obrazie można obliczyć w 
następujący sposób: 

StripsPerlmage = Int((RowsPerStrip+ImageLength -1) 
/ RowsPerStrip) 

Standardowo wartość tego pola jest równa 2132-1, co 
oznacza nieskończoność. Nie jest zalecane zapisywanie 
całego obrazu w jednym pasku. Liczbę RowsPerStrip 
dobiera się tak, aby pasek zawierał ok. 8000 próbek 
(nie bajtów, nie piksli). Wartość tę stosuje się nawet 
dla nie skompresowanych obrazów. 

Zob. ImageLength, StripOffsets, StripByteCounts. 


2.5.5. StripOffsets 


Tag =273 (111) 

Type = SHORT lub LONG (zalecane jest używanie zawsze 
LONG) 

N = StripsPerIlmage, gdy PlanarConfiguration = 1 


= SamplesPerPixel*StripsPerilmage, gdy 
PlanarConfiguration = 2 

Dla każdego paska offset od początku pliku do 
początku tego paska. Położenie danego paska jest 
niezależne od położenia innych pasków. 

Zob. StripByteCounts, RowsPerStrip. 


2.5.6. StripByteCounts 


Tag = 279 (117) 

Type = SHORT lub LONG (zalecany SHORT — paski nie są 
zbyt duże) 

N = StripsPerImage, gdy PlanarConfiguration = 1 


= SamplesPerPixel*StripsPerlmage, gdy 
PlanarConfiguration = 2 
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2.5.7. XResolution 


2.5.8. YResolution 
Tag 
Type 
N 


Dla każdego paska liczba bajtów w tym pasku. Pole to 
jest w zasadzie bezużyteczne, ale może ułatwić 
wczytywanie skompresowanego obrazu. 

Zob. StripOffsets, RowsPerStrip. 


= 282 (11A) 
= RATIONAL 


=|1 Ilość piksli/ResolutionUnit w kierunku X. 
Rozmiary drukowanego obrazu zależą od wartości tego 
pola (a także od YResolution). 

Zob. ResolutionUnit 


= 283 (11B) 
= RATIONAL 


=l 
Ilość piksl;i/ResolutionUnit w kierunku Y. 
Zob. ResolutionUnit 


2.5.9. ResolutionUnit 


3. Klasy TIFF X 


3.1. Klasa TIFF B 


= 296 (128) 
= SHORT 


=l 

Pole używane łącznie z XResolution i YResolution. 
Wartość: 

=] — nie ma jednostki miary 

=2 — cale 

=3 — centymetry 


Pola, które muszą być umieszczone w pliku TIFF, poza opisanymi w pkt.2. 
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3.1.1. SamplesPerPixel 


Tag = 277 (115) 

Type = SHORT 

N =l 
Ilość próbek na piksel. Dla obrazów klasy B wartość 
tego pola = 1. 


Zob. BitsPerSample 


3.1.2. BitsPerSample 


Tag = 258 (102) 

Type = SHORT 

N = SamplesPerPixel 
Ilość bitów opisujących jeden piksel. Dla TIFF B 
wartość = 1. 


Zob. SamplesPerPixel. 


3.1.3. Compression 


Tag = 259 (103) 

Type = SHORT 

N =l 
Wartość: 


=1 — brak kompresji, ale dane są umieszczone w 
bajtach bez nie używanych bitów. Dane są zapisywane 
w tablicy typu: 
BYTE, gdy BitsPerSample <= 8 
SHORT, gdy 8 < BitsPerSample <= 16 
LONG, gdy 16 < BitsPerSample <= 32 
=2 — CCITT (zob.[10]) 
=5 — DLZW Compression (zob.[10]) 
=32773 — PackBits Compression (zob.[10]) 
Dla TIFF B zalecane jest użycie PackBits. 


3.1.4. Photometricinterpretation 


Tag = 262 (106) 
Type = SHORT 
N =l 
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0 — dla obrazów TIFF B i TIFF G. O odpowiada kolorowi 
białemu, 2ABitsPerSample-1 — czarnemu. 


1 — dla obrazów TIFF B i TIFF G. 0 odpowiada kolorowi 
czarnemu, 2ABitsPerSample-1 — białemu. 
2 — RGB. W modelu RGB kolor jest opisywany jako 


kombinacja trzech podstawowych kolorów: 
czerwonego, zielonego i niebieskiego. 


Dla każdej z trzech próbek O odpowiada minimalnej 
jasności, 20BitsPerSample-1 — maksymalnej. Oznacza 
to, że 0,0,0 odpowiada kolorowi czarnemu, a 
2550,250,255 — białemu (zakładając 8-bitowe próbki). 
Gdy PlanarConfiguration = 1, próbki są ułożone w 
kolejności: czerwony, zielony, niebieski. 

Gdy PlanarConfiguration = 2, StripOffsets dla 
poszczególnych planów są umieszczone w identycznej 
kolejności: najpierw StripOffsets dla planu próbek 
czerwonych, potem zielonych, na końcu niebieskich. 


3 — "Palette color". W tym trybie kolor jest opisywany 
przez pojedynczą próbkę, która jest używana jako 
indeks do ColorResponseCurves (zob.[10], 3.3.5). 
Indeks ten wybiera trzy składniki RGB z odpowiedniej 
tablicy. 

4 — Maska Przezroczystości (ang. Transparency Mask). 
Oznacza to, że ten obraz definiuje nieregularny obszar 
w innym obrazie należącym do tego samego pliku 
TIFF. SamplesPerPixel i BitsPerSample muszą mieć 
wartość 1. Zalecana jest kompresja PackBits. Bity o 
wartości 1 definiują wnętrze obszaru, a 0 — zewnętrze. 
Obraz odpowiadający Masce Przezroczystości ma takie 
same wartości ImageWidth i ImageLenght jak obraz 
główny. Piksle w obrazie głównym, które odpowiadają 
bitom 1 w masce, są wyświetlane na ekranie bądź 
drukowane, pozostałe (odpowiadające bitom O) nie 
będą ani drukowane, ani wyświetlane. 
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offset 
(hex) 


HEADER: 
0000 
0002 
0004 


nazwa 


ByteoOrder 
Version 
IFD pointer 


IFD dla pełnego obrazu 


0014 
0016 
0022 
002E 
003A 
0046 


0052 
005E 
006A 
0076 
0082 
008E 
00A6 


EntryCount 
NewSubfileType 
ImageWidth 
ImageLength 
Compression 
Photometric 
Interpretation 
Stripoffsets 
RowsPerStrip 
StripByteCounts 
xResolution 
YResolution 
Software 
wskaźnik do 
nast. IFD 


3.1.5 Przykładowy obraz binarny: TIFF B Image 


wartość 
(hex) 


4D4D 
002A 
00000014 


000D 

O00FE 0004 
0100 0004 
0101 0004 
0103 0003 
0106 0003 


0004 
0004 
0003 
0005 


0111 
0116 
0117 
O1IA 
011B 0005 
0131 0002 
00000000 


00000001 
00000001 
00000001 
00000001 
00000001 


000000BC 
00000001 
000000BC 
00000001 
00000001 
0000000E 


00000000 
000007D0 
00000BB8 
80050000 
00010000 


000000B6 
00000010 
000003A6 
00000696 
0000069E 
000006A6 


Obszary wskazywane przez wskaźniki w odpwiednich polach 
pierwszego IFD 


00B6 
03A6 
0696 
069E 
06A6 


Dane obrazu 


00000700 
XXXXXXXX 
XXXKXKXXXX 


Stripoffsets 


StripByteCounts 


XResolution 
YResolution 
Software 


offset0,...„,offsetl187 
count0,...,count187 
0000012C 00000001 
0000012C€ 00000001 
"PageMaker 4.0" 


skompresowane dane dla paska 10 
skompresowane dane dla paska 179 
skompresowane dane dla paska 2 


Koniec przykładu. 
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Uwagi do przykładu: 
1. Używając 16 linii na pasek, otrzymujemy łącznie 188 pasków. 
2. W przykładzie użyto opcjonalnego pola Software. 


3. Paski w naszym przykładzie są umieszczone w dowolnej kolejności. Nigdy 
nie należy spodziewać się, że za paskiem N znajduje się pasek N+1. Dostęp do 
pasków odbywa się zawsze za pośrednictwem wskaźników (StripOffsets). 


3.2. Klasa TIFF G 


Wymagana jest obecność następujących pól (poza wymienionymi w pkt.2): 


3.2.1. SamplesPerPixel = 1. SHORT 
zob.3.1.1 


3.2.2. BitsPerSample = 4, 8. SHORT 
zob.3.1.2 


3.2.3. Compression =1 lub 5. SHORT 


Zalecane jest użycie 5 ze względu na to, że kompresja DLZW jest jedną z 
efektywniejszych. 


zob.3.1.3, [10] 


3.2.4. Photometricinterpretation =0 lub 1. SHORT 
Zalecana wartość 1. 
zob.3.1.4 


3.3. Klasa TIFF P 


Wymagana jest obecność następujących pól (poza wymienionymi w pkt.2): 
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3.3.1 SamplesPerPixel = 1. SHORT 


Wartość każdego piksla jest indeksem do trzech tablic kolorów w ColorRespon- 
seCurves. 


zob.3.1.1 


3.3.2. BitsPerSample = 1,2,3,4,5,6,7,8. SHORT 
zob.3.1.2 


3.3.3. Compression =1 lub 5. SHORT 
zob.3.1.3, [10] 


3.3.4. Photometricinterpretation = 3 (Palette color). SHORT 
Zalecana wartość 1. 
zob.3.1.4 


3.3.5. ColorResponseCurves 


Tag = 301 (12D) 

Type = SHORT 

N = 2ABitsPerSample (dla próbek czerwonych) + 
24BitsPerSample (dla zielonych) + 2"BitsPerSample 
(dla niebieskich). 


Tu zdefiniowane są trzy krzywe jasności dla kolorów 
czerwonego, zielonego i niebieskiego. Każda krzywa 
ma długość 2ABitsPerSample. Każda pozycja ma 16 
bitów długości (typ SHORT). O odpowiada najmniejszej 
jasności, 65535 — największej. ColorResponseCurve 
jest tablicą transformującą wartości piksli od O do 
2ABitsPerSample-1 na odpowiednie wartości 
intensywności. 
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3.4. Klasa TIFF R 


Wymagana jest obecność następujących pól (poza wymienionymi w pkt.2): 


3.4.1. SamplesPerPixel = 3. SHORT 
Jedna próbka na każdy z trzech podstawowych kolorów. 
zob.3.1.1 


3.4.2. BitsPerSample = 8,8,8. SHORT 
zob.3.1.2 


3.4.3. Compression =1 lub 5. SHORT 
zob.3.1.3, [10] 


3.4.4. Photometricinterpretation = 2 (RGB). SHORT 
zob.3.1.4 


3.4.5. PlanarConfiguration =1 lub 2. SHORT. 


Zalecane używanie 1. 


Tag = 284 

Type = SHORT 

N =l 

1 — próbki są ułożone w sposób ciągły. Na przykład dla 
danych RGB: RGBRGBRGR... 

2 — próbki są podzielone na oddzielne plany "sample 


planes". Wartości w StripOffsets i StripByteCounts są 
rozmieszczone w dwuwymiarowej tablicy o liczbie linii 
równej SamplesPerPixel i StripsPerlmage kolumn. Na 
początku umieszczone są wszystkie kolumny dla linii 
0, potem dla linii 1 itd. Photometricinterpretation 
opisuje typ danych umieszczonych w planach. Pole to 
nie powinno być włączane do pliku TIFF, jeśli 
SamplesPerPixel = 1. Wartością standardową jest 1. 
Zob. BitsPerSample, SamplesPerPixel. 
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Dodatek B 
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Dodatek U 


W tym dodatku znajduje się lista wszystkich funkcji, do których odwoływaliśmy 
się w naszej książce. Każda funkcja jest opatrzona krótkimi komentarzem 
wyjaśniającym jej zastosowanie i sposób działania. 


f/x mmm amamamamnmam— */ 
/* * / 
/* modul listy.c * / 
/* * / 
/* mecza omomonnnomnncnozccnanmancn mmm */ 


// Grupa funkcji rezerwujących pamięć 

// dla odpowiednich struktur. 

char *new text( int size ); 

VERTEX huge *new vertex( int n ); 

EDGE huge *new edge( int n ); 

FEATURE *new feature( double ktR , double ktG , 
double ktB , double krR , 
double krG , double krB , 
double kz , int m ); 

FACE huge *new face( int n ); 

AXES *new _axes( void ); 

OBJECT *new _object( void ); 

OBSERVER *new observer( void ); 

LAMP *new lamp( void ); 

SCENE *new_scene( void ); 


// Grupa funkcji zwalniających pamięć 

// zajmowaną przez odpowiednie struktury. 
VERTEX huge *delete vertex( VERTEX huge *vert ); 
EDGE huge *delete _edge( EDGE huge *edge ); 
FEATURE *delete feature( FEATURE *fe ); 
FACE huge *delete_face( FACE huge *face ); 
void delete _vertices( VERTEX huge *vert ); 
void delete _edges( EDGE huge *edge ); 
void delete features( FEATURE *fe ); 

void delete _faces( FACE huge *face ); 
void delete _object( OBJECT *object ); 
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void delete_axes( AXES *axes ); 

void delete _lamp( LAMP *lamp ); 

void delete observer( OBSERVER * obs); 
void delete_scene( SCENE *scene ); 


// Usunięcie układu lokalnego ze sceny 
// (bez zwalniania pamięci). 
void RemoveAxesFromScene( SCENE *scene , AXES *axegS ); 


// Usunięcie lampy ze sceny (bez zwalniania pamięci). 
void RemoveLampFromScene( SCENE *scene , LAMP *lamp ); 


// Usunięcie obserwatora ze sceny (bez zwalniania 
pamięci). 
void RemoveobserverFromScene( SCENE *scene ); 


// Stworzenie wierzchołka o współrzędnych x,y,z. 
VERTEX huge *add vertex( VERTEX huge *vert, double x , 
double y , double z ); 


// Dołączenie wierzchołka do struktury 

// typu OBJECT opisującej bryłę. 

VERTEX huge *add vertex to object( VERTEX huge *vert, 
OBJECT *object , 
double x , double y , 
double z ); 


// Dołączenie wierzchołka do struktury 

// typu AXES opisującej układ lokalny. 

VERTEX huge *add vertex to axes( VERTEX huge *vert, 
AXES *axes , double x , 
double y , double z ); 


// Dołączenie krawędzi do struktury 

// typu OBJECT opisującej bryłę. 

EDGE huge *add edge to object( EDGE huge *edge , 
OBJECT *object , 
VERTEX huge *vertl , 
VERTEX huge *vert2 ); 

// Dołączenie krawędzi do struktury 

// typu AXES opisującej układ lokalny. 
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EDGE huge *add_edge to axes( EDGE huge *edge, AXES *axes, 
VERTEX huge *vertl , 
VERTEX huge *vert2 ); 


// Dołączenie struktury feature zawierającej 

// parametry powierzchni do struktury 

// typu OBJECT opisującej bryłę. 

void add_feature to object( OBJECT *object, FEATURE *fe); 

// Dołączenie ściany do struktury OBJECT 

// opisującej bryłę. 

FACE huge *add_face to object( FACE huge *face , 
OBJECT *object , 
VERTEX huge *vertl , 
VERTEX huge *vert2 , 
VERTEX huge *vert3 , 
FEATURE *fe , int phong ); 


// Funkcja zwraca najwyższy spośród numerów 

// układów lokalnych na scenie. 

int MaxNumAxes( SCENE *scene ); 

// Utworzenie struktury typu AXES 

// opisującej układ lokalny. 

AXES *add axes( char *name , int number , 
double x local scale , 
double y local scale , 
double z_ local scale ); 

// Utworzenie struktury typu LAMP 

// opisującej źródło światła. 

LAMP *add_lamp( char *name , int number , double x , 
double y , double z , double IR , 
double IG , double IB , double IZ ); 

// Utworzenie struktury typu OBSERVER opisującej 

// położenie obserwatora i ostrosłup widzenia. 

OBSERVER *add observer( double x , double y , double z , 

double modu , double modw ); 

// Utworzenie struktury typu SCENE opisującej scenę. 

SCENE *add_scene( char *name , double „double , double ); 

// Dołączenie obserwatora do sceny. 

void add observer _to scene( SCENE *sc , 

- OBSERVER *observer ); 

// Dołączenie nowego układu lokalnego 

// (wskaźnik *axes) do sceny. 

void add_axes_to_scene( SCENE *scene , AXES *axes ); 
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// Funkcja zwraca TRUE (1), jeżeli układ lokalny 
// *axes jest dołączony do sceny. 
int AreAxesOnScene( SCENE *scene , AXES *axegs ); 


// Dołączenie lampy do sceny. 
void add_lamp to scene( SCENE *scene , LAMP *lamp ); 


// Dołączenie bryły (opisywanej przez strukturę 

// typu OBJECT) do układu lokalnego. 

void add _ object to axes( SCENE *scene , AXES *axes , 
OBJECT *object ); 


// Wyznaczenie liczby wierzchołków na liście 

// wskazywanej przez *vert. 

int CalcNVertices( VERTEX huge *vert ); 

// Wyznaczenie liczby krawędzi na liście 

// wskazywanej przez *edge. 

int CalcNEdges( EDGE huge *edge ); 

// Wyznaczenie liczby ścian na liście 

// wskazywanej przez *face. 

int CalcNFaces( FACE huge *face ); 

// Skopiowanie współrzędnych wierzchołka 

// ze struktury wskazywanej przez *vertsrc 

// do struktury wskazywanej przez *vertdest. 

void CopyVertCoords( VERTEX huge *vertdest , 
VERTEX huge *vertsrc ); 

// Utworzenie kopii struktury typu AXES 

// opisującej układ lokalny. 

AXES *CopyAxes( AXES *axes old ); 

// Utworzenie kopii struktury typu OBJECT 

// opisującej bryłę. 

OBJECT *CopyObject( OBJECT *objold ); 

// Utworzenie kopii układu lokalnego 

// i znajdującej się w nim bryły. 

AXES *CopyAxesAndObject( SCENE *scene , AXES *axesold ); 

// Powielenie n razy układu lokalnego i bryły. 

AXES *ReplicateAxesAndObject( SCENE *scene , int n , 

AXES *axes , double dx , double dy , double dz ); 
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/* * / 
/% modul listy.c * / 
/* * / 
foo zonannnncnnnnnnn nocna nn—— * / 


// Przesunięcie wierzchołków znajdujących się na liście 
// wskazywanej przez *vert 
// o wektor [move_x,move_y,move z]. 
void move( VERTEX huge *vert , int where , double move_X, 
double move_y , double move_z ); 
// Skalowanie współrzędnych każdego wierzchołka 
// na liście *vert. 
void scale( VERTEX huge *vert , int where , 
double scale _x , double scale _y , 
double scale z ); 


// Obrót wokół osi x wierzchołków znajdujących się 
// na liście *vert. 
void rot_x( VERTEX huge *vert, int where, double rot x); 


// Obrót wokół osi y wierzchołków znajdujących się na 
liście *vert. 
void rot_y( VERTEX huge *vert, int where, double rot_y ); 


// Obrót wokół osi z wierzchołków znajdujących się 
// na liście *vert. 
void rot_z( VERTEX huge *vert, int where, double rot z); 


// Funkcja oblicza współrzędne wektora normalnego 
// do ściany *face. 
void vector_norm to _face( FACE huge *face ); 


// wywołanie poprzedniej funkcji dla wszystkich ścian. 
void CalculateVectorsNorm( SCENE *scene ); 


// Rzutowanie perspektywiczne wierzchołków bryły 

// wskazywanej przez *object w układzie obserwatora. 

void observer _view( OBJECT *object , OBSERVER *observer , 
int x res , int y_res , 
double xaspect , double yapsect ); 
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// Rzutowanie perspektywiczne wszystkich brył. 

void observer _view scene( SCENE *scene , int x res , 
int y_ res , double xaspect , 
double yapsect ); 


// Wyznaczenie współrzędnych wszystkich wierzchołków 
// w układzie obserwatora na podstawie ich wartości 
// w układzie globalnym. 

void CoordsInobserversorigin( SCENE *scene ); 


// wyznaczenie współrzędnych wszystkich wierzchołków 
// w układzie globalnym na podstawie ich wartości 

// w układzie obserwatora.. 

void CoordsInGlobalorigin( SCENE *scene ); 


// Skierowanie obserwatore w ten sposób, 

// aby patrzył na układ *target, lub punkt (x,y,Z) 

void LookAtTarget( OBSERVER *obs , AXES *target , 
double x , double y , double z ); 


// Funkcja oblicza dla każdego wierzchołka w każdej 
// bryle znajdującej się na scenie uśredniony wektor 
// normalny zgodnie z algorytmem opisanym 

// w rozdziale 5.5. 

void AverangeVectorsNormTovVertices( SCENE *scene ); 


// Skopiowanie współrzędnych wierzchołków 
// z lokalnych do globalnych. 
void CopyLocalCoordsToGlobal( VERTEX huge *vert ); 


// Skopiowanie współrzędnych wierzchołków 
// z globalnych lokalnych. 
void CopyGlobalCoordsToLocal( VERTEX huge *vert ); 


// Obliczenie kątów przejścia z układu lokalnego 
// do globalnego (patrz rozdział 3.4.) 
void CalculateAlfaBetaFi( AXES *axes ); 


// Funkcja przeliczająca współrzędne lokalne na globalne 
// (patrz rozdział 3.4). 
void FromLocalcoordsToGlobal( AXES *axes , 

VERTEX huge *vert ); 
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// Funkcja przeliczająca współrzędne globalne na lokalne 
// (patrz rozdział 3.4). 
void FromGlobalCcoordsToLocal( AXES *axes , 

VERTEX huge*vert ); 


// Funkcja modyfikująca kąt patrzenia obserwatora. 
void ObserverLens( OBSERVER *observer , double x , 
double y ); 


ft omamamownmomcmnmmn——— */ 
/* * / 
/* modul obiekty.c * / 
/* */ 
f/m cc—— */ 


// Utworzenie struktur opisujących prostopadłościan. 
*add cube( char *name , int number, double a , double b , 
double h , FEATURE *fe , int phong ); 


// Utworzenie struktur opisujących 

// prostokątny wycinek płaszczyzny. 

*add _plane( char *name , int number, double a, double b , 
FEATURE *fe , int phong ); 


// Utworzenie struktur opisujących stożek. 

*add_cone( char *name , int number, double h , double r, 
int bottom , int npoints , FEATURE *fe , 
int phong ); 


// Utworzenie struktur opisujących walec. 

*ąadd tube( char *name , int number , double h , 
double r _ bottom , double r _top , int bottom , 
int top , int npoints , int stagger , 
FEATURE *fe , int phong ); 


// Utworzenie struktur opisujących bryłę obrotową. 
*add figure(  VERTEX huge *curve , char *name , 
- int number , int bottom , int top , 
int npoints , int stagger , FEATURE *fe , 
int phong ); 
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// Utworzenie struktur opisujących sferę. 
*add sphere( char *name , int number , double r , 
int horizontal sections , 
int vertical sections , int stagger , 
FEATURE *fe , int phong ); 
// Utworzenie struktur opisujących pierścień. 
*add_torus( char *torus_ name , int number , double rx 
double ry , double rtorus, 
int horizontal sections , 
int vertical sections , int stagger , 
FEATURE *fe , int phong ); 


// Grupa funkcji realizujących algorytm scanline 
// opisany w rozdziale 5.2. 

LLINE huge *new line( int n ); 

// Dołączenie krawędzi do listy krawędzi LK. 


LLINE huge *add line _to_1ll( LLINE huge *l, 
FACE huge *face , 
VERTEX huge *vertl , 
VERTEX huge *vert2 , 
int x resolution , 
int y resolution ); 


// Usunięcie listy LK. 
void delete _ lines _list( void ); 


// Utworzenie listy krawędzi LK. 


int make _lines_list( SCENE *scene , int x resolution , 


int y resolution ); 


// Inicjalizacja listy krawędzi aktywnych LKA. 
int InitALL( void ); 

void AddLinesToALL( LLINE huge *lin ); 

void MakeFacesUnvisibleALL( void ); 
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void RemoveLinesFromaLL( int row ); 

void AddCALL( void ); 

void SortALL( void ); 

int InitAFL( void ); 

void Reverse( FACE huge *face ); 

double Dist( FACE huge *face , double x , double y ); 
FACE huge *FindvVisibleFace( double xd , int yscr ); 
void FaceVisibilityAFL( FACE huge *face ); 


// Główna funkcja. 

int ScanLine( SCENE *scene, int x resolution , 
int y resolution , double xaspect , 
double yapsect ); 


f* mna ———— */ 
/* * / 
/* modul rysuj.c * / 
/* * / 
f/m omazonomamnnnamnamanaemanneoaana=— */ 


// Rysowanie rzutu równoległego. 
void Drawobject( OBJECT *object ) 


// Rysowanie rzutu perspektywicznego. 
void DrawoObserverViewobject( OBJECT *object ) 


// Rysowanie rzutów wszystkich brył znajdujacych sie 

// na scenie. 

void DrawScene( SCENE *scene ) 

void LogicalToScreen(VERTEX huge *vert, double xScreen, 
double yScreen, double Xmin, 
double Ymin, double Xmax, 
double Ymax ) 

void Paletakolorow( void ) 
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Opis dyskietki 


RAY3D.EXE 

SCN/ allfigs.secn 
kielich.secn 
stolik.scn 
czastka.scn 
lampa.scn 
przykll.scn 
przykl2.scn 
szachy.sen 
waza.sen 


source/  scan.c 
rysuj.c 
listy.c 
mat.c 
obiekty.c 
ray.h 
progl.c 
prog2.c 
prog3j.c 
prog4.c 
prog5.c 
progó.c 
prog7.c 
prog8.c 
svga256.bgi 
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Wydanie pierwsze 
TWORZENIE 
REALISTYCZNEJ 
GRAFIKI 3D 


Zawiera: 


Teoretyczne podstawy grafiki trójwymiarowej 
i przykłady ich praktycznego wykorzystania. 


Tworzenie realistycznych obrazów trójwymiarowych 
(takich, jaki jest na naszej okładce). 


Algorytmy (do własnego wykorzystania), napisane 
w języku C, umożliwiające generowanie własnych 
obiektów trójwymiarowych (3D). 


Dyskietkę z programem RAY3D.EXE (dla WINDOWS) 
pozwalającym na tworzenie trójwymiarowych brył, ich 
transformacje, cieniowanie powierzchni i oświetlanie 
różnokolorowymi lampami. Pozwala także tworzyć 
trójwymiarowe, realistyczne sceny. 


EF] Wiele praktycznych przykładów, a na dyskietce - 
bibliotekę procedur napisanych w standardowym 
języku C, co umożliwi uruchamianie przykładów 
na różnych typach komputerów. 
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